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Abstract 

We present the results from the heavy quarks and quarkonia working group. 
This report gives benchmark heavy quark and quarkonium cross sections for 
pp and pA collisions at the LHC against which the AA rates can be compared 
in the study of the quark-gluon plasma. We also provide an assessment of 
the theoretical uncertainties in these benchmarks. We then discuss some of 
the cold matter effects on quarkonia production, including nuclear absorption, 
scattering by produced hadrons, and energy loss in the medium. Hot matter ef- 
fects that could reduce the observed quarkonium rates such as color screening 
and thermal activation are then discussed. Possible quarkonium enhancement 
through coalescence of uncorrected heavy quarks and antiquarks is also de- 
scribed. Finally, we discuss the capabilities of the LHC detectors to measure 
heavy quarks and quarkonia as well as the Monte Carlo generators used in the 
data analysis. 
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1. INTRODUCTION ^ 



Heavy quarks are a sensitive probe of the collision dynamics at both short and long timescales. On one 
hand, perturbative heavy-quark production takes place at r oc 1/ mq . On the other hand, heavy quarks 
decay weakly so that their lifetime is greater than that of the medium created in heavy ion collisions. In 
addition, quarkonium states, bound QQ pairs, have binding energies of the order of a few hundred MeV, 
comparable to the plasma screening mass. In a quark-gluon plasma, interactions with hard gluons can 
overcome this threshold, leading to a large probability for quarkonium breakup. 

However, for heavy quarks and quarkonium to be effective plasma probes, the baseline production 
cross sections should be well established since one would ideally like to normalize the AA rates to pp 
collisions. The following two sections of this chapter specifically address the baseline rates for open 
heavy flavors and quarkonia in turn. We emphasize that the baseline rates in pp, pA and AA collisions 
are large enough for high statistics studies of all quarkonium and heavy quark states. Thus a complete 
physics program can be carried out. 

The heavy flavor cross sections will be calculated to next-to-leading order (NLO) in perturbative 
QCD. No data are available on the cc total cross sections at collider energies to better fix production pa- 
rameters such as the mass, nic, and the factorization and renormalization scales, ^ip and respectively, 
at high energies. However, it should be possible to reliably interpolate from pp measurements at 14 TeV 
if no pp run at 5.5 TeV is available, as discussed in section ID 

There is still some uncertainty in the quarkonium production mechanism. The quarkonium base- 
line rates are calculated using the color evaporation model (CEM) at NLO for the total cross sections. The 
Pt distributions from the CEM are compared to those calculated in nonrelativistic QCD (NRQCD). The 
comover enhancement scenario for quarkonium production is also introduced although no predictions 
for the rates are given. 

Once the pp baselines are discussed, the effects of modifications of the parton distribution func- 
tions in nuclei, referred to here as shadowing, are discussed in sections |3 and |3 At small momentum 
fractions, the nuclear gluon distribution is expected to be substantially reduced relative to that in the 
nucleon. We show how both the total rates and the rapidity and pr distributions could be affected by 
these modifications. Although the nuclear gluon distribution is not well known, it can be measured in 
pA and Ap relative to pp interactions at the LHC. The AA baseline rates without any other nuclear ef- 
fects can then be more reliably predicted. Final state effects on the quarkonium rates are discussed in 
sections I41I6!! 

It is well known that other effects in cold matter can change the expected quarkonium rates. Nu- 
clear absorption and secondary scattering with produced particles (comover scattering) are both effects 
that can cause quarkonium states to break up. In section]^ the absorption cross section is postulated to 
increase with energy. Its value is extrapolated to LHC energies and the effects on the expected J/V' and 
T rates in Pb-i-Pb and Ar-nAr collisions are discussed. The cross sections for quarkonium interactions 
with comovers have been studied in a number of models. Many estimates of the comover cross section 
exist. A discussion of one of the most recent is also presented in section |^ The last topic discussed in 
this section is energy loss in the medium. It is shown that for angles less than niQ / E, soft gluon radiation 
is suppressed, reducing the energy loss of c and h quarks relative to massless partons. The effects on the 
D /tt ratio are considered for both cold matter and 'hot' plasma. 

Quarkonium suppression in a quark-gluon plasma was first predicted by Matsui and Satz fV\. No 
bound states should exist at temperatures T > Td when the screening radius, 1/ ^^{T), is smaller than 
the size of the bound state even though the bound state may still exist above the critical temperature 
for deconfinement, i.e. Tc < T^. Early estimates of the dissociation of quarkonium by color screening 
suggested that the charmonium states above the J/ijj would break up near Tc, as do hadrons made up of 
light quarks, while the J / ip would need a somewhat higher temperature to break up. Similar results were 

'Author: R. Vogt. 
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predicted for the T family except that the tightly bound T{1S) appeared to break up only at temperatures 
considerably above Tc. Later, it was realized that changes in the gluon momentum distributions near 
Tc as well as thermal activation of the quarkonium states can contribute significantly to quarkonium 
suppression, leading to breakup for T < To- More recent lattice calculations suggest that screening may 
not be an efficient mechanism near Tc so that gluon dissociation may dominate quarkonium suppression. 
In this case, there is no specific transition temperature. These results are summarized in section l6l 

The most natural prediction is then that the per nucleon quarkonium rates in AA collisions would 
be substantially reduced relative to those in pp and pA collisions. However, if the large number of 
QQ pairs produced at LHC energies is taken into account, the initially produced uncorrelated heavy 
quark pairs could provide an important source of final-state quarkonium, particularly charmonium. This 
coalescence mechanism could lead to quarkonium enhancement at the LHC instead of suppression. Some 
predictions for the LHC are given in section |3 The RHIC J/V' data will place important constraints on 
this type of mechanism even though the number of uncorrelated cc pairs is much smaller than at the 
LHC. Since the bb rate is high enough for uncorrelated bb production in heavy ion collisions at the LHC, 
a smaller effect on T production may also be expected. 

Both the baseline distributions and the in-medium effects introduced here need to be simulated for 
data analysis. This is typically done with a standard set of Monte Carlo tools. In section^ heavy quark 
production with the pp generators Pythia and Herwig are compared to LO and NLO calculations. 
Calculations in AA generators are also briefly discussed. 

Finally, we turn to the experimental capabilities for measurements of these states in heavy ion 
collisions at the LHC. The ALICE and CMS collaborations have made detailed studies of quarkonium 
detection, shown in sections l8l and l9l respectivelv. ALICE may be able to reconstruct heavy flavor 
mesons. This effort is also discussed in section!^ CMS is less likely to reconstruct heavy flavor mesons 
through their decay channels but can focus on their contributions to the dilepton continuum, as discussed 
in section!^ Since the ATLAS heavy ion effort is more recent, only a short summary of their potential 
capabilities is included in section flOl 

It is clear that systematic efforts in this area will provide a rich array of data. These data will pro- 
vide important information about the gluon distribution in the nucleus in cold matter and the properties 
of the medium in the hot matter produced in AA collisions. 
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2. BENCHMARK HEAVY QUARK CROSS SECTIONS ^ 

In this Section we collect some benchmark results for the rates and kinematics features of charm and 
bottom quark production. We shall consider pp, pPb, Pbp and Pb+Pb collisions. We shall always assume 
a beam energy of 2.75 TeV per nucleon for Pb and 7 TeV for protons. We shall also study pp collisions 
with 2.75 TeV proton beams. 

The main goal of this Section is to establish to which extent pp and pA measurements can be used 
to infer the expected Pb+Pb rates in the absence of high-density effects specific to the nuclear-nuclear 
environment. We shall show that, while the prediction of absolute rates at the LHC energies is affected 
by large theoretical uncertainties, the extrapolation of the rates from one energy to another, or from pp to 
pPb, is under much better theoretical control. Our results stress the importance of pA and lower energy 
pp control runs. 

Our calculations were performed in the framework of next- to-leading order, NLO, QCD ||2||3l|4| 
13. We expect that NLO QCD properly describes the main features of the production mechanism in pp 
and pA collisions (see e.g. Refs. IHEl) and, in particular, can well describe the extrapolation of these 
features from one energy or beam type to others. We note that data are available for bottom and, most 
recently, charm production in pp collisions at energies up to almost 2 TeV. The absolute rates predicted 
by NLO QCD are affected by a theoretical uncertainty of approximately a factor of 2. Large corrections 
also arise from the inclusion of resummation contributions and from the nonperturbative fragmentation 
of the bare heavy quarks into hadrons. Recent studies {W W TO"] indicate nevertheless that, once all 
known and calculable effects (higher-order logarithmically-enhanced emissions, fragmentation effects, 
and initial-state radiation) which complement the plain NLO calculation are included, the picture which 
emerges form the comparison with the Tevatron data is very satisfactory. Indeed, comparison of the 
theory with charm production data 1 1 1 12 1 is very good even though the small charm quark mass might 
have suggested that reliable rate estimates were not possible. While including effects beyond NLO is 
crucial for reliable absolute predictions, accounting for them here would not affect the thrust of our 
analysis. We show that, once we can normalize the production properties against data from pp and pPh 
collisions, the extrapolation to Pb-i-Pb is well constrained. It should thus be possible to safely isolate 
quark-gluon plasma effects. 

The theoretical systematics we shall explore include: 

• the heavy quark masses; 

• the choice of factorization and renormalization scales (/ii? and /ijj); 

• the choice of parton distribution functions (PDFs). 

We explore the uncertainties due to the PDF choice at the level of the proton parton densities, using 
the MRS99 |13l set as a default, and the more recent sets MRST2001 [14| and CTEQ6M [15| for 
comparisons. We use the EKS98 1.16 . , 17.1 parameterization of the modifications of the parton distribution 
functions in nuclei. Although originally fitted with GRV LO proton PDFs, as argued in Ref. [.17 1 , these 
ratios are to a good extent independent of the PDF choice. 

In the case of charm production we consider the following mass and scale input parameter sets V: 



V 


PDF 


rric (GeV) 






I 


MRS99 


1.5 


2 


1 


II 


MRS99 


1.2 


2 


1 


III 


MRS99 


1.8 


2 


1 


IV 


MRS99 


1.5 


2 


2 


V 


MRST2001 


1.5 


2 


1 


VI 


CTEQ6M 


1.5 


2 


1 



^Authors: S. Frixione, M. L. Mangano. 
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log(xi) log(xi) 



Fig. 1: Parton-a;i distributions for cliarm production in tlie full pseudorapidity range (left) and in the central region (right). The 
cross section per nucleon is given. 



where 

/io = ^{pl^,+pl^^)/2 + ml . (1) 
In the case of bottom quarks, we have instead: 



V 


PDF 


nib (GeV) 


/"f//"o 


/^i?//"0 


I 


MRS99 


4.75 


1 


1 


II 


MRS99 


4.5 


1 


1 


III 


MRS99 


5 


1 


1 


IV 


MRS99 


4.75 


2 


0.5 


V 


MRS99 


4.75 


0.5 


2 


VI 


MRST2001 


4.75 


1 


1 


VII 


CTEQ6M 


4.75 


1 


1 



with jjLQ defined as above but now h replaces c. 
2.1 X coverage in charm production 

We first study the range of x values that contribute to charm production. Tliis question is relevant for 
understanding whether the current parameterizations of nuclear densities allow safe extrapolations. We 
also want to establish whether the constraints set by charm measurements in pPh collisions cover x 
ranges which are relevant for Pb+Pb. As a reference, we shall use parameter set "P = I here and limit 
ourselves to LO predictions. 

We shall consider the full pseudorapidity range as well as limited pseudorapidity regions defined 
by the acceptance coverage typical of ALICE. In particular, we consider the central region, \rj\ < 0.9, 
and the forward region, 2.5 < < 4. 

The left (right) plot in Fig. ^ shows the differential rate distribution for cc production over the 
full (central) pseudorapidity range as a function of xi, the momentum fraction of partons in the beam 
traveling in the positive rj direction. The X2 distribution is identical for pp and Pb+Pb collisions, while 
the xi and X2 distributions are interchanged for pPb and Pbp collisions. The same distributions, obtained 
for charm pairs in the forward pseudorapidity regions are given in Fig. |2l where the two peak structures 
refer to the X2 (left peak) and xi (right peak) variables. In the case of central production, we note that 
the bulk of the cross section comes from x values peaked at 10^^. The gluon density of the proton in 
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Fig. 2: Parton-a; distributions for charm production in the forward pseudorapidity range 2.5 < r; < 4. The cross section per 
nucleon is given. 



this region is well known, thanks to the HERA data. No data are, however, available for the nuclear 
densities in this x region, at least not in a range of relevant for charm production. We see that the x 
distributions for Pbp and Pb+Pb collisions are very similar in shape. This suggests that a determination 
of nuclear corrections to the nucleon PDFs extracted from a Pbp run would be sufficient to properly 
predict the Pb+Pb behaviour. 

In the case of forward production, the x ranges probed inside the two beams are clearly asym- 
metric. The selected x ranges have almost no overlap with the domains probed in central cc production. 
In particular, the range relevant for the beam traveling in the negative r] direction peaks below 10~^, a 
region where data for ~ 10 GeV^ are not now available. Assuming a reliable extrapolation of the 
HERA data to this domain for the proton beam, the determination of the nuclear corrections for Pb will 
therefore require a pPh run. A Pbp run would probably have lower priority since the peak of the xi 
distribution is at xi > 10^^ where current data are more reliable. 

2.2 Total cross sections 

We study here the predictions for total cross sections, starting with charm production. Table [2 gives 
results for pp collisions at both = 14 TeV and 5.5 TeV. The rates vary over a large range and the K 
factors, quantifying the size of NLO corrections, are very large. As anticipated previously, these factors 
make absolute rate estimates quite uncertain. Nevertheless, the predicted extrapolation from 14 to 5.5 
TeV is independent of the chosen parameter combination at the level of few percent, as shown in the 
last two columns of Tabled Similar results and conclusions are obtained for bottom quarks, detailed in 
Table El 

We next consider the extrapolation of cross sections within the central. Table |3l and forward. 
Table |4j acceptance regions. We also consider the extrapolation of the pPb, Pbp and pp rates to Pb+Pb. 
As in the case of the total rates, the extrapolations appear to have a limited dependence on the parameter 
set, contrary to the large variations of the absolute rates. We do not exphcitly quote the PDF dependence 
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Table 1: Production cross sections for charm pairs in pp collisions at 14 TeV and 5.5 TeV. The penultimate column gives the 
ratios of cross sections at the two energies for the various parameter sets. The last column gives the ratios, normalized to the 
ratio obtained with V = 



V 


o"i4 (mb) 


K = ctnlo/o'lo 


0-5.5 (mb) 


R{V) = (J^.^/cTii 


R{V)/R{V = I) 


I 


10.4 


1.7 


5.4 


0.52 


1 


II 


16.7 


1.7 


9.2 


0.55 


1.06 


III 


6.8 


1.7 


3.4 


0.50 


0.96 


IV 


7.3 


2.1 


3.7 


0.51 


0.98 


V 


8.57 


1.8 


4.2 


0.49 


0.94 


VI 


10.6 


1.8 


5.3 


0.50 


0.96 



Table 2: Production cross sections for bottom pairs in pp collisions at 14 TeV and 5.5 TeV. The penultimate column gives the 
ratios of cross sections at the two energies for the various parameter sets. The last column gives the ratios, normalized to the 
ratio obtained with V = 



V 


cji4 (mb) 


^ = 0-NLO/0"LO 


0-5.5 (mb) 


^(^) = 0-5.5/0-14 


R{V)/R{V = I) 


I 


0.43 


2.3 


0.17 


0.40 


1 


II 


0.51 


2.4 


0.20 


0.39 


0.98 


III 


0.37 


2.3 


0.15 


0.41 


1.03 


IV 


0.66 


1.4 


0.26 


0.39 


0.98 


V 


0.20 


3.2 


0.088 


0.44 


1.1 


VI 


0.40 


2.4 


0.17 


0.43 


1.08 


VII 


0.45 


2.4 


0.18 


0.40 


1 



since it is very small. 

While in the case of central production there is clearly no difference between Pbp and pPb, the 
forward acceptance has an intrinsic asymmetry. In is interesting to note that the extrapolation of the 
forward rate appears to be more stable when using a Pbp normalization rather than a pPh one. 

We find qualitatively similar results for bottom production. 
2.3 Inclusive pr spectra 

Here we present more detailed studies of the reliability of extrapolating the charm transverse momentum 
spectra. Figure 15] shows the ratios of the charm and bottom quark pT spectra in pp collisions at 14 and 
5.5 TeV. Rather than the differential spectra, we present the rates integrated above a given px threshold, 
PTmin- The curvcs correspond to different parameter sets. In spite of the strong sensitivity of the ratios 
to the pseudorapidity range and to pTmin> we see once again that the dependence on the input parameter 

Table 3: Single inclusive charm production cross sections per nucleon in the central region. The pp cross section is calculated 
at = 5.5 TeV. 



V 


0-pp (mb) 


o-pbp (mb) 


o-pbPb (mb) 


O-pbPb /o-pp 


O-pbPb/o-pbp 


I 


0.88 


1.00 


0.59 


0.67 


0.59 


II 


1.5 


1.6 


0.98 


0.65 


0.61 


III 


0.55 


0.64 


0.38 


0.69 


0.59 


IV 


0.58 


0.66 


0.40 


0.69 


0.61 
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Table 4: Single inclusive charm production cross sections per nucleon in the forward region. The pp cross section is calculated 
atVS ^ 5.5 TeV. 



V 


0-pp (mb) 


o-pbp (mb) 


cjpPb (mb) 


o-pbPb (mb) 


fPbPb/fpp 


CPbPb/c^Pbp 


0'PbPb/<7pPb 


I 


0.75 


0.92 


0.79 


0.50 


0.67 


0.54 


0.63 


II 


1.3 


1.5 


1.3 


0.83 


0.64 


0.55 


0.64 


III 


0.47 


0.60 


0.51 


0.33 


0.70 


0.55 


0.65 


IV 


0.51 


0.65 


0.56 


0.36 


0.71 


0.55 


0.64 




set is very small, ensuring a rather safe extrapolation. 

In Fig. 0] we present similar ratios for charm production, comparing the spectra in 5.5 TeV pp 
collisions with pPb and Fbp interactions. We consider here only predictions obtained with rric = 1.5 
GeV but different values of the renormalization scale, V = 1 and IV. 

In Fig.|5]we consider the ratio of the spectra in pPb and Pbp collisions to Pb+Pb. We note that, 
while the pPb/PbPb ratio shows a strong pT and t] dependence, the Pbp/PbPb ratios are quite similar up 
to an overall normalization factor. This is because production in the forward proton region is enhanced 
by the larger energy of the proton relative to the nucleon energy in the Pb nucleus. When the charm 
is instead detected in the forward region, in the direction of the Pb nucleus, the spectrum is much less 
sensitive to whether the "target" is a proton or another Pb beam. In this case, therefore, a Pbp run 
provides a better normalization benchmark compared to a pPb run. This is opposite the conclusion 
reached in subsection 12. ll for the determination of the small x gluon density in Pb. 

2.4 A(j) azimuthal correlations 

Azimuthal correlations between heavy quark pairs could provide useful information on the mechanisms 
of heavy quark production in central nuclear-nuclear collisions. A large fraction of these correlations 
could be washed out by the fact that the average number of charm pairs produced in Pb+Pb collisions 
is large and the detection of two charm quarks in the same event is no guarantee that they both orig- 
inate from the same nucleon-nucleon interaction. However, the probability of tagging correlated pairs 
increases at large px since these rates are smaller. 

The plots in Fig. |6l represent the distribution of the azimuthal difference, Ac^, between the Q and 
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Fig. 6: Charm azimuthal correlations in pp collisions. The solid lines correspond to \/5 = 5.5 TeV, the dashed lines to 
^/S — 14 TeV, rescaled by the amount indicated. 



Q, defined as: 

AcP = \(j)Q - , (2) 

in pp collisions at 14 and 5.5 TeV. The upper plots refer to distributions for pairs in the central and 
forward regions without any pr cut. The lower plots require a minimum pr of 5 GeV for both the c 
and c. Note both the very strong suppression of back-to-back production over all px, indicating very 
strong Sudakov suppression effects, and the enhancement of production in the same hemisphere at high 
Pt, a result of the dominant contribution from final-state gluon splitting into cc. Also note the similar 
shape predicted for 14 and 5.5 TeV collisions after rescaling the overall normalization as well as the mild 
dependence on the input parameters. 

The results for bottom pair production are given in Fig. In this case, the dominant contributions 
are in the back-to-back region, both without and with the px cut. The distributions are similar when Pb 
beams are considered. 

2.5 Conclusions 

We summarize here the findings of this study. We find that, in spite of the large uncertainties involved 
in the absolute predictions for charm and bottom quarks at the energies of relevance for the LHC, the 
correlation between the rates expected at different energies or with different beams is very strong. Using 
data extracted in pp collisions at 14 TeV it is possible to predict the rates expected in pp collisions at 
5.5 TeV to within few percent. A comparison of these predicted rates with pPb and Php measurements 
will therefore allow a solid extraction of the nuclear modifications of the parton PDFs in Pb. In spite of 
the different energy configurations for the pPb and the Pb-i-Pb runs, it is possible to extrapolate to Pb-i-Pb 
collisions in the absence of strong medium-dependent modifications. Differences with respect to these 
predictions could then used to infer properties of the production and propagation of heavy quarks in the 
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Fig. 7: Bottom azimuthal correlations in pp collisions. The solid lines correspond to \/S = 5.5 TeV, the dashed lines to 
\/S = 14 TeV, rescaled by the amount indicated. 



dense medium resulting from high energy Pb+Pb collisions. We stress once more the value of pPh and 
Php control runs, and possibly also lower energy pp runs. We showed examples indicating the need of 
both pPb and Pbp runs. The pPb runs are useful to determine the small-x gluon density of Pb while 
Php runs provide a better normalization benchmark for inclusive charm pT spectra. A concrete plan of 
measurements to benchmark the predictions against real data and to determine the residual systematic 
uncertainties in the extrapolations remains to be outlined. 
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3. QUARKONIUM BASELINE PREDICTIONS ^ 

3.1 Introduction 

In this section, we discuss quarkonium production in pp, pA, and AB colUsions. Here we only consider 
the effects of nuclear shadowing on quarkonium production. Other possible consequences of nuclear 
collisions such as absorption and energy loss in nuclear matter and finite temperature effects are discussed 
later in this chapter. 

Early studies of high-pj- quarkonium production revealed that direct J/ tp production in the color 
singlet model was inadequate to describe charmonium hadroproduction. However, it may be able to 
describe charmonium production in cleaner environments such as photoproduction. Given this difficulty 
with hadroproduction, several other approaches have been developed. The Color Evaporation Model 
(CEM), discussed in section 13.21 treats heavy flavor and quarkonium production on an equal footing, 
the color being 'evaporated' through an unspecified process which does not change the momentum. 
Nonrelativistic QCD (NRQCD), sections 13.31 and 13.341 is an effective field theory in which the short 
distance partonic interactions produce QQ pairs in color singlet or color octet states and nonperturbative 
matrix elements describe the evolution of the QQ pair into a quarkonium state. The color singlet model 
essentially corresponds to dropping all but the first term in the NRQCD expansion (see sect. 13. 3H . 
Finally, the Comover Enhancement Scenario (CES), section l331 bases its predictions on the color singlet 
model enhanced by absorption of gluons from the medium. 

3.2 Quarkonium production in the Color Evaporation Model ^ 

To better understand quarkonium suppression, it is necessary to have a good estimate of the expected 
yields. However, there are still a number of unknowns about quarkonium production in the primary 
nucleon-nucleon interactions. In this section, we discuss quarkonium production in the color evaporation 
model (CEM) and give predictions for production in pp and AA interactions at the LHC. 

The CEM was first discussed a long time ago n^ TOl and has enjoyed considerable phenomeno- 
logical success. In the CEM, the quarkonium production cross section is some fraction Fc of all QQ 
pairs below the HH threshold where H is the lowest mass heavy flavor hadron. Thus the CEM cross 
section is simply the QQ production cross section with a cut on the pair mass but without any constraints 
on the color or spin of the final state. The produced QQ pair then neutralizes its color by interaction 
with the collision-induced color field — "color evaporation". The Q and the Q either combine with light 
quarks to produce heavy-flavored hadrons or bind with each other to form quarkonium. The additional 
energy needed to produce heavy-flavored hadrons when the partonic center of mass energy, V§, is less 
than 2mH, the heavy hadron threshold, is obtained nonperturbatively from the color field in the interac- 
tion region. Thus the yield of all quarkonium states may be only a small fraction of the total QQ cross 
section below 2mH- At leading order, the production cross section of quarkonium state C in an AB 
collision is 

^CEM = ds dxidX2 fi/Aixi,fJ^'^) /j/B(x2,At^) (Tij(s) 6{s - XIX2S) , (3) 

where A and B can be any hadron or nucleus, ij = qq or gg, &ij{s) is the ij —>■ QQ subprocess cross 
section and fijA{x, ji^) is the parton density in the hadron or nucleus. The total QQ cross section takes 
s ^ s. In a collision where either or both A and 5 is a nucleus, we use the EKS98 parameterization of 
the nuclear modifications of the parton densities fT6lfT7ll to model shadowing effects. Since quarkonium 
production is gluon dominated, isospin is negligible. However, because shadowing effects on the gluon 
distribution may be large, they could strongly influence the results. 

''Section coordinator: R. Vogt. 
''Author: R. Vogt. 
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The fraction Fq must be universal so that, once it is fixed by data, the quarkonium production 
ratios should be constant as a function of y/s, y and px- The actual value of Fq depends on the heavy 
quark mass, mq, the scale, /i^, the parton densities and the order of the calculation. It was shown in 
Ref. I20II that the quarkonium production ratios were indeed constant, as expected by the model. 

Of course the leading order calculation in Eq. ^ is insufficient to describe high pT quarkonium 
production since the QQ pair px is zero at LO. Therefore, the CEM was taken to NLO 120 21.1 using 
the exclusive QQ hadroproduction code of Ref. |5|. At NLO in the CEM, the process gg — > gQQ is 
included, providing a good description of the quarkonium pr distributions at the Tevatron 1211 . In the 
exclusive NLO calculation of Ref. |5 1, both the Q and Q variables are available at each point in the phase 
space over which the integration is carried out. Thus, one is not forced to choose the renormalization and 
factorization scales proportional to the heavy-quark mass. A popular choice is for example /i oc niT = 
\JfnQ + Pj', with pj, = Q.bipj,^ + Pj-_). For simplicity, we refer to oc mq in the following but the 
proportionality to uit is implied. 

We use the same parton densities and parameters that agree with the QQ total cross section data, 
given in Table |5l to determine Fc for J/ip and T production. The fit parameters |22 23 1 for the parton 
densities [24 25 26 1, quark masses and scales are given in Table|5]while the QQ cross sections calculated 
with these parameters are compared to pp QQ and ir^p QQ data in Fig. [8] 



Table 5: Parameters used to obtain the 'best' agreement to the QQ cross sections. 



CC 


bb 


Label 


PDF 


mc (GeV) 


fi/rric 


Label 


PDF 


TTife (GeV) 


fi/rrib 


ijjl 


MRST HO 


1.2 


2 


Tl 


MRST HO 


4.75 


1 


%l)2 


MRST HO 


1.4 


1 


T2 


MRST HO 


4.5 


2 


V'3 


CTEQ 5M 


1.2 


2 


T3 


MRST HO 


5.0 


0.5 




GRV 98 HO 


1.3 


1 


T4 


GRV 98 HO 


4.75 


1 



We now describe the extraction of Fc for the individual quarkonium states. The J/i/j has been 
measured in pp and pA interactions up to ^/s = 63 GeV. The data are of two types: the forward cross 
section, a{xF > 0), and the cross section at zero rapidity, da/dy\y=o. All the cross sections are inclusive 
with feeddown from Xc and tp' decays. To obtain Fj^^ for inclusive J ftp production, the normalization 
of Eq. ^ is fit for the cc parameters in Table |5l The comparison of cr to the xp > data for all 
four fits is shown on the left-hand side of Fig. 15] The ratios of the direct production cross sections to the 
inclusive J/ip cross section can be determined from data on inclusive cross section ratios and branching 
fractions. These direct ratios, Rc, given in Tabled are multiplied by the inclusive fitted Fj/^ to obtain 
the direct production fractions, F^" = Fj/^Rc- 

Table 6: Direct quarkonium production ratios, Rc = o"c"^/o"c''^ where C' = J/ip and T. From Ref. 1271 . 



J/V 




Xcl 


Xc2 


T 








Xb{2P) 


Rc 0.62 


0.14 


0.60 


0.99 


0.52 


0.33 


0.20 


1.08 


0.84 



The same procedure, albeit somewhat more complicated due to the larger number of bottomonium 
states below the BB threshold, is followed for the bottomonium. For most data below ^/s = 100 GeV, 
the three bottomonium S states were either not separated or their sum was reported. No x i^-integrated 
cross sections were available so that we fit the CEM T cross section to the effective lepton pair cross 
section at y = for the three T{nS) states. The exti^acted fit fraction is labeled Fj^ t- The comparison 

of cr^^^ with Fj2 t ^ the data for all parameter sets in Table|5lis shown on the right-hand side of Fig.|9l 
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Fig. 8: The cc, (a) and (b), and hb, (c) and (d), total cross section data in pp and Tr^p interactions compared to NLO calculations. 
In (a) and (b), we show the results obtained with parameter choices t/;! (solid), ?/>2 (dashed), i/)3 (dot-dashed) and ?/)4 (dotted). 
In (c) and (d), we show those obtained with Tl (solid), T2 (dashed), T3 (dot-dashed) and T4 (dotted). 



Using the individual branching ratios of the T, T' and T" to lepton pairs and the total cross sections 
reported by CDF 128J, it is possible to extract the inclusive T fit fraction, Fx. The direct production 
ratios obtained in Ref. l29l have been updated in Ref. 1271 using recent CDF Xh data. The resulting 
direct to inclusive T ratios, are also given in Table |6l The subthreshold hb cross section is then 
multiplied by Fq"^ = FyRc to obtain the direct bottomonium cross sections. 

The energy dependence shown in Fig. |9l for both states is well reproduced by the NLO CEM. 
All the fits are equivalent for ^/s < 100 GeV but differ by up to a factor of two at 5.5 TeV. Since the 
pr-integrated T cross sections have been measured at the Tevatron, the ^/s range of the extrapolation 
to the LHC is rather small for the T. The high energy T data seem to agree best with the energy 
dependence obtained with the parameter sets Tl and T2. A similar check cannot be made for the J/ip 
because the high lepton px cut excludes J/ijj acceptance for = at the Tevatron. However, the good 
agreement with the lower energy data results in less than a factor of two difference between the four 
cases at ^/s = 5.5 TeV. 



Table 7: The inclusive J/ip production fractions obtained from data for the cases given in Table |5| The direct charmonium 
cross sections for pp collisions at 5.5 TeV are also given. 







Case Fjf^, 


J/i' Xcl Xc2 1p' 


0.0144 
iIj2 0.0248 
V'S 0.0155 
i^A 0.0229 


19.0 18.3 30.2 4.3 
12.4 12.0 19.8 2.8 
22.2 21.6 35.6 5.0 
19.8 19.3 31.8 4.5 
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Vs" (GeV) Vs" (GeV) 

Fig. 9: Forward J/?/; (left) and combined T + T' + T" inclusive (right) cross sections calculated to NLO in the CEM. On the 
left-hand side, we show the results obtained with parameter choices tpl (solid), tp2 (dashed), ip3 (dot-dashed) and tpA (dotted). 
On the right-hand side, we show those obtained using Tl (solid), T2 (dashed), T3 (dot-dashed) and T4 (dotted). 

The pp cross sections obtained for the individual states at 5.5 TeV are shown in Tables and |8] 
along with the values of the inclusive Fc- We give both Fj-^ y and Fy for bottomonium. Only the direct 

cross sections are given in the tables. To obtain cr'^"', we use the production ratios in Table |6l 

Table 8: The inclusive T production fractions obtained from data for the cases given in Table|5| The direct bottomonium cross 
sections for pp collisions at 5.5 TeV are also given. 











(nb) 




Case 




Ft 


T 






Xb{lP) 


Xbi2P) 


Tl 


0.000963 


0.0276 


188 


119 


72 


390 


304 


T2 


0.000701 


0.0201 


256 


163 


99 


532 


414 


T3 


0.001766 


0.0508 


128 


82 


49 


267 


208 


T4 


0.000787 


0.0225 


145 


92 


56 


302 


235 



The range of the fit parameters allows us to explore the dependence of Fj/^ and Fx on uiq and 
fi. The range of fit parameters is limited because we only choose parameters that are in relatively good 
agreement with the QQ total cross sections. The Fj/^ obtained with the parameter sets tpl and ip3, 
employing the same mass and scale, are rather similar but a'^^^ differs by 20% at 5.5 TeV. For the same 
PDFs but with a lower scale, case 1^2 relative to tpl, Fj/^, is about a factor of two larger. However, the 
cross sections at 5.5 TeV differ only by 50%. Changing the mass and PDF at the same scale, cases ip3 
and V'4, does not change Fj substantially. 

Since the b quark is more massive, the scale dependence can be more sensibly explored. The bb 
cross sections for cases Tl, T2 and T3 are essentially equivalent. However, the resulting Fx differs by 
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a factor of 2.5 with the highest mj, giving the largest Fr and the lowest a . There is a factor of two 
between the corresponding T cross sections. For different PDFs but the same mass and scale, Tl and 
T4, the fitted Fr's differ by only ~ 10%. 

Table 9: The direct cross section per nucleon pair and the dilepton cross section per nucleon multipMed by for the minimum 
bias lepton pair cross section. The resuhs are given for tpl. We compare pp to pA and AA interactions. 





cj'^"'/nucleon pair (/ib) 




(ah) 


System 


V^(TeV) 




Xcl 






J/i> 




PP 


14 


32.9 


31.8 


52.5 


7.43 


3.18 


0.057 


PP 


8.8 


25.0 


24.2 


39.9 


5.65 


2.42 


0.044 


pPh 


8.8 


19.5 


18.9 


31.1 


4.40 


392.3 


7.05 


PP 


7 


21.8 


21.1 


34.9 


4.93 


2.11 


0.038 


0+0 


7 


17.6 


17.0 


28.1 


3.98 


436.2 


7.84 


PP 


6.3 


20.5 


19.9 


32.8 


4.63 


1.99 


0.036 


Ar+Ar 


6.3 


15.0 


14.5 


23.9 


3.38 


2321 


41.7 


PP 


6.14 


20.2 


19.6 


32.3 


4.56 


1.96 


0.035 


Kr+Kr 


6.14 


13.7 


13.2 


21.8 


3.08 


9327 


167.6 


PP 


5.84 


19.6 


19.0 


31.3 


4.42 


1.90 


0.034 


Sn+Sn 


5.84 


12.8 


12.4 


20.4 


2.89 


17545 


315.2 


PP 


5.5 


18.9 


18.3 


30.2 


4.26 


1.83 


0.033 


Pb+Pb 


5.5 


11.7 


11.3 


18.7 


2.64 


48930 


879 



We now show the effects of nuclear shadowing on the total cross sections for one particular set 
of parameters. We choose parameter set ipl for charmonium and set Tl for bottomonium. The results 
are given in Tables l9l and ITOl In the middle of each table, the direct cross sections per nucleon pair are 
given for all states. For every applicable energy, the pp and AA, or for 8.8 TeV, the pp and pA cross 
sections are compared to directly show the shadowing effects. We can see that the effects are largest for 
charmonium and for the heaviest nuclei even though these are at the lowest energies and thus the highest 
X. Shadowing effects on the gluon distributions do change significantly at this low x L16. ,IH . In all 
cases, the effect is less than a factor of two. 

On the right-hand side of the tables, the inclusive cross sections are multiplied by the lepton pair 
branching ratios. They are also multiplied by A"^ to reproduce the minimum bias cross sections. The 
reduction due to shadowing is then the only nuclear dependence included. We have not added in nuclear 
absorption effects in cold matter, discussed in section 1411 

The J/ip cross sections reported here are a factor of two or more lower than those calculated in 
Ref. |20|. This should not be a surprise because the PDFs used in those calculations were available 
before the first low-x HERA data and generally overestimated the increase at low x. At x ~ 10~^, the 
MRS D-' gluon density [30| is nearly a factor of five greater than the MRST gluon density |24l| based on 
more recent HERA data that better constrain the low-x gluon density. The differences between the GRV 
HO II3TI and GRV 98 HO ll26l are smaller. The gluon densities used in this study are compared to those 
used in Ref. ||^ in Fig.lTOl 

The direct J/V' and T rapidity distributions in Pb+Pb interactions at 5.5 TeV/nucleon for all the 
parameter choices are compared in the left-hand side of Fig. The rapidity distributions reflect the 
differences in the total cross sections quite well. The 'corners' in the J/?/; rapidity distributions at \y\ ~ 4 
occur at X ~ 10^^, the lowest x for which the MRST and CTEQ5 densities are valid. The behavior of 
the parton gluon densities for x < 10~^ varies significantly, as shown in Fig.QOl The minimum x of the 
GRV 98 densities is 10^^ so that no problems are encountered for this set. The MRST densities below 
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Table 10: The direct cross section per nucleon pair and the dilepton cross section per nucleon multiplied by for the minimum 
bias lepton pair cross section. The results are given for Tl. We compare pp to pA and AA interactions. 









/nucleon pair (ph) 




Bcr'°M^ (^b) 


System 




T 










T 






Tin 
Ft' 


14 


0.43 


0.27 


0.16 


0.89 


0.69 


0.020 


0050 


0030 


Tifi 


8.8 


0.29 


0.18 


0.11 


0.60 


0.47 


0.014 


0.0040 


0020 


pPb 


8.8 


0.25 


0.16 


0.097 


0.52 


0.41 


2.51 


0.65 


0.37 


pp 


7 


0.23 


0.15 


0.090 


0.48 


0.38 


0.011 


0.0029 


0.0016 


0+0 


7 


0.21 


0.13 


0.081 


0.44 


0.34 


2.57 


0.66 


0.38 


pp 


6.3 


0.21 


0.14 


0.082 


0.44 


0.34 


0.010 


0.0026 


0.0015 


Ar+Ar 


6.3 


0.18 


0.12 


0.070 


0.38 


0.29 


13.8 


3.59 


2.02 


pp 


6.14 


0.21 


0.13 


0.080 


0.43 


0.33 


0.0099 


0.0026 


0.0014 


Kr+Kr 


6.14 


0.17 


0.11 


0.066 


0.35 


0.28 


57.4 


14.8 


8.38 


pp 


5.84 


0.20 


0.12 


0.076 


0.41 


0.32 


0.0094 


0.0024 


0.0014 


Sn+Sn 


5.84 


0.16 


0.10 


0.062 


0.33 


0.26 


108.1 


28.0 


15.8 


pp 


5.5 


0.19 


0.12 


0.070 


0.39 


0.30 


0.0090 


0.0020 


0.0013 


Pb+Pb 


5.5 


0.15 


0.094 


0.057 


0.31 


0.24 


304 


78.8 


44.4 



X = 10^^ are fixed to the density at this minimum value and are thus constant for lower values of x. On 
the other hand, the CTEQ5M distributions turn over and decrease for x < 10^^, causing the steep drop 
in the rapidity distributions at high \y\. Only the GRV98 HO distributions are smooth over all y. 

On the right-hand side of Fig. the AA/pp ratios are compared for all the AA combinations 
given in Tables l9l and fTOl Since the rapidity distributions are not smooth due to the Monte Carlo integra- 
tion of Ref. (5 1, the ratios enhance the fluctuations. The 'comers' do not appear in the ratios because the 
change in slope occurs at the same point for pp and AA at the same energy. The biggest effect of shad- 
owing is at midrapidity when both x values are small. As the rapidity increases, the ratios also increase 
since, e.g. at large y, xi is large and in the antishadowing region, reducing the shadowing effect in the 
product. The effect on the J ftp is the largest, from a ~ 45% effect on Pb-i-Pb at y = to a 20% effect 
at y = for O-i-O. The overall effect on the T is lower since the x values probed, as well as the scale, 
are larger. The evolution of the shadowing parameterization decreases the effect. Thus the T result is 
~ 25% for Fh+Fh and only ~ 10% for 0+0. 

The rapidity-integrated pT distributions of direct quarkonium production are compared for all fit 
parameters in Pb-i-Pb collisions at 5.5 TeV in Fig.[^ The cross sections are not shown all the way down 
to pt = because we have not included any intrinsic fc-r broadening. Broadening effects on quarkonium 
are discussed in the pA chapter of this report. The distributions are all fairly similar but changing the 
mass and scale has an effect on the slope, as is particularly obvious for the T. The highest nib, T'^S, has 
the hardest slope while that of the lowest mt, T2, decreases the fastest with px. 

The AA/pp ratios are also shown in Fig.^J The fluctuations are again large but the general trend 
is clear. The ratios at low pT are similar to those at midrapidity and increase to unity around px ~ 15 — 20 
GeV for both the J/V' and T. 
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Fig. 10: Gluon distribution functions in the proton at the scale of the charmonium calculations. The lower solid curve is the 
scale independent (1 — x)^, the other solid curve employs the MRST HO distributions with fi — 2.4 GeV, the dashed, GRV 
98 HO with fi = 1.3 GeV, the dot-dashed, MRSD-' with n = 2.4 GeV, the dotted, GRV HO with = 1.3 GeV and the 
dot-dot-dot-dashed, CTEQ 5M with /i = 2.4 GeV. 

3.3 Quarkonium production in Non-Relativistic QCD ^ 

3.31 The NRQCD Factorization Method 

In both heavy-quarkonium decays and hard-scattering production, large energy-momentum scales ap- 
pear. The heavy-quark mass mq is much larger than Aqcd> and, in the case of production, the trans- 
verse momentum pT can be much larger than Aqcd as well. Thus, the associated values of as are much 
less than one: as(mc) ~ 0.25 and as(m{,) ~ 0.18. Therefore, one might hope that it would be possi- 
ble to calculate the rates for heavy quarkonium production and decay accurately in perturbation theory. 
However, there are clearly low-momentum, nonperturbative effects associated with the dynamics of the 
quarkonium bound state that invalidate the direct application of perturbation theory. 

In order to make use of perturbative methods, one must first separate the short-distance/high- 
momentum, perturbative effects from the long-distance/low-momentum, nonperturbative effects — a pro- 
cess which is known as "factorization." One convenient way to carry out this separation is through the use 
of the effective field theory Nonrelativistic QCD (NRQCD) |32l|3ll|34|. NRQCD reproduces full QCD 
accurately at momentum scales of order rnqv and smaller, where v is heavy-quark velocity in the bound 
state in the center-of-mass (CM) frame, with 0.3 for charmonium and 0.1 for bottomonium. 
Virtual processes involving momentum scales of order mq and larger can affect the lower-momentum 
processes. Their effects are taken into account through the short-distance coefficients of the operators 
that appear in the NRQCD action. 

Because QQ production occurs at momentum scales of order mq or larger, it manifests itself 
in NRQCD through contact interactions. As a result, the quarkonium production cross section can be 

^Authors: G.T. Bodwin, Jungil Lee and R. Vogt. 
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Fig. 11: The direct J/tp and T rapidity distributions in Pb+Pb collisions (left) and the AA/pp ratios (right). All the plots are 
normalized per nucleon. On the left-hand side, the J/tp calculations are tpl (solid), ?/>2 (dashed), ip3 (dot-dashed) and ip4 
(dotted) while the T calculations are Tl (solid), T2 (dashed), T3 (dot-dashed) and T4 (dotted). The ratios on the right-hand 
side use ipl for the J/ip and Tl for the T. The ratios are given for the maximum AA energy: ^4 = Pb (solid), Sn (dashed), Ar 
(dot-dashed) and O (dot-dash-dash-dashed). 



written as a sum of the products of NRQCD matrix elements and short-distance coefficients: 

^(^) = E^(0|O^(A)|0). (4) 

Here, H is the quarkonium state, A is the ultraviolet cutoff of the effective theory, the F„ are short- 
distance coefficients, and the are four-fermion operators, whose mass dimensions are dn- A formula 
similar to Eq. Q exists for the inclusive quarkonium annihilation rate l34l . 

The short-distance coefficients F„ (A) are essentially the process-dependent partonic cross sections 
to make a QQ pair. The QQ pair can be produced in a color-singlet state or in a color-octet state. 
The short-distance coefficients are determined by matching the square of the production amplitude in 
NRQCD to full QCD. Because the QQ production scale is of order uiq or greater, this matching can be 
carried out in perturbation theory. 

The four-fermion operators in Eq. Q create a QQ pair, project it onto an intermediate state that 
consists of a heavy quarkonium plus anything, and then annihilate the QQ pair. The vacuum matrix 
element of such an operator is the probability for a QQ pair to form a quarkonium plus anything. These 
matrix elements are somewhat analogous to parton fragmentation functions. They contain all of the 
nonperturbative physics that is associated with evolution of the QQ pair into a quarkonium state. 



22 




Fig. 12: The direct J/ip and T px distributions in Pb+Pb collisions (left) and the AAlpp ratios (right). All the plots are 
normalized per nucleon. On the left-hand side, the J/ip calculations are ipl (solid), ?/>2 (dashed), ip?) (dot-dashed) and ?/>4 
(dotted) while the T calculations are Tl (solid), T2 (dashed), T3 (dot-dashed) and T4 (dotted). The ratios on the right-hand 
side use xpl for the J/i/; and Tl for the T. The ratios are given for the maximum AA energy: A = Pb (solid), Sn (dashed), Ar 
(dot-dashed) and O (dot-dash-dash-dashed). 



Both color-singlet and color-octet four-fermion operators appear in Eq. @. They correspond, 
respectively, to the evolution of a QQ pair in a relative color-singlet state or a relative color-octet state 
into a color-singlet quarkonium. If we drop all of the color-octet contributions in Eq. (|4]i, then we have 
the color-singlet model 1 35 1. In contrast, NRQCD is not a model, but a rigorous consequence of QCD in 
the Umit Aqqd /mg 0. 

The NRQCD decay matrix elements can be calculated in lattice simulations f36l El or deter- 
mined from phenomenology. However, at present, the production matrix elements must be obtained 
phenomenologically, as it is not yet known how to formulate the calculation of production matrix el- 
ements in lattice simulations. In general, the production matrix elements are different from the decay 
matrix elements. However, in the color-singlet case, the production and decay matrix elements can be 
related through the vacuum-saturation approximation, up to corrections of relative order li34J . 

An important property of the matrix elements, which greatly increases the predictive power of 
NRQCD, is the fact that they are universal, i.e., process independent. NRQCD f -power-counting rules 
organize the sum over operators in Eq. Q as an expansion in powers of v. Through a given order in v, 
only a limited number of operator matrix elements contribute. Furthermore, at leading order in v, there 
are simplifying relations between operator matrix elements, such as the heavy-quark spin symmetry [34| 
and the vacuum-saturation approximation L34.1 . that reduce the number of independent phenomenological 
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parameters. In contrast, the CEM ignores the hierarchy of matrix elements in the v expansion. 

The proof of the factorization formula (HJ relies both on NRQCD and on the all-orders perturbative 
machinery for proving hard-scattering factorization. A detailed proof does not yet exist, but work is in 
progress |38|. Corrections to the hard-scattering part of the factorization are thought to be of order 
{mQv/pT)'^, not {niQ/pT)'^, in the unpolarized case and of order rnqv/pT, not itiq/pt, in the polarized 
case. It is not known if there is a factorization formula at low px or for the -integrated cross section. 
The presence of soft gluons in the quarkonium binding process makes the application of the standard 
factorization techniques problematic at low pT- 

In the decay case, the color-octet matrix elements can be interpreted as the probability to find the 
quarkonium in a Fock state consisting of a QQ pair plus some gluons. It is a common misconception that 
color-octet production proceeds, like color-octet decay, through a higher Fock state. However, in color- 
octet production, the gluons that neutralize the color are in the final state, not the initial state. There is 
a higher-Fock-state process, but it requires the production of gluons that are nearly coUinear to the QQ 
pair, and it is, therefore, suppressed by additional powers of v. 

In practical theoretical calculations of the quarkonium production and decay rates, a number of 
significant uncertainties arise. In many instances, the series in ag and in v of Eq. Q converge slowly, 
and the uncertainties from their truncation are large — sometimes of order 100%. In addition, the matrix 
elements are often poorly determined, either from phenomenology or lattice measurements, and the im- 
portant linear combinations of matrix elements vary from process to process, making tests of universality 
difficult. There are also large uncertainties in the heavy-quark masses (approximately 10% for nic and 
5% for mft, for the mass ranges used in the calculations) that can be very significant for quarkonium rates 
proportional to a large power of the mass. 

3.32 Experimental Tests of NRQCD Factorization 

Here, we give a brief review of some of the successes of NRQCD, as well as some of the open questions. 
We concentrate on hadroproduction results for both unpolarized and polarized production. We also 
discuss briefly some recent two-photon, e+e~, and photoproduction results. 

Using the NRQCD-factorization approach, one can obtain a good fit to the high-p^ CDF data [391, 
while the color-singlet model under predicts the data by more than an order of magnitude. (See Fig.fTSl) 
The pt dependence of the unpolarized Tevatron charmonium data has been studied under a number of 
model assumptions, including LO collinear factorization, parton-shower radiation, smearing, and kx 
factorization. (See Ref. i40ll for a review.) 

Several uncertainties in the theoretical predictions affect the extraction of the NRQCD charmoni- 
um-production matrix elements from the data. There are large uncertainties in the theoretical predictions 
that arise from the choices of the factorization scale, the renormalization scale, and the parton distribu- 
tions. The extracted values of the octet matrix elements are very sensitive to the small-p^ behavior of 
the cross section and this, in turn, leads to a sensitivity to the behavior of the small-x gluon distribution. 
Furthermore, the effects of multiple soft-gluon emission are important, and their omission in the fixed- 
order perturbative calculations leads to overestimates of the matrix elements. Effects of higher-order 
corrections in as are a further uncertainty in the theoretical predictions. Similar theoretical uncertainties 
arise in the extraction of the NRQCD production matrix elements for the T |41] states, but, owing to 
large statistical uncertainties, they are less significant for the fits than in the charmonium case. 

At large pT (jpr ^ 4mc for the J / Tp) the dominant quarkonium-production mechanism is gluon 
fragmentation into a QQ pair in a '^5i color-octet state. The fragmenting gluon is nearly on mass shell 
and is, therefore, transversely polarized. Furthermore, the velocity-scaling rules predict that the color- 
octet QQ state retains its transverse polarization as it evolves into 5-wave quarkonium (^T\, up to cor- 
rections of relative order v'^. Radiative corrections, color-singlet production, and feeddown from higher 
states can dilute the quarkonium polarization [43j|33|351|46,,47J. Despite this dilution, a substantial po- 
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Fig. 13: J/ijj cross section as a function of px- The data points are from tlie CDF measurement [39,J. The solid curve is the 
NRQCD-factorization fit to the data given in Ref. |40| . The other curves give various contributions to the NRQCD-factorization 
fit. From Ref. |40|. 

larization is expected at large pT- Its detection would be a "smoking gun" for the presence of color-octet 
production. In contrast, the color-evaporation model predicts no quarkonium polarization. The CDF 
measurement of the J/ if) and ip' polarization as a function of pT |48| is shown in Fig. [M] along with 
the NRQCD factorization prediction ^44l l45ll46l . The analysis of ijj' polarization is simpler than for the 
J/ip, since feeddown does not play a role. However, the statistics are not as good for the if)' . The degree 
of polarization is a = (1 — 3^) /(I + C), where ^ is the fraction of events with longitudinal polarization. 
a = I corresponds to 100% transverse polarization, and a = — 1 corresponds to 100% longitudinal 
polarization. The observed polarization is in relatively good agreement with the prediction, except in the 
highest Pt bin, although the prediction of increasing polarization with increasing pT is not in evidence. 

Because the polarization depends on a ratio of matrix elements, some of the theoretical uncertain- 
ties are reduced compared with those in the production cross section, and, so, the polarization is probably 
not strongly affected by multiple soft-gluon emission or K factors. Contributions of higher order in ag 
could conceivably change the rates for the various spin states by a factor of two. Therefore, it is impor- 
tant to carry out the NLO calculation, which involves significant computational difficulties. It is known 
that order-w^ corrections to parton fragmentation into quarkonium can be quite large l49ll . If spin-flip 
corrections to the NRQCD matrix elements, which are nominally suppressed by powers of v, are also 
large, perhaps because the velocity-scaling rules need to be modified, then spin-flip contributions could 
significantly dilute the J/ip polarization. Nevertheless, in the context of NRQCD, it is difficult to see 
how there could not be substantial charmonium polarization for pT > ^nic- 

Compared to the J/ V'-polarization prediction, the T-polarization prediction has smaller v-expan- 
sion uncertainties. However, because of the higher T mass, it is necessary to go to higher pT to insure that 
fragmentation dominates and that there is substantial polarization. Unfortunately, the current Tevatron 
data run out of statistics in this high-p^ region. CDF finds that a = —0.12 it 0.22 for 8 < < 20 GeV 
ll50l . which is consistent with both the NRQCD-factorization prediction lISTl and the zero-polarization 
prediction of the CEM. There are also discrepancies between the polarizations observed in fixed-target 
experiments and the NRQCD predictions. 
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Fig. 14: Left-hand side: J/t/; polarization at the Tevatron. The band is the total NRQCD-factorization prediction. The other 
curves give the contributions from feeddown from higher charmonium states. Right-hand side: t/;' polarization at the Tevatron. 
The bands give various NRQCD-factorization predictions. The data points are from the CDF measurement 1481 . From Ref . 1461 . 

Calculations of inclusive J/ and T production in 77 collisions l52l l53l have been compared 
with LEP data I54ll55l l56l. Both the J /ij; and T measurements favor the NRQCD predictions over those 
of the color-singlet model. 

Belle L57| and BaBar have also measured the J /if) total cross sections in e+e~ J/ijjX. 
The results of the two experiments are incompatible with each other, but they both seem to favor NRQCD 
over the color-singlet model. A surprising new result from Belle IWI is that most of the produced J/?/;'s 
are accompanied by an additional ccpair: cr(e"'"e~ — > J cc) j o[€^ e~ — > J/'4)X) = 0.591q ± 0.12. 
Perturbative QCD plus the color-singlet model predict that this ratio should be about 0.1 1.60 J . There 
seems to be a major discrepancy between theory and experiment. However, the order-a^ calculation 
lacks color-octet contributions, including those that produce Jj^cc. Although these contributions are 
suppressed by 0.1, it is possible that the short-distance coefficients are large. In other results, 

the angular distributions favor NRQCD, but the polarization measurements show no evidence of the 
transverse polarization that would be expected in color-octet production. However, the center-of-mass 
momentum is rather small, and, hence, one would not expect the polarization to be large. 

Quarkonium production has also been measured in inelastic photoproduction 1611 l62l and deep- 
inelastic scattering (DIS) 163^, at HERA. The NRQCD calculation deviates from the data near large 
photon-momentum fractions, owing to the large LO color-octet contribution. The NLO color-singlet 
result agrees with the data over all momentum fractions, as well as with the data as a function of py. 
See Ref. l40l for a more complete review. In the case of deep-inelastic scattering, the and pT depen- 
dences are in agreement with NRQCD, but the results are more ambiguous for the dependence on the 
longitudinal momentum fraction. 

J.Ji Quarkonium Production in Nuclear Matter 

The existing factorization "theorems" for quarkonium production in hadronic collisions are for cold 
hadronic matter. These theorems predict that nuclear matter is "transparent" for J/V' production at large 
•pT- That is, at large p^, all of the nucleai" effects are contained in the nuclear parton distributions. The 
corrections to this transparency are of order {rnQv/pT)"^ for unpolarized cross sections and of order 
mQv/pT for polarized cross sections. 

The effects of transverse-momentum kicks from multiple elastic collisions between active partons 
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and spectators in tiie nucleons are among tiiose effects tiiat are suppressed by {rngv / pr)'^ . Nevertheless, 
these multiple-scattering effects can be important because the production cross section falls steeply with 
Pt and because the number of scatterings grows linearly with the path length through nuclear matter. 
Such elastic interactions can be expressed in terms of eikonal interactions l65l or higher-twist matrix 
elements [66|. 

Inelastic scattering of quarkonium by nuclear matter is also an effect of higher order in {rnqv/pT)'^. 
However, it can become dominant when the amount of nuclear matter that is traversed by the quarkonium 
is sufficiently large. Factorization breaks down when 

min(zQ,ZQ)P^ 
~ M^A;2(tot) ' ^' 

where L is the length of the quarkonium path in the nucleus. Ma is the mass of the nucleus, z is the 
parton longitudinal momentum fraction, Ph is the momentum of the quarkonium in the parton CM 
frame, and A;r(tot) is the accumulated transverse momentum "kick" from passage through the nuclear 
matter. This condition for the breakdown of factorization is similar to "target-length condition" in Drell- 
Yan production [67, 681. Such a breakdown is observed in the Cronin effect at low pT and in Drell-Yan 
production at low Q^, where the cross section is proportional to A", and a < 1. 

It is possible that multiple-scattering effects may be larger for color-octet production than for 
color-singlet production. In the case of color-octet production, the pre-quarkonium QQ system carries 
a nonzero color charge and, therefore, has a larger amplitude to exchange soft gluons with spectator 
partons. 

At present, there is no complete, rigorous theory to account for all of the effects of multiple 
scattering and we must resort to "QCD-inspired" models. A reasonable requirement for models is that 
they be constructed so that they are compatible with the factorization result in the large-p-p limit. Many 
models treat interactions of the pre-quarkonium with the nucleus as on-shell (Glauber) scattering. This 
assumption should be examined carefully, as on-shell scattering is known, from the factorization proofs, 
not to be a valid approximation in leading order in {niQv/pT)'^. 



3. 34 NRQCD Predictions for the LHC 

In this section, we shall use the formalism of NRQCD to give predictions for quarkonium production 
in the LHC energy range. We rewrite the cross section in Eq. © for the inclusive production of a 
charmonium state H as follows: 

a{H) = (O^) , (6) 

n 

where a'^^^)" = Fn{A)/m'Q'^, {O^) = (OlO^lO), and n runs over all the color and angular mo- 
mentum states of the QQ pair. The cross sections a^QQ)" can be calculated in perturbative QCD. All 
dependence on the final state H is contained in the nonperturbative NRQCD matrix elements (O^). 

The most important matrix elements for J/iJj = ip{lS) and tj;' = ip{2S) production can be re- 
duced to the color-singlet parameter {Of^"'^\^ Si)) and the three color-octet parameters {Oq^"'^\^Si)), 
{Og^"'^\^ Sq)) , and {Og^'^^\^ Pq)) . Two of the three color-octet matrix elements only appear in the 
linear combination 

Mfi-S) ^ ^k/ml){ot^''''\'Po)) + (of """V^o)) • (7) 

The value of k is sensitive to the pr dependence of the fit. At the Tevatron, k ^ 3. Fits to fixed-target total 
cross sections give larger values, k 6-7 |69|. The most important matrix elements for Xcj production 
can be reduced to a color-singlet parameter (O^^"" (^Pq)) and a single color-octet parameter (Og"" (^5*1)). 
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Table 11: Matrix elements for charmonium production. Note that here (Of) = (O^^"^' (^Si)) for J/ip and but (Of) = 
(Oj (^-Po)) for XcO- Uncertainties are statistical only. From Ref. 1461 . 



H 


(Of) 






k 




J/^ 


1.3 ±0.1 GeV^ 


(4.4 ±0.7) X 10-3 


GeV3 


3.4 


(8.7 ± 0.9) X 10-2 GeV3 




0.65 ± 0.06 GeV^ 


(4.2 ± 1.0) X 10-3 


GeV3 


3.5 


(1.3 ±0.5) X 10-2 GeV3 


XcO 


(8.9 ± 1.3) X 10-2 GeV^ 


(2.3 ±0.3) X 10-3 


GeV3 







These matrix elements are sufficient to calculate the prompt J/ if) cross section to leading order in ag and 
to order relative to the color-singlet contribution. 

In pp coUisions, different partonic processes for QQ production dominate in different pT ranges. 
If pt is of order mg, fusion processes dominate, and, so, the QQ pair is produced in the hard-scattering 
process. These contributions can be written in the form 

a^,{H) = J2jdxi dx2 fi/A{xi,l^^)fj/B{x2,li^)o^lp^'' (Of) , (8) 

where A and B are the incoming hadrons or nuclei. In Eq. we include the parton processes ij — > 
QQ X, where ij = gg, qq, qg and qg, and q = u,d, s. The relevant partonic cross sections o-j^^'^)" are 
given in Refs. ITUHtTII . 

For pt ^ ruQ, the dominant partonic process is gluon fragmentation through the color-octet "^Si 
channel. This contribution can be expressed as 

aFriH) = J2 I dx,dx2dzh/J,[xl,^?)f,,B(x2,^i^)afpf^^''^'^'\z,^4,) (Of) , (9) 

where ^^\z, /Xpj.) is the fragmentation function for a gluon fragmenting into a QQ pair, P/ z is 

the momentum of the fragmenting gluon, P is the momentum of the QQ pair, and /ipr is the fragmen- 
tation scale. The fragmentation process scales as da/dpj. ~ ^/Pt f72ll73l . The fragmentation process 
is actually included in the fusion processes of Eq. (ISj. In the limit px ^ mq, the fusion processes that 
proceed through g* {QQ)?,(^ Si) are well-approximated by the expression At large Pt, one can 
evolve the fragmentation function in the scale ji-pv, thereby resumming large logarithms of pj^/niQ. Such 
a procedure leads to a smaller short-distance factor 1 45 1 and a more accurate prediction at large pT than 
would be obtained by using the fusion cross section However, in our calculations, we employ the 
fusion cross section ©J which leads to systematic over-estimation of the cross section at large px- 

In order to predict the cross section for prompt J/^p production (including Xc and ^p' feeddown) 
at the LHC, we need the values of the NRQCD matrix elements. There have been several previous 
extractions of the color-octet matrix elements |45l|46l|70l|7ill2lll2lll26l from the CDF J/ifj, Xc and ip' 
Pt distributions IBQIITTII . We use the matrix elements given in Ref. f461, which are shown in Table fTTI 
Our calculations are based on the MRST LO parton distributions |78|. In calculating the cross section 
per nucleon for prompt J/V' production in pA or AA collisions, we take /j/^ = fi/pRf. We employ 
the EKS98 parameterization fT6l ITTl for the nuclear shadowing ratio Rf. We evolve Ug at one-loop 
accuracy, and we set /i = {im^ + PrY^"^ and nic = 1.5 GeV. 

There are several sources of uncertainty in our predictions for the cross sections. There are large 
uncertainties in the NRQCD matrix elements themselves. The errors shown in Table are statistical 
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Fig. 15: Differential cross sections per nucleon multiplied by leptonic branching fractions for prompt J /ijj (upper curves), J /ijj 
from Xc (middle curves), and prompt ipi^S) (lower curves) in pp, pPh, and Pb+Pb collisions at yfs = 5.5 TeV. The EKS98 
parameterization I16II17I is employed forpPb and Pb+Pb collisions. 

only. There are additional large uncertainties in the matrix elements that arise from truncations of the 
series in ag and v in the theoretical expressions that are used to extract the matrix elements. The matrix 
elements (08(^5*0)) and {O^i^Po)) are fixed by the data only in the linear combination M^. In the 
present calculation, we take (08(^50)) = xM^ and {Os(^Po)) /ml = (1 — x)M^ /k, use the values 
of k given in Table fTTl and choose x = 1/2. Variation of x between and 1 affects the cross sections 
at low pT by amounts on the order of 5%. There are additional uncertainties in the predicted cross 
sections that arise from the choices of the parton distributions, the charm-quark mass rric, and the scale 
fj,. Because they affect the matrix-element fits, these uncertainties are highly correlated with those of the 
matrix elements. We have not tried to estimate their effects on the predicted cross sections. 

In Fig.^J we show the px distributions per nucleon multiplied by the dilepton branching fractions 
for prompt J/t/j (upper curves), J/ip from Xc decays (middle curves), and prompt t/j' (lower curves) at 
^/s = 5.5 TeV. For pPb and Pb-i-Pb collisions, we use the EKS98 parameterization fT6lfT7l to account 
for the effect of nuclear shadowing. The pp, pPb, and Pb-i-Pb results essentially lie on top of each other 
in Fig. El owing to the many decades covered in the plot. 

In order to display small differences between the distributions, we define the function Rab- 

daAB/dpT - dupp/dpT 
Rab[pt) = 1 , , ■ (10) 

dapp/dpT 

In Fig. El we present Rab as a function of pT- As is shown in Fig. fT6l a). nuclear shadowing increases 
the cross section at large px and decreases it at small px- The deviation of the Pb-i-Pb cross section from 
the pp cross section is twice as large as that seen in the case of pPb collisions. In order to investigate 
the dependence of the shadowing effect on the short-distance cross sections that arise in hadroproduc- 
tion of S'-wave charmonium states in Pb-i-Pb collisions, we plot Rab for all channels separately. [See 
Fig. fT6l b).1 Even though the px dependence the contribution to the cross section of the color-octet ^Si 
channel is quite different from those of the color-octet and ^Pj channels, all three channels show 
the same nuclear effect. The only channel that shows a slightly different behavior is the color-singlet 
channel, which gives a negligible contribution to the cross section. While the differential cross sections 
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Fig. 16: The pr dependence of Rab [Eq. <10H . (a) We compare the results for pPb and Pb+Pb collisions. (The pr dependence 
is stronger in the Pb+Pb result.) (b) We show the dependence of Rab on the various production channels in Pb+Pb collisions 
at ^ = 5.5 TeV. 

in Fig. are strongly dependent on the nonperturbative NRQCD matrix elements, Rab is almost in- 
dependent of the matrix elements, making it a good observable for studying nuclear shadowing at the 
LHC. 

The T rates are somewhat more difficult to calculate because of the many feeddown contributions. 
The matrix elements are also not particularly well known. Since it is unlikely that all the different 
contributions can be disentangled, we follow the approach of Ref. iHTI and compute the inclusive T{nS) 
production cross section 

J 

+ (^Y.(^J + l)da^'^^^('^'^^ {O^^'''\'P0))inc , (11) 

where the last term makes use of heavy-quark spin symmetry to relate all of the octet ^Pj matrix elements 
to the octet ^Pq matrix element. The "inclusive" matrix elements are defined by 

(07^"''^(n))inc = E BH^r(nS){Of'{n)) , (12) 

H 

where i = 1 or 8 for singlet or octet, respectively. The sum over H includes the T{nS) as well as 
all higher states that can decay to T{nS). The branching ratio for H ^ H' decays is Bjj^jji with 
Bh^h = 1- Only Xbi^P) and XbC^P) decays are included; the possibility of feeddown from the as- 
yet unobserved Xfe(3-P) states is neglected. In the linear combination the color-octet matrix 

element from the '^Pq state is neglected, and, so, AfJ^""^^ = (C'J^""^V'S'o))inc- We use mb = 4.77 GeV 
and the MRST LO parton distributions. The values of the inclusive color-singlet matrix elements are 
given in Table (121 and the values of the inclusive color-octet matrix elements, from Ref. [411, are given 
in Table [ni 
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Table 12: Inclusive color-singlet matrix elements for bottomonium production. The errors on the ^Si matrix elements come 
from estimates of the T{nS) decay rate to lepton pairs. The errors on the "^Pj states come from an average over potential- 
model estimates. The inclusive matrix elements are a linear combination of branching ratios, as in Eq. <12> . The S-state matrix 
elements are in units of GeV"^ while the Pj-state matrix elements are in units of GeV^. From Ref. 1411 . 



H 




(Of(3Po))inc 


i(Of(3A)>inc 


i(OfeP2)>inc 


T(35) 


4.3 ±0.9 











T{2S) 


5.0 ±0.7 


0.12 ±0.06 


0.55 ±0.15 


0.42 ±0.10 


T{1S) 


12.8 ± 1.6 


< 0.2 


1.23 ±0.25 


0.84 ±0.15 



Table 13: Inclusive color-octet matrix elements for bottomonium production. The matrix elements were fit using the MRSTLO 
parton distributions. The first set of error bars is from fits to the T pr distributions in the region pr > 8 GeV. The second 
set is associated with the variation of the scales and corresponds to multiplying fi — y'm^+p^ by 2 (upper error) and 0.5 
(lower error). The matrix elements are in units of 10~^ GeV''. From Ref. [4r|. 



H 




T{3S) 
T{2S) 
T(15) 


3.7±1.7lj:^ 7.5 ±4.9+1^ 

19.6 ±6.3^1^ -8.7 ± 11.1:^1:8 

11.7 ±3.0^1^ 18.1 ± 7.2^"!"^ 
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Fig. 17: Differential cross sections per nucleon multiplied by leptonic branching fractions for inclusive T(15') (upper curves), 
T{2S) (middle curves), and prompt T{3S) (lower curves) in pp, pPb, and Pb+Pb collisions at -^s = 5.5 TeV. The EKS98 
parameterization I16II17I is employed forpPb and Pb+Pb collisions. 

In Fig.^J we show the px distributions per nucleon multiphed by the dilepton branching fractions 
for the 3 T S states at ^/s = 5.5 TeV. The feeddown contributions are included as in Eq. (fTTT l. For pPh 
and Pb+Pb collisions, we use the EKS98 parameterization 1 16 in order to account for the effects of 
nuclear shadowing. The pp, pPb, and Pb+Pb results lie essentially on top of each other in Fig. El 

The unusual relative behavior of the T(2S') and T(35) states at both low and high pT is due to the 
fact that the bottomonium matrix elements are not very well determined. For pr < 10 GeV, the T{2S) 
cross section drops below the T(35) cross section because the T{2S) has a large negative color-octet 
matrix element. (See Table fT3]') The short-distance coefficients multiplying mJ^^^^ are significant at 

n T{nS) 

low pt- Thus, there is a large cancellation between the octet matrix element and , which 

reduces the T(25) cross section in this region, causing it to drop below the T{3S) cross section at low 
Pt- At the high-pj' end of the spectrum, the large value of the ^Si T{2S) color-octet matrix element 
(Table fT3l causes the T(25) cross section to approach that of the T(IS'). In this region, the color-octet 
^5*1 contribution dominates the other channels. Its large matrix element gives the T{2S) an unreasonably 
large cross section relative to that of the T(15). The T(25) rate at px ~ 15 GeV is more reasonable 

n Tf25') 

because the large and positive contribution, and the large and negative contribution nearly 

cancel each other. 

Better determinations of the T matrix elements are required in order to make more accurate pre- 
dictions of the NRQCD T-production rates at the LHC. As is shown in the T-polarization analysis in 
Ref. I5TI . some theoretical predictions have quite large uncertainties even at Tevatron energies, owing 
to our poor knowledge of the matrix elements. However, as is shown in Fig. the ratio Rab is still a 
good measure of the effect of shadowing on T production. The ratio is independent of the T state and is 
quite similar to the J/ip ratio in Fig.[^ The shadowing effect in Pb+Pb interactions may be somewhat 
less for the T at 5 GeV than for the J ftp, but the difference is small. Note also, from Fig. flSl b'). 

that Rab is essentially independent of the matrix elements and is, therefore, largely unaffected by their 
uncertainties. 
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Fig. 18: The pr dependence of Rab [Eq. <10H for T production, (a) We compare the results for pPb and Pb+Pb collisions. 
(The pt dependence is stronger in the Pb+Pb result.) (b) We show the dependence of Rab on the various production channels 
in Pb+Pb collisions at y/s = 5.5 TeV. 

3.4 Comparison of CEM and NRQCD Results ^ 

Here we briefly compare the pT distributions of inclusive J /ip and T production in Pb+Pb collisions 
at 5.5 TeV calculated in the CEM and NRQCD approaches. Neither calculation includes any intrinsic 
transverse momentum effects which could alter the slopes of the pT distributions. 

The J/ ifj distributions are compared on the left-hand side of Fig. ^] The NRQCD result from 
Fig- El is given in the solid curve. The branching ratio to lepton pairs has been removed and the cross 
section converted to jJo. The CEM results from Fig. El are shown in the histograms. In this case, the 
direct J/ip cross section has been converted to the inclusive cross section. Considering the difference 

in mass, scale and parton densities in the two approaches, rric = 1.5 GeV and ^ = ^ 4m^ + p^ with 
MRST98LO for NRQCD and the parameters in Table |5l for the CEM, the agreement is rather good 
over the pr range shown. However, the px slopes appear to be somewhat different. Note also that the 
NRQCD calculations are made in the rapidity interval |y| < 1 while the CEM results are integrated over 
all rapidity, affecting the relative normalization. We can expect a comparable level of agreement for the 
other charmonium states which have a px dependence similar to that of the inclusive J/ ijj in Fig. El 

The inclusive T(IS') distributions are compared on the right-hand side of Fig. El Here also we 
have converted the NRQCD result to /ib and divided out the branching ratio to lepton pairs. We again 
see a relatively good agreement of the calculations in the two approaches, despite some differences in 
the masses, scales and parton densities used, as in the case of the J /iIj. The rapidity \y\ < 1 cut on the 
NRQCD result is a smaller relative factor for T production due to the narrower T rapidity distribution, 
see Fig. We note that the agreement of the two approaches for the T{2S) and T{3S) states would 
not be as good, primarily due to the poorly determined matrix elements for these states, as previously 
discussed. 

''Author: R. Vogt. 
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Fig. 19: The pr dependence of the NRQCD and CEM results for inclusive J/ip (left) and T (right) production in 5.5 TeV 
Pb+Pb collisions is compared. The curves are the NRQCD predictions while the histograms are the CEM predictions. On the 
left-hand side, the J/ip CEM predictions are V*! (solid), ■ip2 (dashed), V'3 (dot-dashed) and V'4 (dotted). On the right-hand side 
the CEM T calculations are Tl (solid), T2 (dashed), T3 (dot-dashed) and T4 (dotted). 



3.5 The Comover Enhancement Scenario (CES) ^ 

3.51 The Quarkonium Thermometer 

The production of heavy QQ quarkonia may offer valuable insights into QCD dynamics, complemen- 
tary to those given by open heavy flavor production. In both cases, the creation of the heavy quark pair 
requires an initial parton collision of hardness Oimq). Most of the time the heavy quarks hadronize 
independently of each other and are incorporated into separate hadrons. The QCD factorization theorem 
exploits the conservation of probability in the hadronization process to express the total heavy quark pro- 
duction cross section in terms of target and projectile parton distributions and a perturbative subprocess 
cross section such as a{gg QQ). 

The quarkonium cross section is a small fraction of the open flavor one and is thus not constrained 
by the standard QCD factorization theorems. Nevertheless, it is plausible that the initial QQ production 
is governed by the usual parton distributions and hard subprocess cross sections with the invariant mass of 
the QQ pair constrained to be close to threshold. Before the quarkonium emerges in the final state there 
can, however, be further interactions which, due to the relatively low binding energy, can either "make 
or break" the bound state. Quarkonium studies can thus give new information about the environment of 
hard production, from the creation of the heavy quark pair until its "freeze-out". The quantum numbers 
of the quarkonium state furthermore impose restrictions on its interactions. Thus states with negative 
charge conjugation, '4){nS), or total spin J = 1, Xci^ require the QQ pair to interact at least once after 
its creation via gg QQ. 

Despite an impressive amount of data on the production of several quarkonium states with a va- 
riety of beams and targets we still have a poor understanding of the underlying QCD dynamics. Thus 
quarkonia cannot yet live up to their potential as 'thermometers' of AB collisions, where A,B = ^^*\ 
hadron or nucleus. Rather, it appears that we need simultaneous studies and comparisons of several 

processes to gain insight into the production dynamics. 

^Authors: P. Hoyer, N. Marchal and S. Peigne. 
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We will now summarize the successes and failures of the Color Singlet Model f79ll80l . which we 
consider as a guideline for understanding the nature of the quarkonium production dynamics. 

3.52 Successes and Failures of the Color Singlet Model 

In the Color Singlet Model (CSM), the QQ pair is directly prepared with the proper quantum numbers in 
the initial hard subprocess and further interactions are assumed to be absent. The quarkonium production 
amplitude is then given by the overlap of the non-relativistic wave function with that of the QQ pair. 

This model at NLO correctly predicts the normalization and momentum dependence of the J/ijj 
photoproduction rate BHll^ . While the absolute normalization of the CSM prediction is uncertain by 
a factor of 2 — 3 there appears to be no need for any additional production mechanism for longitudinal 
momentum fractions 0.3 < xp < 0.9 and 1 < p|i < 60 GeV^. The comparison with leptoproduction 
data [ 82.1 is less conclusive since only LO CSM calculations exist. 

The CSM underestimates the directly produced J/ip and xp' hadroproduction rates by more than an 
order of magnitude. This is true both at low p± < rric (fixed target) l83l and at high p±_ » nic (collider) 
(40]. Similar discrepancies for the T states 1281 l84l l85l indicate that the anomalous enhancement does 
not decrease quickly with increasing quark mass. 

The inelastic cross section ratio ct(V'') /o'dir('^/V') is similar in photoproduction l86l and hadropro- 
duction j^lHHl and consistent with the value ~ 0.24 expected in the CSM l89l . The ratio does not 
depend on Xi;' in the projectile fragmentation region and is independent of the nuclear target size in hA 
collisions. The CSM thus underestimates the J/t/j and tp' hadroproduction cross sections, as well as 
that of the Xci 1 89 1, by similar large factors. The quantum numbers of these charmonium states require 
final-state gluon emission in the CSM, gg — > J/ tpg. This emission is not required for the Xc2 where the 
CSM cross section (j{gg — > Xc2) is only a factor ~ 2 below the hadroproduction data l89l . 

In the CSM, Xc photoproduction is suppressed by a power of as compared to the J ftp and i/;' pro- 
duction rates. One indeed observes a smaller value of the a{xc2)/(^{J/4') ratio in photoproduction ll90l 
than in hadroproduction 1.91. .92.1 . 

3.53 Description of the CES and its generic timescales 

The analysis of agreements and discrepancies between the CSM and quarkonium data led to the comover 
enhancement scenario of quarkonium production f 93 941. Hadroproduced QQ pairs are created within a 
comoving color field and form J/ip, ijj' and Xci through gluon absorption rather than emission, enhanc- 
ing the cross section relative to the CSM since the pair gains rather than loses energy and momentum. 
The Xc2 cross section is not as strongly influenced since no gluon needs to be absorbed or emitted. Most 
importantly, such a mechanism is consistent with the success of the CSM in photoproduction since no 
color fields are expected in the photon fragmentation region, xp ^ 0.3. 

The origin of a comoving color field in hadroproduction is illustrated in Fig. |20l Light charged 
particles carry gauge fields which are radiated in high energy annihilations into a heavy particle pair. 
In e+e~ ^^Ai annihilations, the photon fields pass through each other and materialize as forward 
bremsstrahlung. Fig. l20l' a). In gg QQ, on the other hand, the self-interaction of the color field can 
also result in the creation of a gluon field at intermediate rapidities, Fig. l20l b). Hadroproduced QQ pairs 
thus find themselves surrounded by a color field. We postulate that interactions between the QQ pair 
and this comoving field are important in quarkonium hadroproduction. In direct photoproduction, the 
incoming photon does not carry any color field and the QQ pair is left in a field-free environment after 
the collision, Fig. l20l c). The proposed rescattering thus does not affect (non-resolved) photoproduction. 

The importance of rescattering effects in hadroproduction as compared to photoproduction is also 
suggested by data on open charm production. Hadroproduced DD pairs are almost uncorrelated in 
azimuthal angle |95|, at odds with standard QCD descriptions. Photoproduced pairs on the other hand, 
emerge nearly back-to-back |96|, following the charm quarks of the underlying ^g cc process. 
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(b) QCD (hadroproduction) 



Q 
Q 



A 



(c) QCD (photoproduction) 

Fig. 20: Schematic scenarios of gauge field interactions are compared: (a) e^e~ in QED; (b) hadroproduction of a 

QQ pair, e.g. gg QQ; and (c) photoproduction of a QQ pair, 717 ^ QQ. The creation of a comoving color field is specific 
to hadroproduction, (b). 



Since the hardness of the gluons radiated in the creation process increases with quark mass, the 
rescattering effect persists for bottomonium. Due to the short timescale of the radiation the heavy quark 
pair remains in a compact configuration during rescattering and overlaps with the quarkonium wave 
function at the origin. The successful CSM result for a{i/j')/adir{J/ip) Ii89i is thus preserved. 

The QQ pair may also interact with the more distant projectile spectators after it has expanded and 
formed quarkonium. Such spectator interactions are more frequent for nuclear projectiles and can cause 
the breakup (absorption) of the bound state. This conventional mechanism of quarkonium suppression 
in nuclei is thus fully compatible with, but distinct from, interactions with the comoving color field. 

We have investigated the consequences of the CES using pQCD to calculate the interaction be- 
tween the QQ and the comoving field. While we find consistency with data, quantitative predictions 
depend on the structure of the comoving field. Hence tests of the CES must rely on its generic features, 
especially the proper timescales over which the QQ proceeds from production through hadronization. 

The CES distinguishes three proper timescales in quarkonium production: 

• tq ~ l/niQ, the QQ pair production time; 

• TAP, the DGLAP scale over which the comoving field is created and interacts with the QQ pair; 

• TA ~ 1/Aqcd> while rescattering with comoving spectators may occur. 

In the following we will consider quarkonium production at p_\_ < mq. In quarkonium production at 
Pt ^ Tnq, a large px parton is first created on a timescale l/pr, typically through gg —>■ gg*. The 
virtual gluon then fragments, g* — > QQ, in proper time tq. Thus high px quarkonium production is also 
describable with the CES |97j. 



Timescale tq ~ I/tjiq: creation of the QQ pair — The QQ pair is created in a standard parton 
subprocess, typically gg —>■ QQ, at a time scale tq ~ l/mg. This first stage is common to other 
theoretical approaches such as the CEM l981l99lfT00ll and NRQCD OHj^. The momentum distribution 
of the QQ is determined by the product of projectile (A) and target {B) parton distributions, such as 
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9A{xi)gB{x2) where A, B may be a hadron, nucleus or resolved real or virtual photon. (In the case 
of direct photoproduction or DIS, when the projectile is a photon or lepton, the process depends on the 
single distribution 5b(x2).) Such production is consistent with the quarkonium data. 

According to pQCD, the QQ is dominantly produced close to threshold in a color octet, S = L = 
configuration. Such a state can obtain the quarkonium quantum numbers through a further interaction 
which flips a heavy quark spin and turns the pair into a color singlet. The amplitude for processes of this 
type are suppressed by the factor k / uiq where k is the momentum scale of the interaction. The various 
theoretical approaches differ in the scale assumed for k. 
CSM: Here k = 0{mQ). Thus J /if) production proceeds via the emission of a hard gluon in the 

primary process, gg QQ + g. The Xc2 is produced without gluon emission, gg Xc2, through 

a subdominant S = L = 1 color singlet production amplitude. 
NRQCD: The QQ quantum numbers are changed via gluon emission at the bound state momentum 

scale k = 0{asmQ). This corresponds to an expansion in powers of the bound state velocity 

V = k/niQ, introducing nonperturbative matrix elements that are fit to data. 
CEM: Here k = O(Aqcd)- Soft interactions are postulated to change the QQ quantum numbers with 

probabilities that are specific for each quarkonium state but independent of kinematics, projectile 

and target. 

CES: The quantum numbers of the QQ are changed in perturbative interactions with a comoving field 
at scale k = 0{\/tap), as described below. 

Timescale tap- interactions with the comoving field — The scale tap refers to the time in which 
collinear bremsstrahlung, the source of QCD scaling violations, is emitted in the heavy quark creation 
process ll93l . Thus \/tap characterizes the effective hardness of logarithmic integrals of the type 
dk/k where ij.f "mq is the factorization scale. We stress that 1 / tap is an intermediate but 
still perturbative scale, tq <C tap <C ta, which grows with itiq. 

The fact that the QQ pair acquires the quarkonium quantum numbers over the perturbative time- 
scale TAP is a feature of the CES and distinguishes it from other approaches. At this time, the pair is still 
compact and couples to quarkonia via the bound state wavefunction at the origin or its derivative(s). Thus 
no new parameters are introduced in this transition. However, the interactions of the QQ pair depend on 
the properties of the comoving color field such as its intensity and polarization. Quantitative predictions 
in the CES are only possible when the dependence on the comoving field is weak. 

Ratios of radially excited quarkonia, such as (T(V'')/^dir(<^/V')» are insensitive to the comoving 
field and are thus expected to be process-independent when absorption on spectators at later times can 
be ignored, see below. The fact that this ratio is observed to be roughly universal 1 87 , 88J is one of the 
main motivations for the CES. Even the measured variations of the ratio in different reactions agree with 
expectations, see Ref. IIOIII for a discussion of its systematics in elastic and inelastic photoproduction, 
leptoproduction and hadroproduction at low and high pT- 

The ratio cr(xci)/o"(Xc2) is measured to be 0.6 it 0.3 in pion-induced 1 10211 and 0.31 ± 0.14 in 
proton-induced 1103] reactions. The CSM underestimates this ratio as well as that of cf{J /il:) / a{x2) 
l89l . The rescattering contribution increases a{J/ip) and ct(xi), enhancing the above ratios. 

Nuclear target dependence — The quarkonium cross section can be influenced by rescattering effects 
in both the target and projectile fragmentation regions. For definiteness, we assume the charmonium is 
produced in the projectile fragmentation region, xj? > 0. 

The nuclear target dependence is usually parameterized as a{hA — > J /ij; + X) (x A". Data 
show that a < 1 and obeys Feynman scaling: a depends on (and decreases with) xp rather than on the 
momentum fraction X2 of the target parton II1041I1051 . The comparison with lepton pair production in the 
Drell-Yan process shows that the J /ip nuclear suppression cannot be attributed to shadowing of parton 
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distributions in the nucleus |88|. The A dependence is thus difficult to explain in the CSM, NRQCD and 
CEM approaches. 

In the Feynman scaling regime, we may assume that the QQ pair energy is high enough to remain 
compact while traversing the target. The relative transverse momentum of the Q and Q could increase as 
a result of rescattering in the target, thus suppressing the binding probability. However, this explanation 
is unlikely in view of the absence of nuclear suppression in photoproduction 1 106 1 . 

In the CES, the nuclear target suppression is ascribable to absorption of the comoving color field 
in the target nucleus. This field is emitted by a projectile parton with transverse size tap, larger than the 
size, ~ l/niQ, of the QQ pair. Due to Lorentz time dilation, the field is emitted long before reaching 
the target and reinteracts with the QQ long after passing the target. Absorption of the comoving field in 
the target implies suppression of J /ip production in the CES. At high energies, we have xi ~ xp, which 
explains the Feynman scaling of this effect. Moreover, xp increases, less energy is available to be 
radiated to the gluon field which therefore becomes softer, further increasing its absorption in the target 
and thus the nuclear suppression. 

This explanation is consistent with the fact that the nuclear target suppression of the J/V' and the 
ijj' is found to be the same for xp > 0.2 I105II . It also predicts that the suppression will be similar for 
Xci production. On the other hand, the nuclear target suppression should be reduced for Xc2 since a sub- 
stantial fraction is directly produced without gluon absorption from the comoving field. A measurement 
of (T{hA Xci,Xc2 + X) in the projectile fragmentation region would thus constitute an important test 
of the CES. 

In a Glauber picture of the nuclear suppression, a relatively large value for the absorption cross 
section is required, (Tabs ~ 5 mb. We interpret this value as the joint cross section of the QQ pair and 
the comoving field, thus of order t\p » l/w-g- Since 1/tap scales with mq, we expect less nuclear 
absorption for the T states than for the J/^p, as observed experimentally L107..108J . 

Timescale ta ~ 1/ Aqcd : interactions with comoving spectators — By the time, ta ~ 1 fm, that 
the QQ pair encounters comoving projectile spectators, the pair has already expanded and is distributed 
according to the quarkonium wave function. The spectator rescattering effects are thus independent of 
the quarkonium formation process. Larger and more loosely bound charmonia are more easily broken 
up by secondary scattering. Hence the Xc and tf;' cross sections should be depleted compared to that of 
the J/V'- Likewise bottomonium is generally less affected by spectator interactions than charmonium. 

Spectator interactions at large enough xp are likely to be unimportant for hadron projectiles judg- 
ing from the approximate universality of the a{'ip')/adir{J/'ip) ratio in photo- and hadroproduction. The 
lower ratio seen in nucleus-nucleus collisions [ 109 1, on the other hand, is most naturally explained by ab- 
sorption on spectators. It would be important to confirm this by also measuring the ratio in hA scattering 
in the nuclear fragmentation region. 

3.54 Summary 

Data on quarkonium production have proved challenging for QCD models. The richness of the observed 
phenomena indicates that quarkonium cross sections aie indeed sensitive to the environment of the hard 
QCD scattering. We may, however, only be able to decipher its message through systematic experimental 
and theoretical studies of several species of quarkonia produced with a variety of beams and targets in a 
range of kinematic conditions. 

The rates of quarkonium production at the LHC are large enough for high statistics studies of the 
production mechanism. Differences in pp rates and pT distributions may help distinguish between the 
CEM, NRQCD and CES approaches. 

Acknowledgments S. R would like to thank Cristina Volpe for very stimulating discussions. 
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4. ABSORPTION AND ENERGY LOSS ^ 

In this section, we first discuss J ftp absorption in matter by nucleons (sect. I4.lt and comoving secon- 
daries (sect. l4.2t . We then present a brief review of heavy quark energy loss (sect. 14.31) . 

The effective nucleon absorption cross sections in section 14.11 are obtained by extrapolating the 
values extracted from the fixed-target measurements to higher energies. The resulting charmonium sup- 
pression can be used as a baseline in the search for deconfinement. 

In section l4!2l an evaluation of the couplings of the J/ip to Z)(*)Z)(*) and D^*^D^*^tt are presented 
in the Constituent Quark Model. These couplings are a crucial ingredient in the calculation of cross 
sections for the processes ttJ/^/; D^*^D^*\ an important background for the J/ip suppression signal 
in quark-gluon plasma. While final results relevant to the cross section for J/ i/j absorption by comovers 
are not yet available, the finished calculation should provide some of the most reliable results for these 
processes. 

Energy loss has been much discussed in the context of jet quenching. While some calculations 
have addressed quarkonium, this aspect of energy loss is still too unsettled to be summarized. Therefore, 
only energy loss by heavy quarks in medium is evaluated in section |431 

4.1 Charmonium Suppression in Nuclear Collisions at the LHC: the baseline ^ 

Suppression of the J/ in nucleus-nucleus collisions was proposed as a signal of color deconfinement 
a long time ago Q. Shortly after, such a suppression was observed in 0-i-U and, subsequently, in S-i-U 
collisions I1101llll1lll2l . It was promptly noted that at least part of the observed effect was already 
present in pA collisions Fl 131 11 14l . A quantitative analysis of charmonium suppression in nuclear matter 
has been performed in Ref. fll51. Here we extend that analysis to RHIC and LHC energies, ^/s = 200 
GeV and 5.5 TeV respectively. In this extension, we assume that the increased energy does not change 
the nature of the nuclear absorption process. In principle, interference between successive scatterings 
could lead to modifications at higher energies. The forthcoming p{d)A data from RHIC could clarify 
this point. 

Apart from the independently-determined nuclear matter distributions, the only parameter entering 
the calculation is the effective absorption cross section of the quarkonium precursor in nuclear matter, 
determined from pA data. The latest data indicate cjabs =^ 5 mb 11161 11171 . We assume that (Jabs 
increases with charmonium-nucleon center of mass energy as Sj/^^, reflecting the growth of the small 

X gluon density in nucleons. A massive QQ state produced at zero rapidity will interact with nucleons 
at center of mass energy squared sjj^j^ ~ 2mj/^(\/s/2) where y/s is the center of mass energy per 
nucleon pair of the heavy ion collision. Therefore, the effective absorption cross section will increase 
with ^/s as 



We choose ~ 17.3 GeV (the CERN SPS energy), cTabs(\/so) = 5 ± 0.5 mb for charmonium, 
and A = 0.125. Note that since cjabs is the absorption cross section of the precursor charmonium 
state, it is assumed to be the same for all charmonium resonances. With these parameters, the values of 
the absorption cross sections at RHIC and LHC are given in Table The energy dependence of the 
absorption cross section is shown in Fig. 1211 

^Section coordinator: D. Kharzeev. 
'Authors: D. Kharzeev, M. Nardi and H. Satz. 




(13) 
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Table 14: Effective charmonium-nucleon absorption cross sections at SPS, RHIC, and LHC energies. The middle row is the 
result for the central value of (7abs(\/so) while the upper and lower rows show the errors on crabs- 



0-abs (mb) 


SPS 


RHIC 


LHC 


4.5 


6.1 


9.3 


5.0 


6.8 


10.3 


5.5 


7.5 


11.3 




The quarkonium survival probability in an AB collision is given by, see Eq. (6) of Ref. 11151 . 

Sj/ipib) = J(fsJdzAdzBPA{s,ZA)pBib-s,ZB) 

X exp |-(^ - 1) J dz'^pA{s, z'A)cFah^ 

xexp|-(5-l) j dz'^pB{h - s,z'Q)a3i,^ . 

The resulting J /ip survival probability as a function of impact parameter in a Pb+Pb collision is given 
in Table [13] for SPS, RHIC and LHC energies. The LHC result is also shown in Fig.l^ 

To illustrate the A dependence, we have also calculated the survival probability in Ar+Ar collisions 
at the LHC. The results are shown in Table [T^ and Fig.^H 

Since the impact parameter of the collision cannot be measured directly, it is convenient to choose 
the number of participants, those nucleons that have participated in at least one inelastic collision, as a 
measure of centrality. At least in principle the number of participants, A'part, can be inferred from the 
measurements of the number of spectator nucleons, A'^gpoct, going forward since A^part = A — A^spoct- 
The survival probabilities in Pb+Pb and Ar+Ar collisions at the LHC are shown in Fig.|23]as a function 
of N^.rt/N^f where N^^ = 2A. 

Figure|24]the compares charmonium survival probabilities in Pb+Pb as a function of A^part at SPS, 
RHIC, and the LHC. The suppression due to nuclear absorption alone at the LHC is expected to be much 
stronger than at the SPS. 

Finally we note that at the LHC, bottomonium production will be abundant enough for study of 
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Table 15: J/V* survival probability in Pb+Pb collisions as a function of impact parameter for SPS, RHIC, and LHC energies. 
The central results are calculated with the value of (Tabs given in the middle row of Table fT4l The errors represent the uncertainty 

in CTabs- 







J /yj 




6(fm) 


SPS 


RHIC 


LHC 


(b) 


0.52 ± 0.03 


0.42 ± 0.03 


0.29 ± 0.03 





0.47 ± 0.03 


0.37 ± 0.04 


0.25 ± 0.03 


1 


0.47 ± 0.03 


0.37 ± 0.04 


0.25 ± 0.03 


2 


47 ± 03 


37 ± 04 


25 ± 03 


3 


0.48 ± 0.03 


0.38 ± 0.04 


0.25 ± 0.03 


4 


0.48 ± 0.03 


0.38 ± 0.04 


0.25 ± 0.03 


5 


0.49 ± 0.03 


0.39 ± 0.04 


0.26 ± 0.03 


6 


0.50 ± 0.03 


0.40 ± 0.04 


0.27 ± 0.03 


7 


0.52 ± 0.03 


0.42 ± 0.04 


0.29 ± 0.03 


8 


0.54 ± 0.03 


0.44 ± 0.03 


0.31 ± 0.03 


9 


0.57 ± 0.03 


0.47 ± 0.03 


0.34 ± 0.03 


10 


0.60 ± 0.03 


0.51 ±0.03 


0.38 ± 0.03 


11 


0.65 ± 0.03 


0.57 ± 0.03 


0.44 ± 0.03 


12 


0.71 ± 0.02 


0.63 ± 0.03 


0.51 ±0.03 



their suppression in hot matter to be feasible. Because of their smaller size and larger binding energy, the 
T states will provide valuable information complementary to the study of charmonium. The available 
E772 data on T production at = 39 GeV 1 107 118| imply an effective 66 absorption cross section 
of (Tabs = 2.5 ± 0.5 mb. Assuming the energy dependence of this effective absorption cross section is 
the same as that of Eq. (IT3l . we find Cabs = 4.6 ± 0.9 mb at 5.5 TeV. The corresponding T survival 
probabihties are presented in Fig. 125) 

4.2 A Note on the Jj^ Strong Couplings '° 

In this section, we present a preliminary report on a study of J/V' absorption due to its interaction 
with the hot hadronic medium formed in relativistic heavy-ion collisions. We will give the full analysis 

'"Authors: A. Deandrea, G. NarduUi, A. D. Polosa. 



Table 16: The J/t/j survival probability in Ar+Ar collisions at the LHC as a function of impact parameter. The central results 
are calculated with the value of tJabs given in the middle row of Table fT?! The errors represent the uncertainty in (Tabs- 



6(fm) 






0.52 ± 0.03 





0.45 ± 0.03 
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0.45 ± 0.03 
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0.46 ± 0.03 


3 


0.48 ± 0.03 


4 


0.51 ±0.03 


5 


0.56 ± 0.03 


6 


0.62 ± 0.03 
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0.69 ± 0.02 
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0.77 ± 0.02 



41 



o 

X 



£1 

o 



o 

C/3 

to 
to 
o 




10 12 



Fig. 22: The impact parameter dependence of Sj/^ in Pb+Pb and Ar+Ar collisions at the LHC. We assume = 5.5 TeV for 
both systems. In both cases, the survival probabilities in the central curves are calculated with the median aabs extrapolated to 
LHC energies while the upper and lower curves give the uncertainty m Sj/^ due to the absorption cross section. 



elsewhere 11191 . Here we concentrate on the strong couplings of the J/ijj to low-mass charm mesons 
and pions. In the calculation of the relevant absorption cross sections, tree-level diagrams such as those 
depicted in Fig.|26lare encountered. Previous studies of these effects can be found in 1 120"121"122"1231. 
Besides the DD*n couplings, for which both theoretical 11 24. 125, 1261 and experimental [127| results 
are available, in Fig.|26lthe JD^*^D^*^ and JD^*^D^*^tt couplings also appear (in this section, we shall 
often use J to indicate the J/ip). These couplings have been estimated by different methods, that are, 
in our opinion, unsatisfactory. For example, the use of SU(4) symmetry puts the charm quark and the 
light quarks on the same footing, at odds with the results obtained within the Heavy Quark Effective 
Theory (HQET), which treats the charm as infinitely massive, rUc ^ Aqcd (see Ref. I126II ). Similarly, 
approaches based on Vector Meson Dominance (VMD) should be considered critically given the large 
extrapolation involved, = ^ ^'j/ip- ^ different evaluation, based on QCD Sum Rules, presents the 
typical theoretical uncertainties of this method 1 128 . 129J . Here, we take another approach, based on the 
Constituent Quark Model (CQM), a quark-meson model which explicitly takes into account the HQET 



symmetries. For more details on the CQM see Ref. I130II131I . 

The CQM is particularly suitable for studies of exclusive heavy meson decays. Since its La- 
grangian contains the Feynman rules for the heavy-light vertices formed by a heavy meson, transition 
amplitudes are computable via simple constituent quark loop diagrams where mesons enter as external 
legs. The model is relativistic and incorporates the chiral SU(2) symmetry of the light quark sector as 
well as the heavy quark symmetries. The calculation of the DD*it coupling in the CQM can be found 
in Ref. 1 130'.'T3r|. Here we shall consider the calculation of the JD^*^D^*^ and JD^*'^D^*\ vertices 
within the CQM. This will be done by treating the J / ip with VMD, as depicted in Fig.|^ and computing 
the Feynman diagrams prescribed by the CQM. 

In the CQM, the evaluation of the loop diagram depicted on the left-hand side of Fig. l27l amounts 



(14) 



to the calculation of the Isgur-Wise function I130II131I . The result is 



2 / 2z^ 



i + ijj V 1 + ^ 



where m is the light constituent mass, Zh is a constant arising from the coupling to their constituent 
quarks so that the coupling constant is v^ZhttTd II130II1311 . and the /j integrals are listed in sect. 14.2 11 
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Fig. 23: The centrality dependence of Sj/^, at the LHC in Pb+Pb and Ar+Ar collisions, both at ^/s = 5.5 TeV. In both cases, 
the survival probabilities in the central curves are calculated with the median (Jabs extrapolated to LHC energies while the upper 
and lower curves give the uncertainty in Sj/^, due to the absorption cross section. A^^rt = 2 A is the maximum number of 
participants. 



Here uj = v ■ v' , the product of the two heavy quark velocities, v and v' , assumed to be equal to the D^*'^ 
velocities in the infinite heavy quark mass limit. The hadron-quark mass difference is A// = mo — rric ~ 
0.3 — 0.5 GeV I130II1311 . This form of ^(cj) for the JDD coupling arises from the loop integral 

7 f ^4. Tr [(7 • ^ + ?7t)75(l + 7 • ^07m(1 + 7 • ^)75] 

""levr^i ' - m^){v ■ e + Ah){v' ■ i - Ah) ' ^ ' 

where Nc = 3 and (1 + 7 • f ) /2v ■ k is the heavy quark propagator in HQET. To calculate couplings for 
D* rather than D, the factor (—75) in Eq. (I15t must be replaced by 7 • e where e is the D* polarization. 

The Isgur-Wise function obeys Luke's theorem, ^(1) = 1, arising from the flavor symmetry of 
HQET, corresponding to the Ademollo-Gatto theorem 1132111331 for light flavours. The definition of the 
Isgur-Wise form factor is 



{H{v')\cj^c\H{v)) = -a^)Tr [Hj^H) (16) 
where H is the multiplet containing both the D and the D* mesons 11261 . 

^ = ^^^(-^575 +7-^') , (17) 

where P5 and are charm meson annihilation operators. For example, the transition between two 
pseudoscalar D mesons is 

{D{v')\cj^c\Div)) = mDau;)iv + v')^. (18) 

The Isgur-Wise function can be calculated in the CQM for any value of lo, not only in the region uj > 1, 
which is experimentally accessible via semi-leptonic B D^*^ decays. We note that uj is related to the 
meson momenta by 

^ = P'+P^P\ (19) 
where pi and p2 are the momenta of the two D's. 
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Fig. 24: The dependence of Sj/^ on A^part for Pb+Pb collisions at SPS, RHIC and LHC energies. 



We now consider the right-hand side of Fig. ^\ We use the matrix element 



(0|c7^c| J(g,??)) = fjmj/^t 



(20) 



for the coupHng of the J/tp to the current with fj = 0.405±0.014 GeV. The strong couphngs JD^*^D^*^ 
are given in the following effective lagrangians, 

^jDD = igjDD {DdyD^ .r , 

CjDD* = igjDD'e^'''''^J^d,Dd^Dl, 



^gjD*D* 



D*^ {d^DD r - D*^' (d^Dl) r - (d*^^,di) r] . (21) 



As a consequence of the HQET spin symmetry, we have 



9jd*d* 
gjDD* 

On the other hand, the VMD ansatz gives 

gjDDipl, pI,p^) 



gjDD : 

gjDD 

niD 



m 



J/4> 



(22) 



(23) 



Since gjDD has no zeros, Eq. d23b shows that ^{lo) must have a pole at 



m 



, as expected from 



dispersion relation arguments. The CQM evaluation of ^ does show a strong peak at (2mc 
even though, due to 0{\/mc) effects, the location of the singularity is not exactly at = m 



. This 



is shown in Fig. l28l where we plot gjooipf, Pq,, P^) for on-shell D mesons as a function of using 
Ah = 0.4 GeV and Zh = 2.36/GeV. For < < 4 GeV^, qjdd is almost constant. 



gjDD 



.0±0.5 



(24) 



For larger values of p"^ the method is unreliable due to the incomplete cancellation between the kinemat- 
ical zero and the pole. The distorted shape around the J/iJj pole suggests that the contribution of the 
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Fig. 25: The centrality dependence of the T survival probability at the LHC in Pb+Pb and Ar+Ar collisions. In both cases, the 
survival probabilities in the central curves are calculated with the median a^^s extrapolated to LHC energies while the upper 
and lower curves give the uncertainty in Sj/^, due to the absorption cross section. N^^^t = 2 A is the maximum number of 
participants. 



nearby '4^{2S) pole could also be relevant. Therefore, we extrapolate the smooth behavior of gjoD in 
the small region up to = fn^jj^ and assume the value of gjDD in Eq. (l24l) for on-shell J/ijj. On 
the other hand, the behavior with p\ and is smooth, compatible with that of a smooth form factor. We 
finally note that Eq. (E^ agrees with the result of the QCD sum rule analysis 1 1281 1129'll . This is not 
surprising since QCD sum rules involve a perturbative part and a suppressed nonperturbative contribu- 
tion. The perturbative term has its counterpart in the CQM loop calculation of Fig. |^ with an overall 
normalization that should agree with the CQM as a consequence of Luke's theorem. 

We now turn to the J D^*'^ D'^*'^ tt couplings. As discussed in Refs. F1 26II134TI . the leading contri- 
butions to the current matrix element, {H{v')T^\c'y^c\H{v)), in the soft pion limit (SPL) are the pole dia- 
grams. Technically, in the SPL, the reducing action of a pion derivative in the matrix element is compen- 
sated in the polar diagrams by the effect of the vanishing denominator in the combined 0, rric — > oo 
limit. Since the effect of the pole diagrams is explicitly taken account in Fig. |26l we do not include any 
further contributions. To be definite, we consider the coupling gjoDw which can be obtained by a VMD 
ansatz similar to Fig.|^ Now the left-hand side is modified by the insertion of a soft pion on the light 
quark line with the couphng q!^/ fnlfilb- We call ^'"{uj) the analogous form factor in the SPL, 



h{AH) - + " ' 



(25) 



1 + Lu \ 1 + Lu I dm? 

The VMD ansatz of Fig. |27] gives 

CjDD. = igjDD.e^'''''^J^d,DdaDdpTT (26) 

with 

9JnnAP. P2, P ) = f^f^^^^^ • (27) 

In Fig.|29lwe plot our result for gjoD-n with on-shell D mesons. By the same arguments used to 
determine gjoD^ we find, with all mesons on mass-shell, 

gjDDn = 125 ± 15 GeV"3 . (28) 
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Fig. 26: Feynman diagrams for J/iji absorption by the pion: (1) J/tpiT — » DD, (2) J/ijiTv — » and J/V'tt D* D* . 




Fig. 27: The Vector Meson Dominance equation for the coupUng of the J/ij) with D and D* in terms of the Isgur-Wise 
function ^. The function ^ on the left-hand side is computed by a diagram with a quark loop. The coupling of each D^*^ meson 
to quarks is given by ^/Zhttid. 
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Fig. 28: The dependence of g = gjooirrio, mjj, p^), showing the almost complete cancellation between the pole of £,{uj) 
and the kinematic zero. The units of p^ are GeV'^. 




Fig. 29: The p^ dependence of — gjDDTrimjj, mfj, p^). As in Fig.|28|there is an almost complete cancellation between 
the pole of the form factor and the kinematic zero. The units on p^ are GeV^ and GeV^"' for gjoDir. 



We now compare this result with the effective JDDtt coupUng obtained by a polar diagram with an 
intermediate D* state. In this case, 



polar 
djDDn 



i9JDD'-)igD*D7T) 
2g,r • PD 



(29) 



Since all the calculations are valid in the SPL, only pion momenta of up to a few hundred MeV should 
be considered. Using the result gD*Dn = '^fno/ firQ, with g = 0.59 ib 0.01 ± 0.07 11271 . we find 
9Tdd-k ~ 393, 196, and 98 GeV"^ foj. ^ 59^ lOQ, and 200 MeV respectively. This analysis 
shows that, within the region of validity of the model, in spite of the rather large value of the coupling 
in Eq. (l28l . the diagrams containing this coupling are in general suppressed. Similar conclusions are 
reached for L'*'s. 

Finally, we discuss the kinematical limits of our approach. To allow the production of a D^*^D^*^ 
pair, as shown in Fig. |26l we must extend the region of validity of the model beyond the SPL since the 
D'^*)'D^*^ threshold is = 700 - 1000 MeV. The CQM, along with other models, is a chiral model 
which puts limits on the pion momenta. Therefore it is necessary to include a form factor enhancing the 
small pion momenta region, for example 



fm) 



1 + (I^Trl/m^ 



(30) 



A similar form factor is considered in Ref . 11201 but with a different motivation. Here we introduce it to 
ensure the vaUdity of our approach. In this sense, the cross sections we can compute may be considered 
as lower bounds. Since the main effect of Eq. d30l should be to reduce contributions from pion momenta 
larger than a few hundred MeV, we expect 400 < < 600 MeV. This choice implies that the direct 
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couplings in Fig.|26l diagrams Ic, 2d and 3e, should not dominate the final result since their contribution 
is larger where the form factor is more effective. 



4.21 Definitions of the integrals 

We list the expressions used to compute the integrals Ij used in the text. The ultraviolet cutoff A, the 
infrared cutoff, /i, and the light constituent mass, m, are A = 1.25 GeV, = 0.3 GeV and m = 0.3 GeV 
L13a.l31il . 



^3(A) 



l4(A) 



l5(Ai,A2,u;) 



where 



16^7 (A;2-m2)(ti./fc + A + ie) 
A^^ /-V/^' ds 



167r3/2 7i/a2 §3/2 

16^ 



[1 + erf (A^/i)] 



(31) 



- m'^y{v ■ k + A + ie) 



N, 



V/^' ds 



-s(m2-A2) 



167r3/2 7i/A2 
iNc /""^^s 
16^ 



[1 + erf (A Vi)] 



(32) 



d'^k 



dx- 



- im?){v ■k + Ai + ie){v' ■k + A2 + ie) 
1 



1 + 2x2(1 -a;) + 2x(a;-l 
1/^2 



X 



+ 



6 



167r^/2 J1/A2 

1/^2 



ds (TA e-^(™ -'^a) s-1/2 [1 + erf (c7A\/i)] 



6 



167r2 7i/a: 



ds e-'"^ s-^ 



CTA = cr(x, Ai, A2,w) 



Ai (1 -x) + A2X 



v/1 + 2 (w- 1) x + 2 x2 



(33) 



(34) 



4.3 Heavy Quark Energy Loss in QCD Matter 

The study of heavy quark production in nuclear collisions will allow extraction of valuable information 
on the gluon densities of the colliding nuclei. In addition, since the heavy quarks propagate through 
the quark-gluon matter created in such collisions, they can be used to probe the properties of the dense 
matter. The practical importance of heavy quark energy loss to the charmed hadron and lepton spectra 
has been clearly established 1135lll36lfTT7lfnHira . 

The idea of heavy quark energy loss has attracted considerable attention. In particular, collisional 
energy loss in quark-gluon plasma has been evaluated 1 140, 141. 142 1 using finite-temperature QCD, see 
Ref. [ 143 1 for a comprehensive review. However, the energy loss of fast partons in medium is dominated 
by gluon radiation Ill44lll45lll46ll . Radiative energy loss by heavy quarks was evaluated in Ref. I147II . It 
was shown 1 147 1 that the "dead cone effect" in the radiation of a heavy quark significantly suppresses its 
energy loss while making the calculation more reliable by reducing the sensitivity to the infrared region. 

In fact, the first measurements of charm production in Au+Au collisions at RHIC 11481 11491 
indicated very small, if any, nuclear effect on the shape of charm transverse momentum distributions. 
The smaller loss predicted in Ref. | provides a plausible and economical explanation of the observed 



Author: D. Kharzeev. 
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phenomenon r i50[ll51 :|. The "dead cone" effect has been taken into account in the studies of charm and 
bottom production at the LHC 1 152|. 

In the following we give a brief and qualitative explanation of the results of Ref. 11471 . We first 
recall the basic features of gluon radiation caused by propagation of a fast parton (quark) through a QCD 
medium. 

As pointed out in Ref. fl45^, the accompanying radiation is determined by multiple rescattering 
of the radiated gluon in the medium. The gluon, during its formation time, 

u 

^form — T2" ' (35) 

accumulates a typical transverse momentum 

, 2 2 ^form /T/'N 

- M — ^> (36) 

where A is the mean free path and fi^ the characteristic momentum transfer squared in a single scattering. 
This is the random walk pattern with an average number of scatterings given by tform/ 
Combining Eqs. (l35l and (l36l we obtain 

the number of scattering centres which participate, coherently, in the emission of the gluon with a given 
energy uj. For sufficiently large gluon energies, lo > /i^A, when the coherent length exceeds the mean free 
path, A'^coh > 1. In this situation, the standard Bethe-Heitler energy spectrum per unit length describing 
independent emission of gluons at each scattering centre is suppressed: 



dW 1 / dW asCR hi^X OsCr q 



dujdz Nf-oh \dujdz 



J TTUJ \ y UJ TTUJ V 



Here Cr is the "colour charge" of the parton projectile (Cr = Cp = (A^ — l)/2Nc = 4/3 for quarks). 
In Eq. d38l we have substituted the characteristic ratio /i^/A by the gluon transport coefficient 

9' dq^ , (39) 

proportional to the density p of the scattering centres in the medium. The transport coefficient character- 
izes the typical momentum transfer in gluon scatterings off these centres. 

An important feature of medium-induced radiation is the relation between the transverse momen- 
tum and the energy of the emitted gluon. Indeed, from Eqs. d35l and d36b . see also Eq. ( l39l . we derive 

~ y/J^ . (40) 

Thus the angular distribution of gluons with a given energy u is concentrated at a characteristic energy- 
and medium- dependent emission angle. 



UJ \UJ^ 



(41) 



Gluon bremsstrahlung off a heavy quark differs from the case of a massless parton produced at 
the same hardness scale in one respect: gluon radiation is suppressed at angles smaller than the ratio of 
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the quark mass iriQ to its energy E. Indeed, the distribution of soft gluons radiated by a heavy quark is 
given by 

where the strong coupling constant should be evaluated at the scale determined by the denominator of 
Eq. (l42l) . Equating kx with uj6 in the small-angle approximation, we conclude that Eq. (l42l) differs from 
the standard bremsstrahlung spectrum. 



dP^ ^ __L = ri^ZIl , (43) 



as Cp du) dk?p Ug Cp duj dO 

\irp 



by the factor 



^0 



dPHQ = dPo 1 + ^ . (44) 



This effect is known as the "dead cone" phenomenon. Suppression of small-angle radiation has a number 
of interesting implications, such as the perturbative calculability of, and nonperturbative A/mg correc- 
tions to, heavy quark fragmentation functions 1154111551 . multiplicity and energy spectra of Ught particles 
accompanying hard production of a heavy quark II1561I1571 . 

In the present context we should compare the angular distribution of gluons induced by the quark 
propagation in the medium with the size of the dead cone. To this end, for the sake of a semi-quantitative 
estimate, we substitute the characteristic angle, Eq. ^A\\ . into the dead cone suppression factor, Eq. (144b . 
and combine it with the radiation spectrum, Eq. (l38l . to arrive at 



dW as Cf 1 

I[uj)=uj— = J— , , (45) 

cLlI tt \ lo (1 + w)'^/^)^ 

where 

i^r'^' {^y. (46) 

To determine whether the finite quark mass essentially affects the medium-induced gluon yield, 
we need to estimate the product £uj for the maximum gluon energy, u; ~ lji, achievable in Eq. ( 138b . 

^-1 = r^/^ = (^)''' , ^HQ - m^^^ . (47) 

where L is the path length through the medium. Thus the quark mass becomes irrelevant when its energy 
exceeds the characteristic value, E}\q, dependent on the size of the medium and on its "scattering power" 
embodied in the transport coefficient q. 

Which regime is realized in the experiments on heavy quark production in nuclear collisions? 
Taking rric = 1.5 GeV, we estimate 11471 

/ T \3/2 

E^Q = 20GeV(^^j , (48) 

/ T \ 3/2 

i?HQ = v^L3/2m,=. 92GeV(^^j , (49) 

for cold and hot matter, respectively. We observe that, in practice, E <C £^hq for the transverse mo- 
mentum (energy) distributions of heavy mesons, especially in the hot medium. We thus conclude that 
the pattern of medium-induced gluon radiation appeai^s to be qualitatively different for heavy and light 
quarks in the kinematic region of interest. 
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The issue of in-medium quenching of inclusive particle spectra was addressed in Ref. 1 158 1. The 
Pt spectrum is given by the convolution of the transverse momentum distribution in an elementary 
hadron-hadron collision, evaluated at a shifted value pT + e, with the distribution D{e) in the energy e 
lost by the quark to the medium-induced gluon radiation: 



da 



med 



dpj. 



de D{e) 



da'^'^pT + e) _ dfT™=(pT) 



dpj. 



dp^ 



Q{pt) 



(50) 



T 



where Q{pt) is a medium-dependent quenching factor. Because in the region of interest, e ^ pr, when 
the vacuum cross section is a steeply-falling function, the calculation of the quenching factor Q can be 
simplified by approximating the e-integral in Eq. d50b as an exponential, 



Q{P7 



j deD{e) exp|-^£| , C 



d 



dlnpT 



■In 



dpj: 



(51) 
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1 . 25 
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Fig. 30: The ratio of quenching factors Qh{pt)/Ql{pt) for charm and Ught quarks in cold nuclear matter with q = 0.01 
GeV^ and L = 5 fm. From Ref. IT?7l . 

According to Eq. J49t . we expect a larger quenching ratio for a hot medium. Indeed, the D/ir ratio 
should become significantly enhanced relative to pp collisions. Assuming a fixed path length traversed 
by the quarks in the hot medium, L = 5 fm, we find a factor of ~ 2 enhancement at ~ 5 — 10 GeV. 

The following result for the heavy quark quenching factor was derived in Ref. 11471 : 



Qnip 



•T 



exp 



q ml 



2asCF ^ Xh _^ IOosCf^ 

/tF \I'^Pt 9\/3 [m^+pl 



1/3- 



(52) 



The first term in the exponent in Eq. d52l represents the quenching of the transverse momentum spectrum, 
universal for light and heavy quarks, 



Ql{pt) ^ exp 



LJq — 
/TT y Pt 
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modulo the difference in £ determined by the vacuum px distributions. The second term in Eq. d52t is 
specific to heavy quarks. It has a positive sign so that the suppression of the heavy hadron px distributions 
is always smaller than that of the light hadrons. This is a straightforward consequence of the fact that 
the heavy quark mass suppresses gluon radiation. At very high transverse momenta both terms vanish 
in accord with the QCD factorization theorem which states that medium effects should disappear as 
Pt — > oo. How fast this regime is approached depends, however, on the properties of the medium 
encoded in the transport coefficient q and the path length L. 

H - L Ratio 
\ I 



2 n2 - 




1.2 - 
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H - L Ratio 

1.3r 

\ 

\ 

1.2^ - 




7.5 10 12.5 15 17.5 20 ^^ 

Fig. 31: The ratio of quenching factors Qh{pt) /Ql{pt) for charm and hght quarks in hot matter with q = 0.2 GeY^ (L = 5 
fm upper panel, L — 2fm lower panel). Solid lines correspond to unrestricted gluon radiation, while the dashed lines are based 
on the calculation with the cut uj > 0.5 GeV on gluon energies. From Ref. 11471 . 



Constructing the ratio Qh to Ql, the heavy-to-light quenching should be enhanced by 

,2 \ 1/31 



Qh{pt) 
Ql{pt) 



exp 



9V3 



L 



q rrir 



(53) 



This simple expression provides a reasonably good approximation to the more accurate quantitative 
results presented in Ref. 11471 . 
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As discussed above, the quenching of heavy hadron pj- distributions caused by QCD matter is 
much weaker than for pions because the gluon cloud around the heavy quark is "truncated" by the large 
quark mass in a medium-dependent way. This interesting effect can be illustrated by the transverse 
momentum dependence of the ratio of hadrons originating from heavy and light quark fragmentation 
such as the D/tt ratio in heavy ion collisions. 

Fig- Uni shows the ratio of quenching factors, Eq. d53l . for heavy and light quarks in cold nuclear 
matter, L = 5 fm, relevant to high-pT particle production in pA collisions. A small value of the transport 
coefficient, QcoM — 0.01 GeV^, translates into a 15% enhancement, reduced to « 1% for L = 2 fm. 

We only present the ratio of quenching factors because, as pointed out in Ref. [158|, the absolute 
magnitude of the quenching turns out to be extremely sensitive to gluon radiation in the few-hundred 
MeV energy range. Thus the magnitude cannot be quantitatively predicted without a detailed under- 
standing of the spectral properties of the medium. 

The heavy-to-light ratio, however, proves to be much less sensitive to the infrared region since 
gluon radiation off heavy and light quarks is universal in the x ^ limit, see Eq. d45l . To illustrate this 
point, in Fig. |^ we show the ratio of quenching factors calculated with the gluon energies restricted to 
to > 500 MeV. We see that, in the 5-10 GeV range of pr, the ratio is modified by 20 - 30 %P 

Clearly, detailed calculations have to be performed before a reliable estimate of the magnitude of 
the D/tt enhancement can be presented. Nevertheless, the D/tt and B/tt ratios appear to be extremely 
sensitive to the density of colour charges in QCD matter. Of special interest is the B/D ratio, for which 
these calculations become even less sensitive to the infrared region and thus more stable. 



The quenching factors themselves change (increase) by an order of magnitude when the radiation of gluons with energies 
smaller than 500 MeV is vetoed. 
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5. QUARKONIUM FORMATION FROM UNCORRELATED QUARK-ANTIQUARK PAIRS 
5.1 Introduction 

The goal of this section is to assess the possibihty that quarkonium production rates may be enhanced in 
nucleus-nucleus interactions at the LHC relative to that predicted by extrapolation of processes thought to 
be dominant at lower energy. This enhancement could follow from the effects of incoherent recombina- 
tion mechanisms involving uncorrelated pairs of heavy quarks and antiquarks which result from multiple 
pair production. Two different approaches have been considered: statistical hadronization and kinetic 
formation. Updated predictions relevant to Pb-i-Pb collisions at the LHC are given. 

The utility of heavy quarkonium production rates in nuclear colhsions as a signature of color de- 
confinement was proposed more than 15 years ago |1|. Since one expects that the long-range color 
confining potential will be screened in a deconfined medium, the quark and antiquark constituents of 
bound states will be liberated. As the system expands and cools, these constituents will, in general, 
diffuse away from each other to separations larger than typical hadronic dimensions. When the confin- 
ing potential reappears, a given heavy quark will not be able to "find" its heavy antiquark partner and 
form heavy quarkonium. It must then bind with one of the antiquarks within range at hadronization. 
Since these antiquarks are predominantly the lighter u, d, and s flavors, the final hadronic states will 
preferentially be those with "open" heavy flavor. The result will be a decreased population of heavy 
quarkonium relative to that which would have formed if a region of deconfinement had not been present. 
This scenario as applied to the charm sector is known as J/ 1/; suppression. 

At LHC energy, perturbative QCD estimates predict that hundreds of pairs of charm-anticharm 
quarks will be produced in a central lead-lead collision. This situation provides a "loophole" in the 
Matsui-Satz argument jjj since there will be copious numbers of heavy antiquarks in the interaction 
region with which any given heavy quark may combine. In order for this to happen, however, one must 
invoke a physical situation in which quarkonium states can be formed from all combinations of heavy 
quarks and antiquarks. This of course would be expected to be valid in the case that a space-time region 
of color deconfinement is present, but it is not necessarily limited to this possibility. 

One can make a model-independent estimate of how such a "recombination" mechanism would 
depend on nuclear collision observables. For a given charm quark, the probability V to form a. J/ip is 
proportional to the number of available anticharm quarks relative to the number of light antiquarks. 



Nr N, 



u,a,s '-'i 

In the second step, we have replaced the number of available anticharm quarks by the total number of 
pairs initially produced, assuming that the total number of bound states formed remains a small fraction 
of the total cc production. We normalize the number of light antiquarks by the number of produced 
charged hadrons. Since this probability is generally very small, one can simply multiply by the total 
number of charm quarks, Nc, to obtain the number of J/tjj expected in a given event, 

« -j^ , (55) 

where the use of the initial values Ncc = = N-^ is again justified by the relatively small number of 
bound states formed. 

The essential property of this result is that the growth of Nj/^, quadratic in the total number of 
charm quarks, with energy 059 1 is expected to be much faster than the growth of total particle production 
in heavy ion collisions I160II . Without this quadratic mechanism, J/^ production is typically some small 
energy-independent fraction of total initial charm production l20ll . We thus anticipate that the quadratic 
formation will become dominant at sufficiently high energy. Generic estimates of the significance of 
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this type of formation process can be made [ 161 1. Here we look at specific predictions of two models 
- statistical hadronization and kinetic formation, considered in sects . 15 .2] and 15 . 3 1 respectively - which 
share the above properties, and we update the expectations to LHC energies. 



5.2 Statistical Hadronization 

The statistical hadronization model is motivated by the successful fits of relative abundances of light 
hadrons produced in high energy heavy ion interactions according to a hadron gas in chemical and ther- 
mal equilibrium I162II . Extension of the model to hadrons containing heavy quarks underpredicts the 
observed abundances. This effect may be attributed to the long time scales associated with thermal pro- 
duction and annihilation of heavy quarks. The statistical hadronization model as first formulated for 
charm quarks |163| assumes that the cc pairs produced in the initial hadronic interactions survive until 
their subsequent hadronization, at which time they are distributed into hadrons according to the same 
thermal equilibrium parameters that fit the light hadron abundances. Chemical equilibrium abundances 
are adjusted by a factor 7c which accounts for the non-thermal heavy quark density. One power of this 
factor multiplies a given thermal hadron population for each heavy quark or antiquark contained in the 
hadron. Thus the relative abundance of the J/ ip to that of D mesons, for example, may be enhanced in 
this model. 

The value of 7c is determined by conservation of the heavy quark flavor. For the charm sector, 
the conservation constraint relates the number of initially -produced cc pairs A'^cc to their distribution into 
open and hidden charm hadrons, 

1 2 

Ncc = -^IcNopen + 1c Miidden, (56) 

where A'opcn is the number of hadrons containing one c or c quark and A^hiddon is the number of hadrons 
containing a cc pair. For most applications, A'^hidden (and also multi-charm hadrons) can be neglected 
compared with A'open due to the mass differences. Thus the charm enhancement factor is simply 

7c = (57) 

leading directly to the quadratic dependence of the hidden charm hadron population on A'^cc- One can 
then express the total number of J/V' in terms of the various thermal densities, rii, and the total number 
of cc pairs, A^cc- One factor of system volume V remains implicit here. It is generally replaced by the 
ratio of number to density for total charged hadrons, rich/Nch- Then the number of J/iJj produced obeys 
the generic form anticipated in Eq. d55t . 

WjA^ A^l 

For collider experiments such as those at the LHC and RHIC, relating the corresponding central 
rapidity densities will be more relevant. Since Eq. ( I5HI is homogeneous in the total particle and quark 
pair numbers, it will also be valid if these are replaced by their rapidity densities. To get an order of 
magnitude estimate, we choose a "standard" set of thermal parameters, T = 170 MeV and fiB ^ 0, for 
which the thermal density ratio is approximately 0.5. For a specific normalization, we assume dN^h/dy 
= 2000 for a central collision at the LHC and take the initial charm rapidity density to be dNcc/dy = 25, 
roughly corresponding to Ncc = 200 for central collisions (6 = 0). Using these inputs, one predicts 
dNji^/dy = 0.625, indicating that several J /iJj will form through statistical hadronization in a central 
collision. To put this number in perspective, it is revealing to form the J / ip to A'cc rapidity density ratio, 
0.025 with the same assumptions. For comparison, one expects the corresponding hadronic production 
ratio to be of 0.01. This number would then be significantly reduced if placed in a region of color 
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deconfinement. Thus the efficiency of J /if) formation via statistical hadronization at the LHC is expected 
to be substantial. 

These numbers can be easily adjusted to other charm and charged particle densities using Eq. (I58t . 
Variations of the thermal parameters can also be investigated. For example, if the hadronization tempera- 
ture is decreased to 150 MeV, the prefactor combination of thermal densities increases by approximately 
a factor of two. 

The centrahty dependence is controlled by the behavior of N^c and N^y^. The former should be 
proportional to the nuclear overlap function TAA{b) but is generally recast in terms of the dependence 
on the number of nucleon participants, A'part- The calculation of A'part requires a model calculation 
dependent on the total inelastic cross section, (Tin, as well as the nuclear geometry. We parameterize the 
expected behavior as a power-law Ncc oc A'^part^^'^- However, there will be deviations from this behavior 
for the larger values of din expected at the LHC 11641 . The centrality dependence of A'ch at RHIC is 
also consistent with a power-law with exponent w 1.2 [ 1601 11651 . We will use the same dependence 
for our estimates at the LHC. It is clear that for sufficiently peripheral collisions one will encounter 
situations in which the average number of initially produced cc pairs is of order unity or less. At this 
point, one must revisit the assumptions of the original statistical hadronization model which assumed a 
grand canonical ensemble. The grand canonical approach is valid only when Ncc is large enough for the 
fluctuations about the average value to be negligible. Thus for peripheral collisions, one must recalculate 
the statistical results in the canonical approach where the charm number is exactly conserved, as noted 
in I166II . Charm conservation can be implemented via a correction factor I167L 

Ncz = ]:lcNo,J-^p^^^ + Tc'iVhidden. (59) 

In the limit of large 7cA'^opcn, the ratio of Bessel functions /j approaches unity and the grand canonical 
result is recovered. In the opposite limit when 7ciVopen —>■ 0, the ratio of Bessel functions approaches 
^Tc^opon- In this hmit, the dependence on Ncc in Eq. changes from quadratic to linear. At the LHC 
this effect will not be relevant until one reaches very peripheral events, but at lower energies it can be 
significant over a much larger range of centralities 11681 . 

The results for dNj/^/dy as a function of A'^part at the LHC are shown in Fig. |32l The results 
are shown for three different values of dNcc/dy{b = 0), corresponding to A^cc(O) ~ 200, 150, and 100. 
There is a rapid increase with centrality due to the quadratic dependence of Nj/^ on Ncc- 

It is also interesting to look at these results normalized by A^part» shown in Fig.|33l This ratio also 
increases with centrality, providing a signature for the statistical hadronization process that is less depen- 
dent on Ncc for the overall normalization. The corresponding results when normalized by dNcc/dy are 
shown in Fig.|34] The same general behavior is seen but the increase with centrality is less pronounced 
since the dNcc/dy is assumed to vary with a larger power, A'^part^^^- AH of these ratios are at the percent 
level for central collisions and hence are larger than expected if the total J/ip population were due to 
initial production followed by suppression in a deconfined medium. 

The region of very peripheral collisions deserves some separate comments. First, there is a rise 
at low A'part in both Figs. |33l and |34] due to the onset of corrections from the canonical ensemble treat- 
ment. However, the extremely large values of the ratios as A'^part are an artifact of the decreasing 
interaction volume, F — > 0. This calculation must be cut off before A'^part = 2, i.e. only one interacting 
pair. It is also in this region where one must take into account the remaining J/ tp from initial production. 
Since the survival probability is maximum for very peripheral collisions and the statistical hadronization 
process is least effective in this same region, there will be a crossover in the relative importance of the 
two mechanisms. Some studies have already been performed for this situation at SPS and RHIC energies 
II69II1701[T7B . 

Finally, there is another lower cutoff in centrality for the statistical hadronization results, needed 
to avoid a contradiction with the ipi /{J/ifj) ratio at the SPS. Since both of these states receive identical 
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Fig. 32: Statistical hadronization results for J/tp production as a function of A'^part at the LHC. 



factors of 7c, their ratio must be that predicted for chemical equihbrium in the absence of any charm 
enhancement or suppression. Although the measured ratio appears to be consistent for more central 
collisions [1721, there is an indication that it rises sharply for more peripheral collisions. Most treatments 
have thus inserted a cutoff of A'^part = 100, below which model predictions become inconsistent 1 163 1. 

The numerical values for dNj/^/dy are tabulated as a function of impact parameter in Table \T7\ 
for the three choices of initial charm multiplicity density and the default values of all other quantities. 



5.3 Kinetic Formation in a Deconfined Region 

The kinetic model has been developed 1173111741 to investigate the possibility that J/ip may form directly 
in a deconfined medium. This formation takes advantage of the mobility of the initially-produced charm 
quarks in a deconfined region. In order to motivate this view, consider the "standard" physical picture of 
deconfinement in which quarkonium is suppressed by collisions with free gluons in the medium [175 1. 
Then the formation process, in which a c and c in a relative color octet state are captured into a color- 
singlet quarkonium bound state and emit a color octet gluon, is simply the inverse of the the breakup 
reaction responsible for the suppression. This is an inevitable consequence of the suppression picture. 

The proper time evolution of the J/ip population is given by the rate equation 

^ = AF-^-ADiV,/^p,, (60) 

where pg is the gluon number density and V{t) is the time-dependent volume of the deconfined spatial 
region. The reactivities Af,d are the reaction rates, {av^ei), averaged over the momentum distributions 
of the initial participants, i.e. c and c for Af and J /ijj and g for Ad- 

The solution of Eq. (l60l grows quadratically with Ncc, as long as A^j/^ ^ Ncc- In this case, we 
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Fig. 33: The statistical hadronization results of Fig. l32l divided by TVpart as a function of A^part- 



have 

Nj/4V) = <y) K/v(^o) + Nl r drXp [V{t) 6(r) 



TO 



(61) 



The function e(T/) = exp(— /^^ drAD Pg) would be the suppression factor if formation were neglected. 

The quadratic factor N^-^ is present, as expected, for the additional formation process. The normal- 
ization factor of Nch is not immediately evident, but is implicit in the system volume factor. This volume 
is now time-dependent, accounting for the decreasing charm quark density during expansion. Here the 
transverse area of the deconfined region is determined not just by the nuclear geometry but by the dy- 
namics which determine the extent of the deconfined region. This area is modeled by the energy density 
in terms of the local participant density in the transverse plane, npart(&i s = 0). The transverse area is 
defined by the ratio of the participant number to the local participant density. Note that the maximum 
local density is at s = 0. Thus, 

Arib) = AT{0)[Np^rt{b)np^rt{0, s = 0)/Np^,t{0)np^rt{b, s = 0)] (62) 

These area effects will be more explicit when the centrality dependence is considered. 

The numerical results depend on a number of parameters, including the initial volume and temper- 
ature, the time expansion profile, the reaction cross sections, the behavior of the quarkonium masses and 
binding energies in the deconfined region, and the charm quark momentum distributions. For specifics, 
see Ref. 1 176 1. Our previous results have used initial values Ncc = 200, 150, and 100, spanning a rea- 
sonable range of expectations 1 177 1. The results are very sensitive to the initial charm quark momentum 
distributions, as may be expected. We assume the charm px distributions are Gaussian and the charm 
rapidity distributions are flat over a plateau of variable width. Ay. The range 1 < Ay < 7 spans the 
range between an approximate thermal momentum distribution. Ay 1, to a distribution similar to that 
of the initial pQCD production. Ay « 7. 
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Fig. 34: The statistical hadronization results for J/?/> production at the LHC, divided by the open charm multiplicity, dN^-/dy, 
as a function of A'^part • 



The results as a function of the initial number of cc pairs produced in central collisions are shown 
in Fig.|35l There is a rapid decrease in formation with increasing Ay. The quadratic dependence on Ncc 
is evident, but there is also a substantial linear component in some of the curves. This linear contribution 
arises because the final J/ ifj formation by this mechanism is large enough for exact charm conservation 
to reduce the number of c and c quarks available to participate in the formation process. The curve 
labeled "Quadratic Extrapolation" uses a quadratic dependence derived from a fit valid only for low Ncc- 
Note that the result for a thermal distribution is very similar to the assumption Ay = 1. 

The corresponding centrality dependence is presented in Fig.|551 where we give Nj at hadroniza- 
tion for three different initial charm quark momentum distributions, thermal. Ay = 4 and Ay = 7, as 
well as for our three choices of N(^{b = 0). 

Finally, the ratio of final J/t/j to initial charm production is shown in Fig. lousing the same 
parameters as in Fig.|55] These ratios are most easily compared to either initial production or suppression. 
There is a substantial variation in the predictions and it is evident that a simultaneous measurement of 
open charm will be required for an interpretation. However, the centrality dependence is opposite to that 
expected in any pure suppression scenario. 

We have updated the calculations to include the charm quark momentum distribution from a lead- 
ing order pQCD calculation [178 1. The rapidity distribution has a somewhat larger effective Ay and the 
Pt distribution does not fall as fast as a simple Gaussian. As a result, the formation efficiency is further 
reduced. Such distributions may be most relevant, given preliminary results from RHIC 1179111801 . 

The numerical values for Nj are compared with the statistical hadronization model results for 
dNji^/dy in Table ITT] The overall magnitudes are comparable, although the centrality dependences 
differ somewhat. Thus details such as the resulting J /ifj momentum distributions will be required to 
differentiate between these two models 11781 . For completeness, A^j/^ for the thermal distributions and 
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Fig. 35: The J/tp production per central LHC collision in the kinetic model as a function of the initial number of cc pairs. 



the assumption Ay = 4 are presented in Table ITSl 
5.4 Conclusions 

The "smoking gun" signature of the quarkonium formation mechanism is the quadratic dependence on 
the total number of charm quarks. For central collisions at the LHC one expects that this feature will 
lead to a total J /ijj rate greater than that produced by an incoherent superposition of the initial nucleon- 
nucleon collisions, even including suppression due to deconfinement effects. In addition, the centrality 
dependence can be used to identify the quadratic dependence on charm assuming that the initial charm 
production scales with the number of binary collisions. Binary scaling leads to an increase of the ratio 
of J /ip to initial charm as the collision centrality increases, independent of specific parameters which 
control the overall magnitudes. A simultaneous measurement of total charm will be essential for such 
conclusions to be drawn. 

Uncertainties in the absolute magnitude of the formation process are inherent in the model param- 
eters. For statistical hadronization, one can constrain the thermal parameters to within a factor of two 
using the observed hadron populations. There is some additional uncertainty related to the lower cutoff 
on centrality needed to ensure the quarkonium ratios are consistent with an overall thermal picture. There 
is also the possibility that the correction for canonical ensemble effects will involve a thermal volume 
parameter not necessarily equal to the total system volume 11811 . In addition, the formation mecha- 
nism could be limited to those charm quarks whose phase space separation is within some maximum 
value, introducing another as yet unconstrained parameter fTTOl. With kinetic formation, a similar set 
of uncertainties exist. There are uncertainties in the space-time properties of the deconfinement region. 
In addition, possible variations of charmonium binding energies and reaction cross sections in a decon- 
fined region are at present not well understood. There are indications that the efficiency of the formation 
mechanism is considerably reduced when included in a partonic transport calculation II182I . 

The primary uncertainty in both models is still the initial number of charm quarks and their mo- 
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Fig. 36: The centrality dependence of J/i/i production in the kinetic model. 



mentum distributions. The tabulated J /if) results should be regarded in this light. Thus numbers may 
be only an order of magnitude estimate. However, the variation with centrality and total initial charm 
should provide experimental signatures which are largely independent of the overall magnitudes. 
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Schroedter for updates on the Kinetic Model calculations. 



61 



0.40 



0.30 



E 

w 

03 

O 



9- 



— — © N^^(b=0)=100,Ay = 4 
» N^3 (b=0) =1 50, Ay = 4 

Q N^- (b=0) =200, Ay = 4 

►N^-(b=0)=100,Ay = 7 

i N^-(b=0)=150,Ay = 7 

A N^- (b=0) =200, Ay = 7 

^ N^-(b=0) =100, Thermal 
X N^-(b=0) =150, Thermal 
+ N - (b=0) =200, Thermal 



XX 




400 



N 



part 



Fig. 37: The ratio of the number of produced J/'ip's in the kinetic model to the initial number of cc pairs as a function of A'part- 



Table 17: Comparison of J/tp production at the LHC by the statistical hadronization (left-hand side) and kinetic formation 
(right-hand side) models. 





dNj/^/dy (Statistical) 

dNc-ciO)/dy 


Nj/^p (Kinetic, LO Charm) 

Nc-ciO) 


b{fm) 


25 


18.75 


12.5 


200 


150 


100 





0.656 


0.370 


0.165 


4.0 


2.26 


1.03 


1 


0.637 


0.359 


0.160 


3.85 


2.19 


1.00 


2 


0.586 


0.330 


0.147 


3.51 


2.00 


0.91 


3 


0.515 


0.290 


0.130 


3.04 


1.73 


0.79 


4 


0.434 


0.245 


0.109 


2.50 


1.43 


0.65 


5 


0.351 


0.198 


0.088 


1.97 


1.12 


0.51 


6 


0.270 


0.152 


0.068 


1.46 


0.84 


0.38 


7 


0.196 


0.110 


0.050 


1.01 


0.58 


0.27 


8 


0.132 


0.075 


0.034 


0.65 


0.38 


0.18 


9 


0.082 


0.046 


0.021 


0.38 


0.22 


0.10 


10 


0.045 


0.026 


0.012 


0.20 


0.12 


0.057 


11 


0.022 


0.013 


0.0061 


0.087 


0.054 


0.028 


12 


0.0097 


0.0058 


0.0029 


0.034 


0.022 


0.012 


13 


0.0045 


0.0029 


0.0016 


0.011 


0.0075 


0.0041 


14 


0.0028 


0.0019 


0.0012 


0.0021 


0.0013 


6.8 X 10-"^ 


15 


0.0025 


0.0018 


0.0012 


1.8 X 10""^ 


1.0 X IQ-'' 


5.1 X IQ-^ 
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Table 18: Kinetic J/ip formation at the LHC assuming both thermal charm momentum (left-hand side) and Ay = 4 (right-hand 
side). 





Nj/^ (Thermal) 
A^cc(O) 


Nj/^ (Ay = 4) 
A^cc(O) 


6(fm) 


200 


150 


100 


200 


150 


100 





52.7 


32.5 


16.4 


17.5 


10.8 


5.48 


1 


50.5 


31.2 


15.8 


16.8 


10.4 


5.25 


2 


44.8 


27.7 


14.0 


14.9 


9.21 


4.65 


3 


37.0 


22.9 


11.5 


12.3 


7.62 


3.82 


4 


28.6 


17.7 


8.73 


9.54 


5.89 


2.90 


5 


20.7 


12.7 


6.05 


6.90 


4.23 


2.01 


6 


13.8 


8.32 


3.72 


4.61 


2.77 


1.24 


7 


8.28 


4.71 


2.14 


2.76 


1.57 


0.71 


8 


4.10 


2.36 


1.10 


1.36 


0.79 


0.37 


9 


1.78 


1.04 


0.50 


0.59 


0.35 


0.16 


10 


0.65 


0.39 


0.19 


0.22 


0.13 


0.064 


11 


0.19 


0.12 


0.063 


0.065 


0.040 


0.021 


12 


0.048 


0.032 


0.018 


0.016 


0.010 


0.006 


13 


0.011 


0.0078 


0.0049 


0.0037 


0.0026 


0.0016 


14 


0.0026 


0.0019 


0.0012 


8.6 X 10""^ 


6.3 X 10-4 


4.1 X 10-4 


15 


5.9 X 10-4 


4.4 X 10-"^ 


2.9 X 10-* 


2.0 X 10"'' 


1.5 X 10-4 


9.7 X 10-^ 
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6. QUARKONIA IN HOT QCD MATTER: DISSOCIATION RATES 
6.1 Introduction 

Quarkonium suppression was long ago suggested as a signal of deconfinement UJ. Due to their small 
size, quarkonia can, in principle, survive the deconfinement phase transition. However, because of color 
screening, no bound state can exist at temperatures T > To when the screening radius, l//i/)(T), 
becomes smaller than the typical bound-state size |1|. Later it was realized that dramatic changes in 
the gluon momentum distributions at the deconfinement phase transition result in a sharp increase in the 
quarkonium dissociation rates 11751 11831 1184II . Both the magnitude I185II and the energy dependence 
[186 1 of charmonium dissociation by gluons result in a significant suppression of the cc states even 
for T < Tf) but higher than deconfinement transition temperature, Tc. Moreover, close to Td the 
thermal activation mechanism is expected to dominate 1 1 87. 188J . The relative importance of gluon 
dissociation and thermal activation is governed by the ratio of the quarkonium binding energy e{T) and 
the temperature T, X{T) = e{T)/T [189|. At X{T) < 1 thermal activation dominates while for 
X{T) » 1 the dominant mechanism is "ionization" by gluons. 

The study of quarkonium dissociation via color screening provides an upper bound on the tem- 
perature at which quarkonium bound states can exist. This temperature is also referred to as the Mott 
temperature 11901 . We note, however, that resonances may still exist above the Mott temperature 11911 . 

The dissociation temperature due to color screening was studied using potential models with dif- 
ferent parameterizations of the heavy quark potential Il27lll90lll92lll93l . All these studies predicted that 
excited charmonium states (xc, will essentially dissolve at Tc while the ground state J /ijj will dis- 
sociate at 1.1 Tc-1.3 Tc- Some potential models also predicted strong change in the binding energy, see 
e.g. Ref. II192I . Recently, charmonium properties were investigated using lattice calculations II194II1951 
which indicate that the ground states exist with essentially unchanged properties at temperatures around 
l.bTc- This issue will be discussed in section WJ\ Lattice investigations may indicate that at low tem- 
peratures, T < 1.5Tc, screening is not efficient and therefore gluon dissociation may be the appropriate 
source of quarkonium suppression, discussed in section 1631 Finally the equilibrium survival probabilities 
are estimated in section l6!4l 

One should keep in mind that nonequilibrium effects in the very early stages of a heavy-ion col- 
lision, when the energy density is highest, should be also considered for quarkonium suppression. Not 
much is known about these effects. However, they may be an even more important source of quarkonium 
suppression than those of the thermalized system, see e.g. Ref. 11961 . 



6.2 Investigation of In-medium Quarkonium Properties on tlie Lattice 

Temporal meson correlators A(r,p) may allow direct study of in-medium quarkonium properties since 
they are related to the T / spectral functions a{u!,p, T) by 

A/ ^ /""^ , / ^.T.NCOsh(u;(r- 1/(2T))) 

A{t,p)= (L;a{uj,p,T) ! . (63) 

Jo smh{uj / {2T)j 

Because at finite temperature, T 7^ 0, the extension of the time direction is limited by the inverse temper- 
ature, the only way to study meson properties is to reconstruct the spectral function'^. The reconstruction 
of the spectral function with any reasonable discretization of the lu -interval involves reconstruction of 
several hundred degrees of freedom from N^- ~ 10 data points on A(r, where Nr is the number 
of time-slices. In general, this can be done only using the so-called Maximum Entropy Method (MEM) 
L197..198 I. For a review see Ref. L198J and references therein. The method has also been quite successful 



at finite temperature 1 1991 12001 . 

'"'Authors: D. Blaschke, D. Kharzeev, P. Petreczky, H. Satz. 

'^At T — meson properties such as the mass, decay constants etc. can be extracted from the large distance behavior of 
meson correlators. 
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If we consider point sources, the correlators are defined as 

A{t,p) = {O{t,p)O{0,p)), (64) 

where 

0{t,p) = e'^^^q{T, x)Tq{T, x), (65) 

X 

and r = 1, 75, 7^ and 7^75 for the scalar, pseudo-scalar, vector and axial-vector channels. By checking 
the quantum numbers one realizes that the ^^o, '^Si, ^Pq and ^Pi states actually correspond to the 
pseudo-scalar, vector, scalar and axial-vector channels. 

One can also consider correlators of extended operators (sources) 

0^^*(r,p-) =X]e'^^X]c/>(y)g(r,x + y)rg(r,x), (66) 

X y 

where (l){y) is some trial quarkonium wave function which typically has the form exp(— ay^) II194II2011 . 
Extended operators have a better overlap with the ground state than a point source. Though the spec- 
tral functions of extended operators are easier to reconstruct due to their simplified structure, typically a 
single pole, they have the disadvantage of introducing a bias by the choice of some specific trial wave- 
function. They also do not correspond to any physically observable quantity such as the dilepton rate. 
The only information the study of such correlators can provide is the position and width of the reso- 
nance, provided that the corresponding structure of the spectral function can be interpreted as a narrow 
resonance. 

The following discussion is primarily based on the finite temperature charmonium analysis of 
Ref. P 195I . In this study, the nonperturbative clover action for quarks [202] was used with isotropic 
lattices, those with the same spacing in the temporal and spatial directions, with inverse lattice spac- 
ing, a~^, between 4.04 GeV and 9.72 GeV. The quark mass was varied in a range corresponding to 
the pseudo-scalar ijci^So) mass between 2.4 GeV and 4.1 GeV. The parameters of the simulations are 

Table 19: The lattice parameters used in the studies where is the lattice coupling, Tc is the critical temperature for decon- 
finement. The lattice spacing a is obtained from the string tension. The hopping parameter k is related to the quark mass by 
TUqa — {1/2k — 1/2kc)/uo where uo is the average value of the link field and K,c{g^) was calculated in Ref. \2Q2i . 





a-^ (GeV) 


Size 


T/T, 




6.499 


4.042 


48^ X 16 


0.93 


0.1300, 0.1234 


6.499 


4.042 


48^ X 12 


1.25 


0.1300, 0.1234 


6.499 


4.042 


48^ X 10 


1.5 


0.1300, 0.1234 


6.640 


4.860 


48^ X 16 


1.1 


0.1290 


6.640 


4.860 


48^ X 12 


1.5 


0.1290 


7.192 


9.720 


64^ X 24 


1.5 


0.13114 


7.192 


9.720 


48^ X 12 


3.0 


0.13114 



summarized in Table [T9l The results of this study are compared to another study fl94'| using a standard 
Wilson action for quarks with anisotropic lattices and extended operators. The temporal lattice spacing 
in this study was four times smaller than the spatial one, 1 / (2.03 GeV), allowing more points in the time 
direction while keeping the spatial volume reasonably large. On the other hand, the large spacing in the 
spatial direction may result in many unwanted lattice artifacts. Both studies were done in the quenched 
approximation, neglecting the effects of quark loops. This approximation was used because such stud- 
ies need quite small lattice spacings both because of the heavy quark mass and the need to have many 
points in the time direction. A',-, while keeping the temperature l/{aNr) = T fixed. Going beyond this 
approximation would require more than 10 Teraflops of computational power. 
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Fig. 38: The temporal correlators for different channels at = 6.499 and k = 0.13. The left-hand side shows the pseudo- 
scalar and vector channels, the right-hand side, the scalar and axial-vector channels. 



The temperature dependence of the temporal correlators alone can provide some information on 
the change of the meson spectral properties. In Fig.|38]we show the ratio of the correlators at different 
temperatures to the free correlators calculated on the lattice for quark masses corresponding to a ^S'o 
mass of 2.4 GeV. As seen in the figure, the correlators in the pseudo-scalar and vector channels show 
only a very small temperature dependence across Tc while in the scalar and axial-vector channels there 
are quite noticeable changes. This difference may imply that the ground state charmonia are not affected 
by deconfinement while the excited IP states are. This conclusion seems to be supported by the analysis 
of the spectral functions shown in Fig.|39| While the spectral function in the pseudo-scalar channel peaks 
at almost the same position at O.QSTc and 1.25Tc, in the axial- vector channel the peak present below is 
completely absent above Tc. The same results apply to the vector and scalar channels. Thus it seems that 
IS" ground state quarkonia are present above Tc with essentially unchanged masses but the excited IP 
states disappear already at T ~ 1.25Tc. The behavior of the correlators in Fig. |38] suggests that situation 
is similar at 1.5Tc. To see what actually happens at temperatures of l.BTc or higher one has to consider 
simulations with lattice spacings a" = 4.86 GeV and a'^ = 9.72 GeV with > 12. The analysis 
of the spectral functions in the pseudo-scalar and vector channels calculated with = 4.86 GeV 
shows that the ground state peak survives at least until 1.5Tc and the spectral functions are essentially 
unchanged for l.lTc < T < 1.5Tc while the scalar and axial- vector spectral functions are similar to 
those at l.STc. This conclusion is also supported by the analysis of the spectral functions at l.BTc with 
= 9.72 GeV. However, at 3Tc the spectral functions do not show any peak-like structures. On the 
other hand, the correlators show quite sizable changes between 1.5Tc and 3Tc. 

Finally we discuss the results obtained using extended operators 11941 . Here only the pseudo- 
scalar and vector channels were considered. The simulations were done at two temperatures, T = 0.88Tc 
and l.OSTc. The spectral functions were first extracted using the MEM. Due to use of extended sources, 
the spectral functions essentially consist of a single peak with a temperature-independent position. Thus 
the authors of Ref. LI 94 J used a fit to a Breit-Wigner function to determine the thermal width r(T), 
obtaining 

: r(r = l.OSTc) = 0.12 ±0.03 GeV, 
^S"! : r(T = l.OSTc) = 0.21 ±0.03 GeV. 

The picture which emerges from the results of the lattice studies described above can be summa- 
rized as follows. The ground state quarkonia survive up to ~ 1.5Tc. Moreover, its mass does not change 
significantly relative to its T = mass, contrary to potential model predictions which show significant 
mass shifts as well as dissociation of ground state charmonia well below 1.5Tc l27lll92l . Thus screening 
is not an efficient mechanism for charmonium dissociation at temperatures close to Tc. However, there 
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Fig. 39: The spectral function for the pseudo-scalar (right) and axial-vector (left) channels at O.QSTc and 1.25Tc and for 
6/5^ = 6.499 and K = 0.1300. 



are indications that the ground state charmonia acquire a thermal width, implying that collisions with 
gluons may be the dominant mechanism for quarkonium dissociation at low temperatures. As the tem- 
perature increases, screening becomes more and more important and will lead to the disappearance of 
any bound state. The present study only indicates that in the absence of dynamical fermions this disso- 
ciation will happen somewhere between 1.5Tc and STc. The excited IP states seems to dissolve already 
at ~ 1.25Tc, probably because the IP state is larger than the ground state and therefore screening may 
be sufficient to dissolve it. 

We now discuss whether these results can be interpreted within the present understanding of fi- 
nite temperature QCD. First we note that the detailed picture of screening depends on the quantity one 
considers. Universal screening can be understood in terms of a simple (perturbative) picture only at 
astronomically high temperatures due to additional infrared problems arising at finite temperature. For 
a review, see e.g. Ref. |203|. The most relevant quantity as far as the finite temperature properties of 
quarkonium are concerned is probably the color singlet free energy'^, Fi{r,T), studied in Ref. |205|. 
This quantity shows quite complex short distance behavior and cannot be parameterized by any simple 
form. Nevertheless, for Tc < T < l.STc, a rough estimate of the screening length is 0.4 — 0.6 fm. 
From the small value of the screening length, it is possible to conclude that the charmonium ground 
state should also dissolve below l.5Tc. However, the free energy also contains an entropy contribution, 
S{r, T) = —dFi (r, T) / dT and therefore is not directly related to the potential between static charges. In 
fact, at large distances the entropy is positive. Therefore the free energy, Fi (r, T) = Vi (r, T) —TS{r, T), 
is always smaller than the potential energy Vi{r, T) |205|. If this situation persists over all distances it 
would imply that screening in Fi is stronger than in Vi, and could possibly explain the discrepancy 
between the lattice studies of the charmonium spectral functions and potential models. 

6.3 Quarkonium Dissociation Cross Sections and Rate Coefficients 

In this section, we evaluate the dissociation rates of Coulombic bound states by collisions in a quark- 
gluon plasma. For the first step, the plasma is represented by an ideal gluon gas with thermal distribution 
functions ng{uj) = gg[ex.p{iOg{p) /T) — 1]^^, where T is the temperature, ujg{p) is the dispersion relation 
and Qg = 2{N^ — 1) is the gluon degeneracy factor with Nc = 3. We focus here on the quarko- 
nium ground states, J/V' and T, which we consider to be Coulombic bound states. Their masses are 
expected to be rather independent of the plasma temperature, as discussed in section 16.21 These results 
are consistent with experimental results for excitonic states in an electron-hole plasma [191 1 and with 

'*The gauge invariant definition of the color singlet free energy was suggested only very recently 12041 . So far, most lattice 
studies only considered the color-averaged free energy which also includes the octet contribution. 
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Fig. 40: Thermal averaged cross sections for heavy quarkonium dissociation by gluon impact as a function of temperature. 

dynamical screening fl901 as opposed to static Debye screening (193). The effects of lowering the con- 
tinuum threshold as well as temperature dependent modifications of the heavy quark potential will be 
neglected in this first step. 

The quarkonium breakup cross section by gluon impact can be estimated using the Bhanot-Peskin 
formula 11861 . 

2^^ (uj/e(0) — 1)^/^ 

^m)9^^^ = (j/,(o))5 - ^(0)) ' (67) 

where e(0) is the T = binding energy of the 15 quarkonium state with a Coulombic rms radius 
^ (r^)is = \/3ao = 2^/?,/ {asmq), mq is the heavy-quark mass and cj is the energy of the incoming 
gluon. The values of mq and e(0) used in our further calculations are given in Table 1^ 

Table 20: Values of binding energy, e(0), and heavy quark mass, mg, 12061 used in our further calculations. 





Set 


(i) 


Set 


(ii) 


System 


e(0) (GeV) 


mq (GeV) 


e(0) (GeV) 


mq (GeV) 


bottomonium 


0.75 


5.10 


1.10 


5.28 


charmonium 


0.78 


1.94 


0.62 


1.86 



The Bhanot-Peskin formula is derived using the operator product expansion, keeping the leading 
color-electric dipole operators to leading order in \/Nc, so that the final-state colour octet interactions 
of heavy-quark pair are negligible. The corresponding expression for the Wilson coefficients (chromo- 
electric polarizabilities) of the excited IP states is given in Ref. 12071 . Since quarks have been found to 
be less efficient in dissociating quarkonium states 11831 . we only include the effects of gluons here. 

The cross section in Eq. d67l has a steep rise at threshold with a drop after the maximum. We give 
the thermal-averaged cross sections in Fig. |^ A straightforward extension to the pT dependence of the 
quarkonium state in the plasma has been studied in Ref. 1. 184 J . 

The dissociation rate for a quarkonium IS state at rest in a heat bath of massless gluons at tem- 
perature T distributed according to ng{uj) is 

\qQ)^^^ " = • (68) 

Results for the partial widths r^ggsdue to collisions, together with their inverses corresponding to 
partial lifetimes, are shown in Fig. ^2 for gluon breakup. It is interesting to note that there is close 
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Fig. 41: The rate coefficients T^gg^^ — 't^qq^ for heavy quarkonium dissociation by gluon impact as a function of the gluon 
plasma temperature (left panel) and the corresponding lifetime ''"(gQ)g (right panel). For comparison, the Breit-Wigner fits to 
the ^So and '^Si charmonium spectral widths from lattice simulations are shown. 



correspondence between the charmonium decay widths due to gluon impact, the gluonic El transition, 
and the widths extracted from recent lattice simulations in the previous section. 

6.4 Quarkonia Abundancies and Observable Signatures 

In order to study observable signatures, we adopt the Bjorken scenario P2081 of plasma evolution, i.e. 
a longitudinal expansion with conserved entropy: T^t = TqTq = const. The initial temperatures and 
proper times are give in Table 1211 The RHIC and LHC values are set 3 from Ref. I184L assuming 
the highest initial fugacities. The SPS values are typical of those used in Ref. [2081. We calculate the 



Table 21: Initial conditions for plasma evolution. The results for LHC and RHIC are from Ref. flSl . 





LHC 


RHIC 


SPS 


To [GeV] 


0.72 


0.4 


0.25 


To [fm/c] 


0.5 


0.7 


1.0 



quarkonium survival probability neglecting nucleon absorption, hadronic comovers and hadronization 

effects, 

5(r,) = exp(-/Jdrr--^^(r)) . (69) 

At the freeze-out time, tj, the collisions stop changing the number of J/ ^p and T. Using Bjorken scaling, 
Tf can be translated into a freeze-out temperature. In Fig. |^we show the survival probability of J/ip 
(left panel) and T (right panel) due to gluon impact using the initial conditions from Table 1211 Note that 
set (i) gives the lowest T survival probability due to its lower values of the binding energy and bottom 
quark mass while set (ii), with its lower charmonium binding energy and charm mass, has the lower 
survival probability for the J/V'- The simple estimates presented here for the LHC suggest that none 
of the initially produced J/i/j's will survive their passage through the plasma. An analysis of the pT 
dependence of the suppression pattern, as in Ref. 1 184|, should provide a more complete picture. The T, 
however, will be a good probe of the plasma lifetime as well as its temperature. 
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Fig. 42: The survival probabilities of heavy quarkonia in a longitudinally expanding gluon plasma as a function of the plasma 
lifetime using the initial conditions of Table l21l and the heavy quark masses and binding energies of Table|20| On the left-hand 
side are the J/?/) results and on the right-hand side are the T results. 



6.5 Conclusions 

We have used recent lattice data on the mesonic spectral functions to constrain predictions of the temper- 
ature dependence of the two-particle quarkonium spectra. We have evaluated the collisional broadening 
of quarkonia states in a quark-gluon plasma using the Bhanot-Peskin cross section. The J/ip partial 
widths in a gluon gas are in good agreement with lattice results for the spectral widths without dynami- 
cal fermions. 

Neglecting secondary quarkonium production we have evaluated the time dependence of the sur- 
vival probability of J /ij; and T for different initial conditions. The results suggest that measurements of 
T suppression at the LHC can provide rather robust information about the initial temperature or lifetime 
oftheQGP|209|. 

Acknowledgments D.B. would like to thank Yuri Kalinovsky and Gerhard Burau for useful discussions. 
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7. HEAVY QUARK PRODUCTION IN MONTE CARLO GENERATORS ^ ' 

7.1 Introduction 

This section presents a study of the predictions of Monte Carlo parton shower programs for the genera- 
tion of charm and bottom quarks in heavy-ion collisions at LHC energies. In section IT2l the production 
cross sections in nucleon-nucleon collisions, calculated at the next-to-leading order (NLO), and their ex- 
trapolations to Pb-i-Pb collisions are reported. The NLO kinematical distributions are then compared with 
the results obtained with PYTHIA |210 211| and HERWIG L212.. 213. .214.1 in sect. 173] and sect. |T1 re- 
spectively. S ection FOl presents a short description of HIJING I215II216L a generator for nucleus-nucleus 
collisions, used to simulate the charm- and bottom-hadron production background in AA collisions. 

7.2 NLO heavy flavour production cross sections at the LHC 

The cross sections for charm and bottom production in pp collisions at LHC energies, ^/s = 5.5 TeV 
and 14 TeV, have been calculated at NLO using a program by Mangano, Nason and Ridolfi |5 1 (referred 
to as MNR hereafter). Two sets of parton distribution functions, MRST ||24| and CTEQ5M1 have 
been used. All the following MNR calculations have been performed using the same values for the heavy 
quark masses (rric and mi,) and the factorization and renormalization scales p and fiji) as in Ref . 1 177ll : 

rUc = 1.2 GeV , fip = fJ-R = 2 no , (70) 

rub = 4.75 GeV , = = , (71) 

where no = mq in calculations of the total cross section and = \J{Pt q + ^r.g^ + ^'^^ 
differential cross sections. The pp results are reported in Table |22l The difference due to the choice 
of the parton distribution functions is relatively small, ~ 20 — 25% at 5.5 TeV and slightly less at 14 TeV. 
In AA collisions, the cross section per nucleon is modified by shadowing which we implement using the 
EKS98 parameterization 1 16 17| . 

A discussion of the inclusive heavy quark pT and y spectra and the QQ correlations is presented 
in the following section. The NLO predictions are also compared to results obtained using PYTHIA. 

Table 22: NLO calculation fsl of the total cc and bb cross sections in pp collisions at 5.5 and 14 TeV, using the MRST and 
CTEQ5M1 parton distribution functions. 







app (mb) 


(TeV) 


QQ 


MRST 


CTEQ5M1 


5.5 


cc 


5.86 


7.42 


bb 


0.19 


0.22 


14 


cc 


10.28 


12.07 


bb 


0.46 


0.55 



7.3 Heavy flavour production in pythia 

The MNR code is not well suited for simulation purposes since it is not an event generator. Also, at 
LHC energies the NLO corrections are large. Therefore still higher order corrections may be important 
because of, e.g., multiple gluon emission. On the other hand, Monte Carlo event generators like PYTHIA 
and HERWIG offer a different set of approximations. They are based on LO pQCD matrix elements and 
parton showers which simulate real corrections to LO QCD processes. For example, a final state heavy 

'^Authors: N. Carrer, A. Dainese, H. Niemi. 
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quark can further radiate a gluon or an initial-state gluon can split into a quark-antiquark pair. The parton 
shower approach to contributions beyond LO is not exact even at NLO but it catches the leading-log 
aspects of multiple-parton emission. 

In the simulations described here, PYTHIA 6.150 was used. We will indicate light quarks (n, d or 
s) by q and charm or bottom quarks by Q. 

7.31 Comparison with LO MNR results 

On general grounds, one expects that the LO description of QQ production should be the same in PYTHIA 
and MNR when parton showers are turned off in PYTHIA. Thus before comparing the results of the full 
PYTHIA simulation with the ones obtained with the NLO MNR calculation, we first check if the basic LO 
subprocesses are handled correctly in PYTHIA. We compare the PYTHIA calculations without initial and 
final-state showers to the MNR LO results in ^/s = 5.5 TeV pp collisions for the following observables: 

• inclusive pT and rapidity distribution of the Q (Q); 

• pair pt, defined as the projection on the plane normal to the beam axis of the QQ total momentum 

PQ + pq; 

• pair mass, M{QQ) = ^{Eq + - (pg + Pg)^, where Eq = + (pg)^ is the quark 
energy; 

• angle Acf) between the Q and the Q in the plane normal to the beam axis. 

The results for charm production via gluon fusion'^, gg — > cc, are shown in Fig. 03] For a meaningful 
comparison the same input parameters were used in PYTHIA and in MNR. We use the CTEQ4L parton 
densities with a two-loop evaluation of Os, trie = 1-2 GeV, fip = fiR = fi = 2^ {jp]^'^^Y + 

where p^'^^ is the quark transverse momentum in the lab frame. Parton intrinsic momentum smearing 
is modeled by a Gaussian with {k?p) = 1 GeV^. Multiple interactions and initial and final-state parton 
showers were switched off in PYTHIA for this comparison. 

The level of agreement between the two calculations is excellent. Thus the LO calculations in both 
methods gives the same results when the same parameters are used. Since at LO the pair is produced 
back-to-back, the px of the pair should be zero and Ac^ = tt. The deviation from the strict back-to-back 
configuration is due to the intrinsic kx of the incoming partons. The /c^—kick is implemented differently 
in the two calculations. Thus the fact that MNR and PYTHIA give the same results for this value of {kj,) 
is quite remarkable. However, this fact should be considered as accidental, and one should not expect it 
to hold in general, for example for larger values of {kj.). 

7.32 Parton shower and process classification in PYTHIA 

When the parton shower is included, the partonic interaction can be subdivided in different stages. The 
event generation starts with the LO process and then the parton shower is used to generate different 
histories for the incoming and outgoing particles that take part in the LO process. For example, outgoing 
particles from the LO process can emit gluons or an incoming charm quark could come from the splitting 
of a gluon into a QQ pair. Therefore the parton shower will shift the kinematics of the produced particles 
and at the same time can increase the multiplicity of final state. The parton showers do not affect the 
probabilities of the LO processes but increase the total rate of heavy quark production through gluon 
splitting into QQ pairs. 

In PYTHIA the processes are divided in the following three classes, according to the number of 
heavy quarks in the final state of the hard process, where the hard process is defined as the one with the 
highest virtuality: 

'*In PYTHIA, gg cc is called by ISUB=82. 
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pair creation: The hard process produces two heavy quarks in the final state tlirough the reactions gg — > 
QQ and qq QQ; these are the reactions which constitute the LO contribution in MNR; 

flavour excitation: A heavy quark in one incoming hadron is put on mass shell by scattering off a 
parton from the opposite hadron, qQ — > qQ or gQ —>■ gQ. The incoming heavy quark is assumed 
to come from a.g QQ splitting in the initial-state shower. This hard process is characterized by 
one heavy quark in the final state; 

gluon splitting: No heavy quark is involved in the hard scattering. Instead a QQ pair is produced in the 
initial or final-state showers from a. g ^ QQ branching. 

This classification might be misleading since all three classes produce pairs at g ^ QQ vertices when 
showers are turned on. However, the requirement that the hard scattering should be more virtual than 
the shower avoids double-counting L217.I . Fig. l44l shows some topologies belonging to the processes 
specified above. 

Flavour excitation and gluon splitting contributions are calculated in the approximation of mass- 
less quarks. As a consequence, these cross sections diverge as p^?'^'^ vanishes. These divergences are 
usually regularized by putting a lower bound on the allowed values of pq?^'^. However, the choice of the 
minimum pJp^'^ is arbitrai^y and is usually considered to be a parameter. The lower bound on p^^'^'^ has 
a strong influence on the heavy flavour cross section at low px- Experiments with acceptances starting 
at several GeV, such as CDF and D0 at the Tevatron and ATLAS and CMS at the LHC, are not very 
sensitive to p^^'^'^. On the other hand, ALICE has been designed for good acceptance down to very low 




Fig. 43: Comparison between PYTHIA results (solid histogram) for gg — > cc, without parton showers, and the MNR calculation 
of the same process at LO (dashed histogram). The centre of mass energy is ^/s = 5.5 TeV. 
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Pt since the measurement of the total heavy flavour cross section in nuclear reactions is of great inter- 
est. Therefore, the choice of the minimum pj?'^'^ in the simulations of heavy quark production has to be 
carefully evaluated. The criteria used in this choice are explained in the next section. 

7.33 Comparison with MNR at NLO and tuning of PYTHIA parameters 

Because of the divergence at low px, PYTHIA does not offer a valid baseline for the calculation of the total 
cross section and of the kinematical distributions in the low-p^ region. For this reason our approach has 
been to adopt the NLO calculation as the baseline for our simulations and to tune the PYTHIA parameters 
to reproduce the NLO predictions as closely as possible. The main parameter tuned is the lower bound 
of 

phard jj^g simulations of Pb+Pb collisions at = 5.5 TeV/nucleon includes shadowing using the 
EKS parameterization I 16..17.I . as noted in section IT2l The complete list of parameters used is shown 
in Table l23l See Ref. 12101 12111 for a description of the PYTHIA parameters. The larger value of the 
intrinsic fc-r-kick in Pb+Pb relative to pp collisions (fcy broadening) has been taken from Ref. [177|, 
based on a Q^-dependent nuclear broadening. 

The results of this tuning for Pb+Pb collisions at 5.5 TeV/nucleon are shown in Fig.|45lfor charm 
and in Fig. |46l for bottom. The PYTHIA and NLO distributions are compared. The PYTHIA results are 
scaled by factors of 0.62 and 0.57 respectively to give the same total cross section as the NLO calculation. 
Despite the fundamental differences between the two models, the agreement is relatively good. However, 
significant discrepancies are present, especially in the Ac/) distribution of the charm quarks. 

We note that the CTEQ4L parameterization used here is rather outdated. However, we verified 
with the MNR code, that these results lie in between those obtained with CTEQ5 and MRST for all the 
relevant kinematical quantities. The sensitivity to the parton distribution functions does not exceed 30% 
for charm and 20% for bottom. Figure |^ shows the inclusive px and rapidity distributions for charm 
with the CTEQ4L, CTEQ5L and MRST distributions 

To simulate heavy quark production for the 14 TeV pp run, the same parton distribution functions 
and values of mq and p are used, see Table |23] The comparison of the PYTHIA spectra with the MNR 
calculation is shown in Fig.08]for charm and in Fig.09lfor bottom. The largest difference is a rather poor 
description of the charm quark rapidity distribution due to the range of validity of the parton distribution 
functions. Most widely used parton distributions, including CTEQ4, are valid only down to x ~ 10^^. 
Below this x, the behaviour depends on the implementation of the parton density but has no physical 
meaning. For example, the CTEQ4 gluon density is fixed to its value at x = 10^^. The rapidity range 
over which the evolution of the set of parton densities is reliable depends on the energy. For charm 

"a direct comparison withi tire results in Fig. l45l is not possible since a slightly different set of parameters was used. 

Pair Creation Flavour Excitation Gluon Splitting 




Fig. 44: Examples of pair creation, flavour excitation and gluon splitting. The thick lines correspond to the hard process, the 
thin ones to the parton shower. 
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Fig. 45: Comparison between charm production in the NLO calculation by Mangano, Nason and Ridolfi and in PYTHIA with 
parameters tuned as described in the text for Pb+Pb collisions at y^s = 5.5 TeV/nucleon. The triangles show the NLO 
calculation, the solid histogram corresponds to the PYTHIA total production. The individual PYTHIA contributions are pair 
production (dashed), flavour excitation (dotted) and gluon splitting (dot-dashed). 
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Fig. 46: Same as Fig. l45l but for bottom production. 
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production at ^/s = 5.5 and 14 TeV this range corresponds to \y\ < 4.3 and \y\ < 3.4 respectively. 
7.34 Fragmentation and decay 

We use the default pythia quark fragmentation and hadron decay parameters. Tablel24lsummarizes the 
total yield and the rapidity density, dN/dy, averaged over \y\ < 1 for heavy flavour hadrons. 

The default fragmentation scheme in PYTHIA is Lund string fragmentation. The string topology is 
derived from the colour flow in the hard process, determining the partonic partners in the colour singlet 
state forming the string endpoints. Let us consider, for example, one of the LO graphs contributing to 
gg — > bb, see Fig. |50| The colour representation of the beam remnants is always described by a quark 
in a colour triplet state and a diquark in a colour anti-triplet state which thus behaves like an antiquark. 

Table 23: PYTHIA parameters used for the generation of charm and bottom quarks for Pb+Pb collisions at ^/s = 5.5 
TeV/nucleon and pp collisions at -^s = 14 TeV. All unspecified parameters are PYTHIA 6.150 defaults. 



description 


parameter 


Charm 


Bottom 


Process types 


MSEL 


1 


1 


Quark mass 


PMAS(4/5,1) 


1.2 


4.75 


Minimum pl^^"^ 


CKIN(3) 


2.1 


2.75 


CTEQ4L 


MSTP(51) 


4032 


4032 


Proton PDF 


MSTP(52) 


2 


2 


Switch off 


MSTP(81) 








multiple 


PARP(81) 








interactions 


PARP(82) 








Initial/Final parton 


MSTP(61) 


1 


1 


shower on 


MSTP(71) 


1 


1 


2"^^ order 


MSTP(2) 


2 


2 


QCD scales 


MSTP(32) 


2 


2 


for hard scattering 


PARP(34) 


1 


1 


and parton shower 


PARP(67) 


1 


1 




PARP(71) 


4 


1 




MSTP(91) 


1 


1 




PARP(91) 


1.304 (Pb+Pb) 


2.035 (Pb+Pb) 


Intrinsic kr 




1 (pp) 


1 (pp) 




PARP(93) 


6.52 (Pb+Pb) 


10.17 (Pb+Pb) 






5 (pp) 


5 (pp) 



Table 24: Total yield and average rapidity density for |y| < 1 for hadrons with charm and bottom produced in Pb+Pb collisions 
at ^/s = 5.5 TeV/nucleon. 



Particle 


Yield 


{dN/dy) 


Particle 


Yield 


(diV/dy)l^l<i 




68.9 


6.87 




1.86 


0.273 


d' 


71.9 


6.83 


b' 


1.79 


0.262 


D+ 


22.4 


2.12 


B+ 


1.82 


0.251 


D- 


22.2 


2.00 


B- 


1.83 


0.270 


Dt 


14.1 


1.30 




0.53 


0.077 


D- 


12.7 


1.22 


b: 


0.53 


0.082 




9.7 


1.18 


At 


0.36 


0.050 


Ac 


8.2 


0.85 


0.31 


0.047 
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Fig. 47: Inclusive pr and rapidity distributions for cliarm withi the CTEQ4L, CTEQ5L and MRST structure functions. Tlie 
distributions are obtained using the MNR code for Pb+Pb collisions at 5.5 TeV/nucleon. 

Therefore, in this simple situation, there aie only three strings, as shown in Fig. |50| The b quark is 
colour connected to a diquark, the b is connected to a quark and the last string connects the remaining 
quark-diquark system. When we include the parton shower, the number of quarks in the final state 
increases and therefore more colour singlet strings can be produced. The strings then fragment according 
to the longitudinal Lund symmetric fragmentation function. For strings with heavy quark endpoints, the 
Lund fragmentation function is too soft to reproduce the charm and bottom meson fragmentation data. 
The Lund fragmentation function is therefore modified through the PYTHIA default settings to take into 
account the harder spectrum of heavy flavour mesons using the Bowler parameterization. The fact that 
the heavy quark is usually connected with one of the beam remnants shifts the rapidity of the heavy quark 
towards the rapidity of the projectile or target 12171 . see Fig.lSTl 

A detailed description of the heavy quark string fragmentation can be found in Ref. [2171. Here 
we report only how it affects the kinematics. Figure|52lshows the charm pT and rapidity as well as the pT 
and A(/) distributions of the cc pairs. The solid histograms represent quarks and the dashed histograms, 
the charm mesons. In this figure we consider only the mesons into which the quarks hadronize. The 
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Fig. 48: Comparison between chiarm production in 14 TeV pp interactions with the MNR code and PYTHIA with parameters of 
Table l23l The curves are the same as those in Fig. l43l but now the MNR results are shown by the solid curves. 
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Fig. 49: Same as Fig. l48l but for bottom production. 
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decay of heavy flavor resonances is not included. Figure |53] shows the equivalent plots for bottom. 

The most important effect is a reduction of the transverse momenta, directly reflecting the mag- 
nitude of the longitudinal fragmentation function. The fragmentation reduces the momentum of charm 
mesons by 25% and of bottom mesons by 14% on average. The rapidity shift shown in Fig.|^is much 
smaller than the rapidity range of the interactions. Therefore it does not strongly affect the inclusive 
rapidity distribution. The azimuthal correlation is also not significantly affected by fragmentation. 

7.4 A study of heavy flavour production in herwig 

Like PYTHIA, HERWIG 12 121 12 131 12 141 is a general-purpose event generator which uses parton showers 
to simulate higher-order QCD effects. The main differences with respect to PYTHIA are the treatment of 
the parton showers, which accounts more correctly for coherence effects, and the hadronization model, 
based on cluster rather than string fragmentation. The version used for the analysis presented here is 
HERWIG 6.4. Here we investigate the main features of heavy flavour production in HERWIG rather than 
trying to produce agreement between HERWIG and the NLO calculation. The following simulations have 
been done for pp collisions at ^/s = 5.5 TeV. 

7.41 Comparison with LO MNR results 

As with PYTHIA, we first checked the description of the LO processes in HERWIG by comparing to the 
LO results of the MNR code. For the LO calculations, we used the same parameters in both calculations. 
We used the CTEQ4L parton distributions with A^^^ = 236 MeV, = 1.2 GeV, rri}, = 4.75 GeV and 
default HERWIG factorization and renormalization scales. Good agreement was found when the HERWIG 
parton showers were switched off. 

7.42 Parton showers in HERWIG.' gluon splitting and flavour excitation 

Since gluon splitting and flavour excitation cannot be considered separately in the NLO code, as in event 
generators, we cannot compare the two results at this level. Instead we can compare HERWIG with 
PYTHIA to see if there are fundamental differences in the implementation of these processes. 

The gluon splitting results for bottom production are compared in Fig. |54] The shape of the 
distributions from the two generators agree at high px. Since different scales are used in HERWIG and 
PYTHIA we do not expect the result to be exactly the same, especially at low px- 

On the other hand, flavour excitation is clearly different in the two generators. The problem is 
apparent in the single quark px spectrum where a large peak is present at px ~ when HERWIG is used. 




Fig. 50: Sketch of the string structure for the LO graph gg bb. Left: the heavy quark is attached to a beam remnant (quark 
or diquark). Right: the same view in momentum space (not to scale). The dashed lines represent strings. 
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Fig. 51: Average difference between tlie rapidity of the B meson and the one of the corresponding b quark, as a function of the 
quarlc rapidity, for interactions with p!^''"''^ > nib- 




Fig. 52: Distributions for charm quarks (solid histogram) and D mesons (dashed histogram). The normalization is arbitrary. 
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Fig. 53: Distributions for bottom quarks (solid histogram) and B mesons (dashed histogram). The normalization is arbitrary. 
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Fig. 54; Gluon spUtting results in bb production, obtained with herwig (solid) and pythia (dotted). 
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We studied the origin of this peak by comparing the single quark spectrum for the quark coming from 
the hard scattering and for the quark coming from the initial-state gluon splitting, the middle diagram in 
Fig. 1331 We have performed this separation and we present the comparison in Fig. |55] We can see that, 
while the hard scattering part is compatible for the two models, HERWIG fails to describe the initial-state 
evolution of the gluon splitting into a QQ pair. This may be due to the fact that initial-state shower 
variables in HERWIG ai^e chosen for light parton showers and mass effects are not considered |218 1. We 
found the same kind of behaviour for charm production although the peak in spectrum is shifted to 
a value closer to nic. Given that HERWIG does not correctly describe flavour excitation in the low-pr 
region, it should not be used for simulations in such a region. 

7.5 HIJING: event generator for nucleus-nucleus collisions 

HIJING (Heavy Ion Jet INteraction Generator) is widely used to simulate the background for feasibility 
studies of heavy flavour detection. HIJING combines a QCD-inspired model of jet production 112151 
|2161 with the Lund model for jet fragmentation. Hard or semi-hard parton scatterings with transverse 
momentum of a few GeV are expected to dominate high energy heavy ion collisions. The HIJING model 
has been developed with special emphasis on the role of minijets in pp, pA and AA interactions at 
collider energies. 

The Lund FRITIOF f2191 model and the Dual Parton Model II220I have guided the formulation 
of HIJING for soft nucleus-nucleus interactions at low energies, ^/s ~ 20 GeV/nucleon. The hadronic 
collision model has been inspired by the successful implementation of perturbative QCD processes in 
PYTHIA. 

Two important features of nuclear collisions in HIJING are the simulation of jet quenching and 
nuclear shadowing. The charged-particle rapidity-density given by HIJING with jet quenching for central 
Pb-i-Pb collisions at ^/s = 5.5 TeV/nucleon is dNch/dy ~ 6200, about a factor of 2 larger than that pre- 
dicted using recent analyses of RHIC results |221 1. However, such a large multiplicity can be considered 
as a "safety factor" for the feasibility studies. HIJING also provides the relevant background sources such 
as direct leptons and strange particles. 

7.6 Conclusions 

We have studied heavy quark production in nucleon-nucleon collisions at LHC energies with the most 
widely used Monte Carlo generators — pythia and herwig. The resulting distributions were compared 
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to those obtained in perturbative QCD to NLO using the MNR code [5| to estabUsh a baseline simu- 
lation of heavy flavour production. We found that by tuning the main PYTHIA parameters, satisfactory 
agreement with the NLO bare quark transverse momentum distributions was obtained for both charm and 
bottom production in pp collisions at ^/s = 5.5 and 14 TeV. The most important parameters tuned in the 
comparison were the heavy quark masses, the minimum p^^'^ of the heavy quarks in the lab frame, and 
the width of the transverse momentum smearing of the incoming partons, {kj.). However, we found that 
HERWIG does not treat flavour excitation correctly at small transverse momenta, making it unsuitable for 
simulations of heavy quark production. 

This problem is solved by merging the NLO computations with parton showers, as advocated in 
MC @NLO 1 10 222,1 . The implementation of heavy flavour hadroproduction, Ref . 1 10] , became available 
after the completion of the work reported in this section. We point out that, in the context of MC@NLO, 
no additional tuning of the Monte Carlo parameters is needed to reproduce the NLO results in the relevant 
regions since this is guaranteed by the formalism. 
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8. QUARKONIA AND HEAVY FLAVOUR DETECTION WITH THE ALICE DETECTOR ^" 
8.1 Introduction 

ALICE is the dedicated heavy-ion experiment at the LHC. The apparatus will investigate strongly inter- 
acting matter at extreme energy densities by comprehensive measurements of a large variety of observ- 
ables 12231 12241 . Among the most promising observables, heavy flavours are especially relevant since 
they provide an essential probe of the earliest stage of heavy ion collisions. They further give precious 
information on the properties of the medium produced at longer time scales. 

Quarkonium states can be pertinent signatures of the Quark Gluon Plasma (QGP). From the early 
predictions of charmonium suppression by Debye screening in a deconfined medium (2, to the recent re- 
sults from the NA50 collaboration at SPS [225 1, much effort has been devoted to the subject (for reviews 
see Refs. 12261 1227 1). While at SPS energies only charmonium states are experimentally accessible, the 
much higher LHC (and RHIC) energies make bottomonium measurements feasible, thus providing an 
additional probe for QGP studies. In fact, since the T(1S) state only dissolves significantly above the 
critical temperature l27llI88l . at a value which might only be reachable at energies above that of RHIC, 
the spectroscopy of the T family at the LHC should reveal unique characteristics of the QGP |29|. We 
note also that quarkonium measurements at RHIC and LHC are complementary in terms of the accessible 
X range. 

The study of quarkonium states at the LHC is significantly different from those at the SPS and 
RHIC. First of all, in addition to prompt charmonia produced by hard scattering, significant indirect 
charmonia can be produced by B meson decays |228|, DD annihilation |229|, and by coalescence 
mechanisms which could result in enhancement rather than suppression, see section |^ Then, in the 
environment of a heavy-ion reaction, in-medium effects such as shadowing and heavy quark energy 
loss may substantially modify the final yields and spectra, see section ^ It is therefore obvious that 
an understanding of the QGP requires systematic investigations. More precisely, quarkonium must be 
studied: 

• as a function of centrality to identify suppression/enhancement patterns; 

• as a function of system size to vary the energy densities and thus disentangle normal and anomalous 
suppression; 

• for all species because the survival probability reflects the temperature; 

• as a function of pT to disentangle models; 

• with good vertex resolution to distinguish between prompt and secondary production; 

• with respect to the orientation of the reaction plane to unravel Glauber and comover absorption; 

• together with other QGP signals to explore correlations; 

• together with open charm and bottom, the most natural normalization of the quarkonium signals. 
This physics program should be achievable with the ALICE detector. Indeed, quarkonium states will be 
identified both in the dielectron and in the dimuon channels. 

The study of open charm and bottom in heavy ion collisions at the LHC is another important issue 
addressed by ALICE since it probes the mechanisms of heavy quark production in the hot and dense 
medium formed in the early stage of the collision. In this sense, comparison with pp and pA interactions 
will be essential for establishing the baseline production cross sections and nuclear shadowing effects. 

If a QGP is formed, secondary parton scattering may provide an additional source of charm 
quarks P 2301 12311 . In such a scenario, the number of heavy quarks produced would strongly depend 
on the initial temperature and lifetime of the plasma, providing information on the early stage of the 
collision. 

Plasma formation would not only modify the total production cross section but would also affect 
the kinematic distributions of the produced heavy quarks. Elastic collisions of the heavy quarks with 

^"Authors: P. Crochet, A. Dainese, E. Vercellin. 
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partons in the plasma and radiative loss in the medium may reduce the heavy quark momenta, see e.g. 
Ref. I1137II . Energy loss by D and B mesons, originating from c and h quarks, is particularly relevant 
because heavy flavor energy loss is expected to be significantly less than for light hadron production, 
dominated by light parton fragmentation at LHC energies, see Section 14.31 and the jet chapter of this 
report for more discussion. 

The present experimental picture of heavy flavor production in heavy ion collisions is quite unclear. 
The dimuon spectrum measured by NA50 I232II at the SPS may indicate enhanced charm production 
in central Pb+Pb collisions. On the other hand, PHENIX measurements of D production through the 
electron yield at RHIC suggest no charm enhancement in Au+Au collisions Iil48i . 

ALICE is equipped with dedicated subdetectors for the identification of secondary vertices with a 
displacement of ~ 100 /im from the primary interaction vertex. These displaced vertices are the primary 
signature of heavy flavour decays. Inclusive D and B meson production will be measured through their 
semi-leptonic decays. Moreover, the exclusive reconstruction of hadronic D meson decays will provide 
a direct measurement of the D pT distribution. 

The expected detector performance in the detection of hidden and open heavy flavours is presented 

here. 

8.2 The ALICE Detector 




Fig. 56: Schematic view of the ALICE detector. 

The ALICE detector, shown in Fig.|56l is designed to cope with large particle multiplicities which, 
in central Pb+Pb collisions, are expected to be between 2000 and 8000 per unit rapidity at midrapidity. 
It consists of a central part, a forward muon spectrometer and forward/backward small acceptance de- 
tectors. The central part of ALICE consists of four layers of detectors placed in the small, < 0.5 T, 
solenoidal field provided by the LEP L3 magnet. From the inner side to the outer side, these detectors 
are the Inner Tracker System (ITS), the large Time Projection Chamber (TPC), the Transition Radiation 
Detector (TRD) and the Time of Flight system (TOE). They provide charged particle identification in the 
pseudorapidity range |r/| < 0.9, with full azimuthal coverage and a broad px acceptance. These lai^ge 
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area devices are complemented by smaller acceptance detectors: the High Momentum Particle Identi- 
fication (HMPID), the PHOton Spectrometer (PHOS) and the Photon Multiphcity Detector (PMD). In 
the forward/backward region, the charged multiplicity and the zero degree energy will be measured by 
additional detectors (TO, VO, FMD, ZDC) which will allow fast characterisation and selection of the 
events. Finally a forward muon spectrometer covering the pseudorapidity range 2.5 < 77 < 4 is placed 
on the right side of the central part. It makes use of the usual techniques of muon identification at small 
angle and consists of a front absorber, a dipole magnet, ten high-granularity tracking chambers, a muon 
filter and four large area trigger chambers. 



8.3 Accessible x Range 

The LHC can probe both the nucleon parton distribution functions, in pp collisions, and their modifi- 
cations in the nucleus, using heavy ion collisions, down to unprecedented low x values where x is the 
fraction of the nucleon momentum carried by the interacting parton. In this section, the x range of heavy 
flavour production accessible to ALICE is quantitatively assessed. 

The X range probed depends on the center of mass energy per nucleon pair ^/s, the QQ pair 
invariant mass, Mqq, and on the pair rapidity, Vq-q- We consider the leading order process, gg — > QQ 
in the collision of two ions with mass and charge {Ai, Z\) and (^2, Z2). The square of the pair invariant 
mass is equal to the square of the center of mass energy of the initial gluons, s. 



Zn 



X\ X2 S 



Ai A2 



(72) 



where xi and X2 are the gluon momentum fractions and Spp = (14 TeV)^ is the square of the pp center 
of mass energy at the LHC. The longitudinal rapidity of the pair in the laboratory is given by: 
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From the pair mass and rapidity, the dependence of xi and X2 on A,Z, Mqq and Dq-q can be derived, 
obtaining 

xi = — ^exp^yggj :,2 = —^e^v[-yQQ)- (74) 



'pp 



We first focus on Pb-i-Pb and pp reactions. At central rapidities, xi ~ X2 with the magnitude 
determined by the ratio MqqI The x values for QQ production at threshold, Mcc = 2mc = 2.4 GeV 
and M,T = 2m;, = 9 GeV, are given in Table l25l 



Table 25: The x values for QQ pair production at threshold for central rapidities in Pb+Pb and -pp collisions. 



system (y^) 


X (cc) 


X (bb) 


Pb-FPb (5.5 TeV) 
pp (14 TeV) 


4.3 X 10-4 
1.7 X 10-4 


1.6 X 10-^ 
6.4 X 10-4 



Because of the lower mass, charm probes smaller x values than bottom. The capability for charm 
and bottom measurements in the forward rapidity region, y ^ 4:, gives access to x regimes about 2 orders 
of magnitude lower, down to x ~ 10^^. 

In Fig. 1^ the regions of the (xi, X2) plane for charm and bottom measurements in the ALICE 
acceptance are shown for Pb-i-Pb collisions at 5.5 TeV and pp collisions at 14 TeV. In this plane, points 
of constant invariant mass he on hyperbolae, xi = Mqq/ {x2 s), which are straight lines on the log-log 

scale. The solid and dashed lines correspond to cc and bb pair production at threshold. Points with 
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Fig. 57: ALICE heavy flavour acceptance in the (xi, X2) plane for Pb+Pb (left) and pp (right) collisions. 




I I Central barrel: hi < 0.9 
Q Muon arm: 2.5 < r) < 4 




Fig. 58: ALICE heavy flavour acceptance in the {xi, X2) plane for pPb (left) and Pbp (right) collisions. 
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Fig. 59: Compiled ALICE acceptance in the {x\, xi) plane for charm (left) and bottom (right) in pp, Pb+Pb, pPb and Pbp 
colUsions. 
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constant rapidity lie on straight lines, xi = X2 exp(— 2?/qq). The shaded regions show the acceptance of 
the ALICE central barrel, \r]\ < 0.9, and of the muon arm, 2.5 < r/ < 4. 

In the case of asymmetric collisions, Zi / Z2 and Ai A2, the QQ rapidity is shifted by: 



as seen in Eq. dTSb . For pPb (Pbp) collisions. Ay = +0.47 (—0.47). Therefore, combining pPb and Pbp 
runs will allow the largest xi and X2 coverage with the central barrel and the muon arm. FigurelSSlshows 
the pPb and Pbp acceptances while Fig. |59l compares the pp, Pb+Pb, pPb and Pbp coverages for charm 
(left) and beauty (right). 

Note that these figures only give a qualitative idea of the x regions accessible to ALICE since 
the leading order definitions of x were used and the rapidity cuts were applied to the QQ pair and not 
the detected particles. In addition, no minimum px cuts were applied. Such cuts increase the minimum 
accessible value of Mqq, correspondingly increasing the minimum accessible x. These approximations, 
however, are not too drastic since there is a very strong rapidity correlation between the initial QQ pair 
and the produced heavy flavour hadrons and, as shown in the following sections, the minimum px cut 
will be quite low, lower than the hadron mass, for most of the channels studied at ALICE. 

8.4 Muons 

The ALICE forward muon spectrometer is designed to detect heavy quarkonia in the decay chan- 

nel. The continuum will be measured together with the resonances allowing the study of open 

charm and open bottom. 

8.41 Quarkonia 

Muons are detected in the pseudo rapidity interval 2.5 < r/ < 4. To eliminate the huge low-pT back- 
ground from vr and K decays, a muon pT threshold is applied to single muons at the trigger level. The 
resulting geometrical acceptance is shown in Fig.|60lfor J /ijj and T. 

The projection of the 3-D plots on the y and px axes are shown in the middle and right panels of 
the same figure. Together with these projections, the ones obtained after applying a sharp px cut of 1 
(2) GeV/c for J/ijj (T) are also shown to give an idea of the trigger effects. As can be seen, only the J /ijj 
acceptance is slightly affected by the trigger px cut, implying that J /ijj and T can be detected down to 
px = 0. 

An important characteristic of the spectrometer is that its mass resolution is expected to be of 
about 70 (100 MeV/c^) for J/iJj (T). As shown in Fig. |^ such resolution allows separation of the T 
states. It should then be possible to disentangle scenarios of quarkonium suppression by studying the px 
dependence of quarkonium ratios [29|. 

The statistics expected in a 10^ s run, roughly corresponding to one month of data taking, are 
summarized in Table l26l for different colliding systems. (Note that since Ca has been replaced by Ar 
in the hst of accepted nuclei, the Ca-i-Ca results are indicative of the expected rates for an intermediate 
mass system.) The signal-to-background ratios and the significances for detecting the various resonances 
are also given. We note that these numbers were obtained by applying analysis cuts which are not fully 
optimized. Therefore they should be considered as indicative numbers only. All numbers refer to primary 
production of the resonances and do not include any specific in-medium effects such as enhancement or 
suppression mechanisms. The total quarkonium production cross sections in proton-proton collisions 
were taken from Ref. [20|. These cross sections have been scaled to central (10%) nucleus-nucleus 
reactions according to /{2 ■ 2) where the factors of 2 represent the reduction for central relative to 
minimum bias collisions and nuclear shadowing 12331 . The px distributions have been parameterized 
according to CDF results (see Ref. 112331 for more details). 




(75) 
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Fig. 60: Geometrical acceptance for J/?/; (top) and T (bottom). Tiie left panels show the geometrical acceptance versus 
rapidity and transverse momentum. The middle and right panels show the acceptance as a function of rapidity and transverse 
momentum, respectively. In the middle and right panels, the acceptance is shown without (solid histograms) and with (dashed 
histograms) the trigger cut on single muon pr- From Ref. 12241 . 

As shown in the table, the high statistics expected for J/ tp will permit detailed studies as a function 
of both pt and centrality. It is then also possible to investigate the azimuthal distribution of J/ tp with 
respect to the reaction plane |235 1. For the T, the statistics are obviously reduced but are nevertheless 
sufficient to investigate the px and centrality dependence. Note particularly the large expected signal-to- 
background ratio (about 10) and good significance. 

8.42 Open Charm and Open Bottom 

The muon spectrometer should allow measurements of the open bottom cross section. This is possible 
because the background muons from vr and K decay surviving the 1 GeV/c trigger px cut, as well as 
those coming from charm decay, can be efficiently removed with a higher software pr cut. In fact, 
as illustrated in Fig. |62l when applying a pT cut of 3 GeV /c on each muon, the correlated signal from 
bottom decays represent a significant fraction of the total dimuon yield both at low (1 < m < 3 GeV) and 
high (6 < m < 9 GeV) invariant mass. Therefore, after subtraction of the uncorrelated combinatorial 
background (evaluated, for instance, by event mixing) from the total dimuon yield, one should be able to 
extract the unlike-sign correlated bottom decay signal in these two invariant mass regions. 

It is important to note that the unlike-sign correlated signal from bottom decay has two different 
origins. In the high invariant mass region (Fig.|62j right) each muon comes from the direct decay of a i? 
meson. In the low invariant mass region (Fig. |62j left), the correlated signal from bottom is dominated 
by the so-called S-chain channel where both muons come from the decay of a single B via a D. This 
implies that the bottom cross section can be estimated independently (at least to some extent) by a study 
of the low and high invariant mass regions. Note that similar measurements of the charm cross section 
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Fig. 61: Opposite sign dimuon mass spectra in the region 2 < m < 4 GeV (left) and 8 < m < 12 GeV (right) for the 
10% most central Pb+Pb reactions. The spectra were obtained by means of fast simulations including acceptance cuts and 
detector efficiencies and resolutions. They correspond to a luminosity of 5 x 10^" cm^^s^^ and a running time of 10^ s. From 
Refs. I233II234I . 
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Fig. 62: Opposite sign dimuon mass spectra in the region < m < 5 GeV (left) and 5 < m < 15 GeV (right) for central 
Pb+Pb reactions. A pr threshold of 3 GeV /c is applied on each single muon. The upper and lower histogram correspond to 
the total combinatorial background and to the correlated signal from bottom decay, respectively. The spectra were obtained by 
means of fast simulations. From Ref. L236J . 



In addition, it is important to remark that the bottom cross section can be also determined from 
the like-sign dimuon yield. In fact, B meson oscillation as well as specific B meson decay chains 
represent a source of correlated like-sign muon pairs. This component can be extracted by subtracting 
the uncorrected signal from the like-sign spectrum II237I . 

Note finally that the open heavy flavor cross sections can also be extracted from the single muon 
spectra, as recently shown by the PHENIX collaboration with single electron spectra 11481 . 
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Table 26: Expected background and signal rates, signal-to-background ratios and significance in pp and in central Pb-l-Pb and 
Ca-l-Ca collisions for charmonium and bottomonium. The numbers correspond to an interval of ±l(j around the resonance 
mass. All rates and ratios are for a 10® s run. The luminosities for the Pb-fPb, Ca-Ca and pp system are 5 x 10^'', 1 x 10^^ and 
1 X 10^^ cm~^s~\ respectively. FromRef. |233|. 



system 


state 


B (xlO^) 


S (xl03) 


S/B 


S/VS + B 






320 


230 


0.72 


310 




-0' 

T 


150 


4.6 


0.03 


12 


Pb-i-Pb 


T 


0.25 


1.8 


7.1 


39 






0.22 


0.54 


2.5 


19 






0.18 


0.26 


1.5 


12 






760 


2000 


2.7 


1200 






360 


41 


0.11 


64 


Ca-i-Ca 


T 


1.5 


16 


10.7 


120 






1.3 


4.8 


3.6 


61 






1.1 


2.3 


2.0 


39 






64 


850 


13 


890 






32 


17 


0.53 


76 


VP 


T 


0.49 


6.5 


13 


78 






0.45 


2.0 


4.4 


40 






0.38 


0.94 


2.5 


26 



8.5 Electrons 

The electrons, measured in the central region, will also give access to heavy quarkonia together with 
open charm and open bottom from the dielectron continuum. The key detector is the TRD which pro- 
vides electron identification for pT > 1 GeV/c and an electron trigger for > 3 GeV/c. It is operated in 
conjunction with the other ALICE central detectors for precise momentum measurements. In addition, 
the vertex capabilities of the ITS can distinguish between primary and secondary J/t/j. While the identi- 
fication of the primary J/t/j is of crucial importance for QGP studies, the identification of secondary J/ijj 
can lead to a direct measurement of the B meson production cross section. Furthermore, single high-p^ 
electrons with displaced vertices give access to the inclusive heavy flavor cross sections. 

8.51 Quarkonia 

The J/t(j and T acceptances in the dielectron channel are displayed in Fig.|63] The effects of the TRD- 
Ll trigger, simulated by requiring a sharp px threshold of 3 GeV/c on both decay electrons, are shown 
in the middle and right panels of the figure. 

The effects of the pT cut on the T are similar to those observed for the J/V' in the dimuon channel, 
see Fig.|60l The cuts have a much stronger effect on the J/ip distribution. It is only at rather large px 
that the decay kinematics allow both electrons to have pT > S GeV/c. Therefore, while the T can 
be measured down to pT = 0, there is no J/t/j acceptance below pT ~ 5.2 GeV/c because of the 
trigger condition. As in the case of the dimuon channel, the expected mass resolution is of the order 
of 100 MeV/c^ for the T when the L3 magnet is operated at B = 0.4 T, allowing the T states to be 
separated. 

The expected statistics are strongly dependent on the trigger efficiency which decreases with par- 
ticle multiplicity. In the worst case anticipated for Pb-i-Pb, dNch/dr] = 8000, the expected number of T 
per month is ~ 2600 in minimum-bias collisions. This number is expected to be significantly larger for 
lower charged particle multiplicities 12381 . 
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The current J/ ijj simulations indicate that the statistics will be significantly lower than that achieved 
in the dimuon channel, mainly due to the high pT cut of the electron trigger. In a scenario with a smaller 
charged particle multiplicity, the electron trigger px threshold could be lowered, leading to better J /ip 
detection efficiency. 

8.52 Open Charm and Open Bottom 

The same studies described in section I02l for the muon channel will be also carried out in the electron 
channel using similar methods. In addition, the electron identification in the TRD can be used with the 
very good vertexing capabilities provided by the ITS. Due to the finite cr of D and B mesons, their decay 
electrons have impact parameters^' dO ~ 100 — 300 jivn and dO ~ 500 /um, respectively. Figurel64lshows 
the resolution of the track position at the vertex achieved in ALICE^^ (left) and the dO distributions for 
electrons coming from D and B mesons and from other parent particles (right). The expected rates (per 
unit rapidity at midrapidity) for semileptonic decays of D and B mesons in central Pb+Pb collisions (5% 
cJtot) at y/s = 5.5 TeV are 2 and 0.09 per event, respectively. 

The performance of ALICE for open bottom detection in the B ^ e + X channel has been studied 
using a detailed simulation of the apparatus and including all relevant background sources (direct-charm 
decays, pair production due to photon conversion in the detector materials, Dalitz decays of light and 
strange mesons, pions misidentified as electrons) I239II . A very pure electron sample is obtained using 
the transition radiation technique (TRD) in conjunction with the dE/dx measurement in the TPC. The 
fraction of pions misidentified as electrons is ~ 10^"^. The electrons from B decays have harder px 

^'The impact parameter dO is defined as the distance of closest approach of the particle trajectory, projected in the plane 
orthogonal to the beam axis, to the interaction vertex. 

^^To obtain the impact parameter resolution, the track position resolution must be quadratically combined with the error on 
the primary vertex position, negligible in Pb+Pb collisions, ~ 15 ^m. 
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spectra and, as already mentioned, broader dO distributions. Therefore, a sample of electrons from 
open bottom decays can be selected using suitable pT and dO thresholds. The left-hand side of Fig. |65] 
presents the signal/(signal+background) ratio for three values of the electron px cut as a function of the 
dO threshold. As an example, a condition of > 2 GeV/c and dO > 180 /^m gives a sample of 90 % 
purity. The expected number of B mesons detected in 10^ central Pb+Pb events (1 month of ALICE 
data taking) is ~ 80000 without any specific electron trigger. A more complete view of the attainable 
statistics is shown on the right-hand side of Fig.|65]as a function of the dO threshold for three different 
Pt thresholds. 

If the Pt and dO thresholds are lowered to ~ 1 GeV/c and ~ 100 /im, the contribution of electrons 
from open charm decays becomes comparable with that of open beauty. A method for disentangling the 
two contributions is currently under study. 

The second way to measure the open bottom cross section is based on the measurement of sec- 
ondary J/?/;'s from B decays. These secondary J/ip's are produced at large distances from the primary 
vertex and can therefore be selected by identifying e'^e~ pairs with displaced vertices, as shown on the 
left-hand side of Fig. |66l The resulting invariant mass spectrum is presented on the right-hand side of 
Fig.lMl 




transverse momentum (GeV/c) dO (|IITl) 

Fig. 64: Left: the transverse momentum dependence of the bending projection (dO) of the impact parameter resolution. Right: 
the dO distribution of electrons originating from different parent particles. From Ref. 12401 . 

8.6 Electron-Muon Coincidences 

The correlated cc and bb cross sections can be measured in ALICE using opposite sign electron-muon 
pairs. The electron is identified in the central part and the muon is detected in the forward muon spec- 
trometer. This channel is the only leptonic channel that gives direct access to correlated cc and bb pairs. 
Indeed, in contrast to the e+e^ and fi^iJ.^ channels, resonance decays, direct dilepton production and 
thermal production cannot produce correlated e/i pairs. Within ALICE, the e/i channel has the additional 
advantage that the rapidity distribution of the corresponding signal extends in 1 < y < 3, bridging the 
central and the forward acceptances 11371 . Measurements of e/x coincidences were successful in pp in- 
teractions at ^/s = 60 GeV [241 1 and in pN interactions at -y/i = 29 GeV [242|. They are planned for 
heavy ion collisions with the PHENIX detector at RHIC 1243.1 . Preliminary simulations have shown that 
such measurements might be possible with ALICE 12401 12441 1^?51 . 
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do threshold [ ^m] dO threshold [^m] 



Fig. 65: Detection of B — > e + X: S/(S + B) ratio (left) and signal statistics for 10^ central Pb+Pb events (right) as a function 
of the dO threshold and for three values of the pr threshold. From Ref. f239|. 



8.7 Hadrons 

In the central part of ALICE, heavy mesons and baryons can be fully reconstructed from their charged 
particle decay products in the ITS, TPC and TOP. Thus, not only their integrated yields, but also their 
Pt distributions can be measured. The most promising open charm decay channel is — > K~tt'^ (and 
its charge conjugate) with a branching ratio of ~ 3.9% and cr = 124 fim. The expected rates (per unit 
rapidity at midrapidity) for (and D^) mesons, decaying into K^ir^ pairs, in the 5% most central 
Pb+Pb coUisions at ^/s = 5.5 TeV and in pp collisions at ^/s = 14 TeV are 5.3 x 10^^ and 7.5 x 10~^ 
per event, respectively. 

This decay channel allows direct identification of the particles by computing the invariant 
mass of fully-reconstructed topologies originating from displaced secondary vertices. The left-hand side 
of Fig. [^sketches the decay. The main feature of this topology is the presence of two tracks with impact 
parameters of order 100 /xm. 

The capabilities of the ALICE central barrel for D*^ -meson searches have been investigated both 
for Pb-i-Pb collisions, see Refs. 12241 1246. ,247.1 . and for pp collisions. All the relevant background 
sources have been included in these studies, along with a detailed simulation of the detector response. 
In particular, for pp, the position of the interaction vertex will have to be reconstructed event-by-event 
using the tracks^^ . 

In Table |27] the statistics for signal, S, and background, B, per event and the signal-to-background 
ratio after reconstruction are presented. The only selection is on pair invariant mass, | Mktt — M^o \ < 
3 cr = 36 MeV. 

The most effective selection in order to extract the charm signal out of the large combinatorial 
background of opposite-charge track pairs is based on the requirement to have two tracks with large 
impact parameters and a good pointing of the reconstructed momentum to the collision point {i.e. 
the pointing angle @p between the momentum and its flight-line should be close to 0, as shown on 
the left-hand side of Fig.l67b. The selection strategy is described in detail in Refs. 12241 12461 1247II . In 
Table 128) the statistics after selection, together with the signal-to-background ratio and the significance, 
are presented. The expected statistics are ~ 13000 reconstructed in 10^ central Pb-i-Pb events and 

^^^The transverse size of the LHC Pb beams will be only 15 ^m. In the pp runs, the beams will be defocused at the ALICE 
interaction point in order to bring the luminosity down to ~ 10^'^ cm^'^s"^, see Section lS^ Thus, the transverse size of the 
beams may be increased up to ~ 100 — 200 fim. The resolution of the interaction vertex position in the transverse plane 
achieved using the tracks is ~ 50 /xm. 
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Fig. 66: Left: dependence of the J/ip signal on the distance to the primary vertex. Right: invariant mass distribution of pairs 
with displaced vertices. From Ref. 12401 . 

Table 27: Statistics for signal and background and signal-to-background ratio including only selection on pair invariant mass, 
\Mk-^ - Moo\ < 3 cr = 36 MeV. 



system {^/s) 


S/event 


B/event 


S/B 


Vh+Vh (5.5 TeV) 
pp (14 TeV) 


1.3 X 10"^ 

2.4 X 10"'' 


2.8 X 10^ 
1.1 X 10-1 


4.5 X 10"^ 
2.3 X 10-3 



~ 19000 in 10^ pp events. The (relatively) poorer quality of the pp results is due to the larger error on 
the primary vertex position with respect to Pb-i-Pb. 

Table 28: Statistics for signal and background, signal-to-background ratio and significance after selections. The significance is 
computed for the number of events collected in a one month run, 10 for Pb-FPb and 10^ for pp. 



system {y/s) 


S/event 


B/event 


S/B (%) 


S/VS + B 


Pb-FPb (5.5 TeV) 


1.3 X 10"-^ 


1.2 X 10"^ 


11 


37 (lO'' events) 


pp (14 TeV) 


1.9 X 10"^ 


1.7 X 10-4 


11 


44 (10^ events) 



A typical example of the D candidates invariant mass distribution in central Pb-i-Pb reactions, 
after selections, is shown on the right-hand side of Fig, [^before and after background subtraction. 

Figures |68l and |69l show the pT distribution of the signal and the background (normalized to one 
event) and the corresponding significance in Pb-i-Pb and pp collisions respectively. Note that for px > 
7 — 8 GeV/c, even if the simulated background statistics are insufficient, the significance is dominated 
by the signal statistics. The significance is larger than 10 for up to pT ~ 10 GeV/c, both in Pb-i-Pb 
and in pp collisions, if a 1 GeV/c pr bin size is used. As shown in Ref. [2471, the significance at low 
PT can be further improved by selecting sub-samples of candidates for which the kaon is identified. 
The values obtained with this method are shown with different symbols in Figs. I68I69I In this case the 
production cross section can be measured down to ~ 1 GeV/c in Pb-i-Pb collisions and down to 
almost pt = Oin pp collisions. 
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Invariant Mass [GeV] 
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Invariant Mass [GeV] 

Fig. 67: Left: schematic representation of tlie Ktv decay showing the impact parameters of the kaon (dOif) and of the 

pion (dO,r) and the pointing angle (6p). Right: Ktv invariant mass distribution for lO'^ central Pb+Pb events. The background- 
subtracted distribution is shown in the inset. From Ref. |247|. 
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Fig. 68: Left: the transverse momentum distribution of selected Z)° candidates in Pb+Pb collisions at y/s — 5.5 TeV. The 
signal and background are normalized to one central Pb+Pb event. Right: the corresponding significance for 10^ events. At 
low pt the values obtained requiring kaon identification are shown by the stars. 



8.8 Proton-nucleus Interactions 

Proton-nucleus collisions are an important part of the ALICE physics program [2241. The main motiva- 
tion is to bridge the results obtained in pp and AA interactions to unravel initial and final state medium 
effects. 

More precisely, pA studies of heavy flavour production will provide essential measurements of 
cold nuclear matter effects such as nuclear shadowing effects on heavy quarks I177II and nuclear ab- 
sorption of quarkonium. Nuclear absorption effects have been demonstrated to be of crucial impor- 
tance for disentangling "normal nuclear absorption" from "anomalous suppression" of quarkonium at 
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Fig. 69: Left: the transverse momentum distribution of selected candidates in pp collisions at ^/s — 14 TeV. The signal 
and background are normalized to one minimum bias pp event. Right: the corresponding significance for 10^ events. At low 
Pt the values obtained requiring kaon identification are shown by the stars. 



the SPS r225' '2261 12271 . Similarly, the pA open charm and bottom yields will provide valuable ref- 
erences for the onset of thermal charm production in nuclear collisions and for the modification of the 
heavy flavour spectra due to energy loss in the QGP. 

The ALICE pA program, foreseen as dedicated pPh, dPh or aPb runs, is expected during the first 
few years of LHC operation 12481 . As in pp and A A collisions, quarkonia and heavy flavor mesons 
will be measured through their (semi-)electronic decay channels in the ALICE central barrel. Also, D 
mesons will be fully reconstructed using exclusive hadronic decay channels. The detector performance 
for K^-K^ (and charge conjugate) decays in pPb collisions is expected to be similar to that of 

pp since the multiplicity of background tracks will not be much larger than in pp interactions. The 
asymmetric muon spectrometer allows all these studies in both direct and inverse kinematics since both 
pA and "Ap" operation is, in principle, possible. Detailed simulations investigating the expected detector 
performance in pA collisions are underway. 

8.9 Run Plan 

The present ALICE plans for data taking during the LHC heavy ion runs, in addition to participation in 
the 14 TeV pp runs, are now briefly described, see Ref. |248 1 for more details. In the first 5-6 years of ion 
runs, we propose: 1-2 years Pb-i-Pb; 1-2 years Ar-i-Ar; and 1 year pPb-like collisions (pPh, dPh or aPb). 
Later options depend on the outcome of the early runs: dedicated pp or pp-like {dd or aa) collisions; 
other AA systems such as O-i-O, Kr-i-Kr or Sn-i-Sn; other d{a)A systems; or a low energy Pb-i-Pb run. 

The maximum Pb-i-Pb luminosity for the ALICE central detectors is L = 10^'' cm~^s~^, imposed 
by the pile-up event rate in the TPC and the DAQ bandwidth, 1 kHz, of the muon spectrometer. Never- 
theless, to achieve better muon statistics, dedicated high luminosity runs of up to L = 2 x 10^^ cm^^s^^ 
involving the muon spectrometer without the TPC could be foreseen with a high pT muon trigger, pri- 
marily for T decays. 

To study the results as a function of energy density, at least one intermediate mass system, most 
likely Ar-i-Ar, will be studied in the first years. If more energy density points are needed, other systems 
such as Sn-i-Sn, Kr-i-Kr and O-i-O may be studied at a later stage. Luminosities similar to those used for 



97 



the Pb+Pb runs will be used for the intermediate mass ion runs. 

Although ALICE will take pp data at 14 TeV in parallel with the other experiments, the luminosity 
must be reduced to 3 x lO^'^ cm~^s~^ to keep the pile up rate in the TPC and the ITS at an acceptable 
level. Higher luminosity runs with the muon spectrometer alone are possible with luminosities up to 

5 X 10^^ cm~^s^^. 

In addition, lower energy pp collisions, closer to the Pb+Pb center of mass energy, 5.5 — 7 TeV, 
might be needed for reference data. 

8.10 Summary 

The capabiUties of ALICE for heavy flavour physics at the LHC will provide a comprehensive under- 
standing of open and hidden heavy flavour production at low x where strong nuclear gluon shadowing 
is expected. The large number of accessible channels will permit detailed investigations of the proper- 
ties of the nuclear medium. Open heavy flavour measurements will provide information on the primary 
production mechanisms, parton energy loss and possible secondary production of heavy quarks. The 
complete spectroscopy of the quarkonium states will be used to probe the deconfined medium and to pin 
down its temperature. The simultaneous measurement of open and hidden heavy flavours will provide a 
powerful tool to study quarkonium suppression/enhancement in the QGP. Thanks to the large acceptance 
of the detector, the signals will be reconstructed over a broad transverse momentum range. This is par- 
ticularly important for establishing the reference production rates with high accuracy and to investigate 
the characteristics of the QGP. The measurements will be performed for a variety of heavy and light 
nucleus-nucleus systems as well as for pp and pA reactions. The latter are mandatory since, besides 
probing shadowing effects, they will provide the normaUzation for understanding the properties of the 
nuclear medium in A A collisions. 

In addition to the measurements discussed here, further exciting possibilities should be opened 
with, for example, the reconstruction of hadronic decays of charged D and B mesons, reconstruction of 
double charm baryons by combining hadronic and leptonic chaimels, and dilepton measurements at very 
high invariant mass. 
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9. QUARKONIA AND HEAVY FLAVOUR DETECTION WITH THE CMS DETECTOR 

9.1 Introduction 

The interest in quarkonium production at the LHC emerged from the CERN SPS results 12491 12251 
which showed a strong anomalous suppression in J /ip production in Pb+Pb collisions at ^/s = 17.3 GeV. 
RHIC will study J/ production and suppression in detail at y/s = 200 GeV since the temperature of the 
system produced at RHIC should be high enough for direct J /ip suppression. Although the T production 
cross section is large enough to be observed at RHIC, albeit with limited statistics, its suppression is not 
expected until the higher initial temperatures at the LHC are reached. CMS is particularly well suited to 
study the T family, the continuum up to the mass and, to a lesser extent, the J and ip' . 

Open heavy flavour production is important for studying the behaviour of massive colour charges 
in a dense medium at the LHC. The charm and bottom production cross sections are much larger than 
at RHIC. Systematics studies of heavy flavours, especially B mesons, can be performed with CMS. 
In particular, in-medium gluon radiation and collisional energy loss of heavy quarks can result in the 
suppression and modification of the high-mass dilepton 111 371 12501 125 II and secondary B J /ijj decay 
spectra I 250ll25ll . Finite quark mass effects can suppress medium-induced radiation of heavy quarks, 
enhancing e.g. the iJ/vr ratio [147|. 

We briefly describe the CMS detector in section |9^ Then we provide the resonance cross sections 
and the background description used in our simulations in section 1931 Section l94l describes the detector 
response, the acceptance and the dimuon reconstruction algorithm. The results of the simulations are 
presented in section l931 including the signal/background ratios and the dimuon invariant mass distribu- 
tions. In the last section, 19.61 the CMS capabilities for open heavy flavour and measurements in the 
mass range above the T are discussed. 

9.2 The CMS Detector 

The CMS detector is designed to identify and measure muons, electrons, photons and jets over a large 
energy and rapidity range. It offers the widest muon acceptance centered at midrapidity. CMS is par- 
ticularly well suited to study the T family, the continuum up to the Z'^ mass and, to a lesser extent, 
the J/ijj and ip' . The CMS dilepton capability allows systematic studies of heavy flavour physics. The 
impact parameter (centrality) of the collision can be determined from measurements of transverse energy 
production over the range |r/| < 5. 

CMS has a high-field solenoidal magnet with a uniform 4 T field, leading to a compact detec- 
tor. The first absorber, the electromagnetic calorimeter, is 1.3 m from the interaction point, eliminating a 
large fraction of the hadronic background. The powerful tracking system provides good track reconstruc- 
tion efficiency for dimuons even for very large charged particle densities, dA'ch /dy 8000, providing 
excellent dimuon mass resolution. 

A detailed description of the detector elements can be found in the corresponding Technical Design 
Reports I252ll253ll254ll255"l . A longitudinal view of the detector is shown in Fig.^Ol The central element 
of CMS is the magnet, a 13 m long, 6 m diameter, high-field solenoid with an internal radius of Ri3 m. 
The tracker and muon chambers cover the pseudorapidity region \rj\ < 2.4, while the ECAL and HCAL 
calorimeters reach \rj\ = 3. A pair of quartz-fiber very forward (HE) calorimeters, located ±11 m from 
the interaction point, cover the region 3 < |t/| < 5 and complement the energy measurement. The tracker 
is composed of pixel layers and silicon strip counters. The barrel part of tracker, |t/| < 0.8, consists of 
three pixel layers 4, 7, and 1 1 cm radially from the beam line and 9 layers of silicon strip counters. The 
endcap part of tracker, rj > 0.8, is more complex. There are two pairs of pixel layers symmetrically 
located 34 cm and 45 cm from the geometrical center of the detector. Outside this, two symmetric sets 
of 12 silicon layers are situated up to 265 cm from the center of the detector. The electromagnetic 

^''Authors: M. Bedjidian, O. L. Kodolova, R. Kvatadze, I. P. Lokhtin. 
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Fig. 70: A longitudinal view of the CMS detector. 



calorimeter is made of almost 83000 scintillating PbW04 crystals. The hadronic calorimeter consists of 
scintillator inserted between copper absorber plates. 

The CMS muon stations consist of drift tube chambers (DT) in the barrel region (MB), |ry| < 1.2, 
cathode strip chambers (CSCs) in the endcap regions (ME), 0.9 < I??! < 2.4, and resistive plate chambers 
(RPCs) in both barrel and endcaps, for |t/| < 1.6. The RPC detector is dedicated to triggering, while the 
DT and CSC detectors, used for precise momentum measurements, also have the capability to self-trigger 
up to \ri\ < 2.1. The muon system can thus reconstruct muons in the range \r]\ < 2.4. 

In each of the 12 sectors of the barrel return yoke of the magnet (YB), covering 30° in azimuthal 
angle (/), there are four muon barrel 'stations', MBi, i = 1,4 in Fig. E0| Each station consists of one 
DT and one (or two, in the two innermost stations) RPC chambers. In each DT chamber, there are two 
"superlayers" of four DT planes measuring the (r, 0) coordinates in the bending plane and another super- 
layer measuring the z coordinate. In the outermost station, MB4, there are only two (r, 0) superlayers. 
Thus a muon traversing the whole barrel has 44 measurement points ('hits') in the DT systems and 6 hits 
in the RPC. There are a total of 240 stations in the barrel system. Each endcap region has four muon 
stations of CSCs. Each chamber gives six measurements of the 0-coordinate (strips) and six measure- 
ments of the r-coordinate (wires). RPCs will be added in the barrel and endcaps to provide an additional 
trigger. Six layers of RPCs will be mounted in the barrel chambers and four in the endcaps. 

9.3 AA Collisions 

9.31 Resonance cross sections 

We use the quarkonium cross sections per nucleon calculated in section IT2l in our simulations. The cal- 
culations, performed in the CEM, include nuclear shadowing but do not include absorption by nucleons 
and secondaries since the interplay between shadowing and absorption as a function of ^/s is unknown. 
For convenience, the inclusive A A dimuon cross sections given in Tables 9 and 10 are reproduced in 
Table |29l 
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Table 29; Total cross sections for quarkonium production in minimum bias AA collisions. 



System 


Pb+Pb 


Sn+Sn 


Kr+Kr 


Ar+Ar 


^/s (TeV/nucleon) 


5.5 


5.84 


6.14 


6.3 


Resonance 


B^f, OAA (/^b) 


J/i, 


48930 


17545 


9327 


2321 


^' 


879 


315 


167.6 


41.7 


T 


304 


108.1 


57.4 


13.8 




78.8 


28 


14.8 


3.6 


^11 


44.4 


15.8 


8.4 


2.0 



932 Background and kinematic distributions 

Soft hadronic background. Soft hadron production is the main dimuon background. This background 
is estimated from the charged particle multiplicity and the shape of the pion and kaon transverse momen- 
tum distributions. 

The charged particle multiplicity, dNc\i/dri, is unknown for the LHC energies. Early estimates 
of this multiplicity in central collisions (6 = 0) and midrapidity (ry = 0) were as high as 8000 for 
Pb+Pbcollisions |256|. In contrast, the first extrapolations from RHIC suggested that the Pb+Pb mul- 
tiplicity could be as low as 2000 1165J. We choose two values for our simulations: 2500 (low) and 
5000 (high). Lower multiplicities reduce the dimuon background. The track reconstruction algorithm 
was designed to function for a multiplicity of 8000, so that the current, lower, estimates should result in 
higher reconstruction efficiencies. The high and low values of the multiplicities used in the simulations 
are shown in Table l30l The upper numbers for each set correspond to the 5% most central collisions, 
assuming the multiplicity scales as A^^^ from Pb-i-Pb. The lower numbers for each set correspond to 
collisions at 6 = 0. In this case, the multiplicity can be higher. 

Table 30: The high and low multiplicity sets, dNch/dr]\rf^o, used in the simulations. In both cases, the results are shows for 
the 5% most central collisions and impact parameter 6 = 0. 





dNc\,/dri\r,=Q 


Set 




Vh+Vh 


Sn-FSn 


Kr+Kr 


Ar+Ar 


High 


5% most central 


5000 


2400 


1500 


550 




6 = 


7500 


4000 


2600 


850 


Low 


5% most central 


2500 


1200 


750 


280 




6 = 


3800 


1900 


1300 


450 



The soft hadronic background is assumed to entirely consist of pions and kaons produced accord- 
ing to the ratio K/tt = 0.12, independent of 6. We use the SHAKER algorithm [257] to simulate the 
pion and kaon momentum distributions where {pt)-k = 0.47 GeV/c and {pt)k = 0.67 GeV/c. We 
note that HIJING 1 215 216 1 produces pions and kaons with lower momenta, {pt)tt = 0.43 GeV/c and 
{pt)k = 0.52 GeV/c. Due to the magnetic field, only muons with pT > 3.5 GeV/c will reach the barrel 
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muon chambers. This minimum momentum requirement serves as the primary background rejection cri- 
teria. More than 99.8% of the SHAKER hadron background is rejected. The lower average momentum 
of HIJING reduces the probability for acceptance of decay muons by a factor of two with respect to 
SHAKER. Track reconstruction will also be an effective method of background rejection. 

Heavy flavour decays. The other important source of quarkonium background is QQ pair production 
and decay. The number of pairs produced in AA collisions as a function of impact parameter b is 

N{QQ) = a{QQ)TAA{b) , 

directly proportional to the nuclear overlap function r^^(6) where TpbPb(O) = 30.4/mb, TsnSn(O) = 
13.3/mb, TKrKr(O) = 8.5/mb and rArAr(O) = 2.9/mb. The NLO QQ production cross sections in 5.5 
TeVpp interactions are a{cc) = 6.3 mb and a{hb) = 0.19 mb (with large uncertainties - see section 121. 
Since the center of mass energies are higher for lighter ions, the QQ cross sections per nucleon increase 
slightly from Pb+Pb to Ar+Ar interactions. The number of QQ pairs at 6 = for the collision systems 
considered is shown in Table 1311 



Table 31: Number of QQ pairs produced at 6=0. 





Pb+Pb 


Sn+Sn 


Kr+Kr 


Ar+Ar 


cc 


192 


86 


57 


20 


bb 


6 


3 


2 


1 



Although N{cc) » N{bb), dimuons from cc decays are not the largest contribution to the back- 
ground. Among all background sources taken into account, bb decays produce muons with the highest 
average transverse momentum, making a significant contribution to the background. The muon pT and 
r] distributions are extracted from PYTHIA UTTl with {pi^)c = 0.51 GeV/c and {pi^)b = 1.2 GeV/c. 

Only 6% of single b decays produce two muons and, of these, 2/3 are of opposite sign. When the 
CMS kinematic cuts, pt^ > 3.5 GeV/c and \ri\ < 2.4, are included, more than 98% of the single b decay 
muons do not survive the cuts. The probability of obtaining a dimuon from these decays after cuts is 
« 10~^, shown in Table 1^ About 80% of the accepted pairs are of opposite sign so that the correlated 
like-sign muon rate is negligible. Thus, to simplify the simulations, we assume that all muon pairs from 
bb decays are of opposite sign. 



Table 32: The probability of n-muon signals in b quark decays. 





fi 


1/i 


2 fj, 


3 fi 


before CMS cuts 


0.626 


0.308 


0.060 


0.006 


after CMS cuts 


0.987 


0.0125 


0.00018 






9.4 Detector Response 

9.41 Acceptances 

The pt and rj distributions of the quarkonium sstates used in the simulations were obtained from several 
sources. The pT distributions of the three T(S') states were extrapolated from CDF data I258II while 
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PYTHIA was used to generate the t] distributions. PYTHIA was also used to obtain the J/tp and Tp' 
distributions. 

We assume that the primary vertex is at the geometrical centre of the detector. Each resonance 
decay muon is tracked through a 4 T magnetic field by GEANT 12591 using the CMSIM package [260 1. 
A dimuon is accepted if each decay muon of the pair passes through at least one muon chamber. We 
also introduce a px cut of 3.5 GeV/c for each muon, just above the minimum pi^ needed to reach the first 
muon chambers in the barrel. This cut on pf^ leads to very different resonance acceptances. The J/ijj 
decay muons have (p^) = 1.7 GeV/c so that only 0.8% of the J/i/j decays give leptons above the muon 
threshold. On the other hand, the T decay muons have (pj,) = 4.6 GeV/c . Thus 40% survive the pT cut. 
The acceptances discussed below and shown in Fig.^Jare geometrical only. The dimuon reconstruction 
efficiencies, the trigger efficiency and the muon cut are not taken into account. 

The T acceptance is fiat, 27% in the full t] range of the detector before the muon pT cut, 
dropping to 16% once the muon pT cut is applied. With the barrel alone, the T acceptance is 13% before 
the pt cut and 9.5% after. The px dependence of the accepted T in Fig.^Jdoes not change between the 
full detector and the barrel. The p-r distribution starts from zero and is only statistically limited in the 
high-pT region, allowing a comparison between the high- and \ow-pT events. 

The situation is different for the J/ip. The integrated J/ip acceptance is about ten times lower 
than the T acceptance. There is an additional factor of 10 reduction when the J/ip acceptance is limited 
to the barrel due to the natural pt^ cut provided by the material located in front of the muon chambers. 
Including a 3.5 GeV/c cut diminishes the acceptance very strongly. Only muons from high px J/ip 
decays (j)t J/ip i 4 GeV/c) can reach the barrel muon chambers, as clearly shown in Fig.^J However, 
when the endcap regions are included can the J/ ip be detected over the whole pT range. The acceptance 
is small but the expected rate is large. 




p^'' (GeV/c) 

Fig. 71: The pT dependence of the dimuon acceptance for J/ip (upper) and T (lower) resonances. In each case 
both the full CMS detector (solid line) and the barrel only (dashed line) have been considered. 

9.42 Track reconstruction 

Pattern recognition and track reconstruction determine, to a large extent, the feasibility of heavy ion mea- 
surements with CMS. The particles were tracked through the detector using a full GEANT simulation. 
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The tracks where followed down to 100 KeV for electromagnetic processes and 1 MeV for hadronic pro- 
cesses, an order of magnitude below the default GEANT energy threshold cuts. The quarkonium studies 
have has been carried out for a range of dNch/dr], up to an extreme value of 8000, far beyond the present 
upper limit of 5000 for central Pb+Pb collisions. The upper limit of dNch/drj = 8000 led to occupancies 
of up to 18% for the 5 outermost Si-layers of the barrel, located more than 70 cm from the beam, and of 
50% or more for the innermost Si-strip detectors. Occupancies of 2-4% 126 11 are expected for the pixel 
layers due to their very high granularity. Thus, the occupancies determine the number of layers used in 
the track-reconstruction algorithm. 

The goal of the track-finding algorithm is to select pairs of muon tracks coming from the interac- 
tion point. In the barrel, the algorithm uses the two pixel layers at 7 and 1 1 cm and the four outermost 
Si-layers. This is the minimum combination necessary to get good efficiency and purity for high occu- 
pancy rates in the tracker. The usage of the third pixel layer at 4 cm is under investigation. 

The vertex is determined from a combination of clusters in the pixel layers assuming the interaction 
occurs at the geometrical centre of the detector, (0, 0, 0). Figure ^2] shows the results of the interaction 
point vertex determination for a simulated collision at z = 5 cm for dNch/drj = 2500 and 8000. The 
resolution is az = 140 /xm. 

The track finding algorithm starts from the muon chambers since they are the least populated 
planes. The tracks are propagated within roads in transverse {R/z, (p) and longitudinal (R, z) planes with 
simple parameterizations d(p = k dr/p^i for the barrel, dcf) = k dz/pi^ for the endcaps, and a straight line 
in the (i?, z) plane. 

The track candidates are then fitted and, after vertex constraints, the best quality track is selected. 
Several dimuon events have been studied: homogeneous dimuons from T decays as well as mixed pairs 
originating from t: /K/h decays where one muon comes from a vr or K decay and the other from a b 
decay. For the moment, only barrel tracks have been investigated. In the region < 0.8, the efficiency 
of T reconstruction is « 90% for dNch/dr] = 2500 and 85% for 5000 with a purity above 99%. The T 
reconstruction efficiency is 76% for dNch/dr] = 8000. The efficiency for vr and K decays is about six 
times lower, leading to good background rejection. For 0.8 < \r]\ < 1.3, the T efficiency is 50%. The 
high precision tracking leads to an excellent dimuon mass resolution of 46 MeV at the T peak when both 
tracks have \r]\ < 0.8, increasing to 60 MeV when one track has 0.8 < \r]\ < 1.3. 

9.5 Invariant Mass Spectra 

9.51 Signal/background and mass distributions 

The simulations are based on the weight method with acceptance tables. They consider combinations of 
muons with pj, > 3.5 GeV/c. Either both muons originate from vr and/or K decays, hh, or from cc and 
bb decays, or from a "mixed source", one from it or K decays and one from heavy quark decays, hb and 
he. 

The invariant mass of each pair is calculated. Two mass regions are then defined, the T region, 
8.5 < M^^ < 11 GeV/c^ and the J/ip region, 2 < M^^ < 4.5 GeV/c^ . The resulting invariant mass 
distribution is smoothed by a Gaussian with a width determined by the expected mass resolution. 

The opposite sign dimuon invariant mass distributions obtained in Pb-i-Pb, Sn-i-Sn, Kr-i-Kr and 
Ar-i-Ar collisions are shown in Fig. ^31 for the J/t/j and in Fig. ^4] for the T. Using the high and low 
multiplicity sets given in Table l30l we present results for a one month run, 1.3 x 10^ seconds, assuming 
a 50% machine efficiency. The J/^p and T statistics are given in the figures and summarized in Table l33l 
assuming the average luminosities of Ref. 1.262 1 . The 18000 T measured in a one month Pb-i-Pb run are 
statistically significant enough to bin with impact parameter. Limiting detection to the barrel reduces the 
statistics, but not unacceptably for intermediate-mass ions. 

The signal-to-background ratio, S/B, between the number of detected resonances, J/ip or T(1S), 
and the number of opposite sign dimuons under the peak in the mass window Mi.es±50 MeV, is calculated 
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z (cm) z (cm) 

Fig. 72: Distribution of the z position of the primary vertex for dNch/drj of (a) 2500 and (b) 8000. The initial 
collision was generated at x = y = 5 and z = 5 cm. 



from Fig.l73]for the J /ij; and from Fig.E4]for the T. This ratio depends on the multiplicity. Tablel34lgives 
the S/B ratio for the J/ifj and T and the significances, S/VS + B, for the three T{S) states. The results 
are shown for both the high and low multiplicity sets of Table BUI for each collision system. The lower 
signal-to-background ratios and significances correspond to the larger background, high multiplicity set. 
The higher signal-to-background ratios and significances result from the reduced background of the low 
multiplicity set. 

Table 33: The number of resonances detected in a one month, 1.3 x 10^ s, run in the full detector and in the barrel alone (in 
parentheses). 



iV,es/103 


Vh+?h 


Sn-FSn 


Kr-FKr 


Ar-i-Ar 




24 


160 


360 


1450 




(4.7) 


(31) 


(81) 


(275) 


T 


18 


116 


260 


1020 




(7.9) 


(51) 


(114) 


(450) 




5.4 


35 


80 


310 




3.4 


22 


50 


195 



Figure corresponding to Pb-i-Pb collisions with the high multiplicity set, shows the expected 
opposite sign dimuon mass distribution in the T mass region in more detail. The most important con- 
tributions to the background are given. The background is typically dominated by vr and K decays in 
combination with each other, hh, or with a heavy quark decay, hb and he. 

Although N{bb) <^ N{cc), in an AA collision the bb contribution is more important. The back- 
ground is primarily from uncorrelated muon pairs. This background can be subtracted from the opposite 
sign spectrum using the 'like-sign' spectra, as done by the SPS NA50 collaboration 12491 12251 . In this 
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Table 34: The signal-to-background ratios, S/B, and the significances, S/y/S + B in the 5% most central collisions for the two 
multiplicity sets given in Table l30l The smaller numbers correspond to the high multiplicity set while the larger numbers are 
obtained with the low multiplicity set. 





Vh+Vh 


Sn-FSn 


Kr+Kr 


Ar-i-Ar 


S/B 




0.2-0.5 


0.4-1.1 


0.7-1.8 


2.0-6.8 


T 


0.4-0.9 


0.7-1.9 


1.5-4.3 


5.3-15.6 




T 


69-93 


220-276 


396-460 


925-978 


S/VS + B 




24-38 


84-123 


165-218 


447-512 




16-26 


55-86 


113-157 


325-391 



case, the signal is 

jySig ^ _ 2^/N++ N— 

where , A^++ and are the combinations of opposite sign pairs and positive and negative hke- 
sign pairs respectively in a given mass interval. Figure EH presents the mass distributions resulting from 
this subtraction. 




Fig. 73: The opposite sign dimuon invariant mass distributions in the J /ip mass region obtained in the 5% most 
central collisions in a one month run with the high multiplicity set (left) and low multiplicity set (right) as defined 
in Table|30| 



9.6 Open Heavy Flavours 

We briefly discuss some preliminary open heavy flavour results in CMS. The heavy flavour production 
cross sections in minimum bias nucleus-nucleus collisions, as well as those of Drell-Yan and Z", W, 
WW, WZ^ and Z'^Z'^ decays, were obtained from those in interactions at -y/s = 5.5 TeV assuming 
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Fig. 74: The opposite sign dimuon invariant mass distributions in the T mass region obtained in the 5% most 
central collisions in a one month run with the high multiplicity set (left) and low multiplicity set (right) as defined 
in TableEOl 
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Opposite sign dimuon invariant mass (GeV/c^) 



Fig. 75: The invariant mass distribution in the T mass region of opposite sign dimuons, including the largest 
background contributions, with the high multiplicity set. For hh, both muons come from tt and K decays. For hb 
and he, one muon is from tt or if decay while the second is from a 6 or c decay. For bb, both muons come from bb 
decays. 



aAA = -^(^pp- The pp cross sections were evaluated using PYTHIA with the CTEQ2L parton distri- 
bution functions. Since PYTHIA uses leading order matrix elements, a K factor of two was included 
for cc and bb production. Uncorrelated muon pairs from hadronic -k/K decays in minimum bias Pb+Pb 
collisions were obtained from HIJING. Figure ^Jpresents the invariant dimuon mass spectra for pj,>b 
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Fig. 76: The signal invariant mass distribution, after background subtraction, in J/ tp (left) and T (right) mass 
regions obtained in one month with the high multiplicity in Pb+Pb collisions. 




Fig. 77: Invariant mass distribution of fi^fi pairs for muons with pr > 5 GeV/c in one month Pb+Pb run. 



GeV/c and |?7| < 2.4, corresponding to a one month Pb+Pb run, 1.3 x 10^ s, assuming emachine = 0.5 
and {Caa) = 10^^ cm^^s^^. The dimuon detection efficiencies have been taken into account. In the 
mass range 10 < M^j^ < 70 GeV/c^ the dominant contribution comes from 66 fragmentation |217|. 
Hence, this mass region can be used to estimate the in-medium bottom quark energy loss in heavy ion 
collisions 1 137 .. .250. .25 1 1. The mixed-origin contribution, when one muon is from b ^ B fragmentation 
and the other from tt/ K decays, is about 16%. The contribution from cc fragmentation and vr/K decays 
are 6% and 5% respectively. Note that the charm energy loss could be significantly larger than the bottom 
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loss due to the lower charm quark mass, resulting in an additional suppression of the cc —>■ fi^fJ,'' yield. 
A clear signal from —>■ /i^/i^ decays, ~ 11000 events within Mz ± 10 GeV/c^, has a background 
of less than 5%. Since production is unaffected by final state interactions and is assumed to be pro- 
portional to the number of nucleon-nucleon collisions, it can be used as reference process to normalize 
the jet, quarkonium and heavy flavour rates in AA collisions relative to pp collisions. The dimuon rates 
from other massive sources (ti, WW, WZ'^, Z^Z^) are negligible. 

If we consider dimuons from bb decays as a signal, two kinds of background can be extracted: 
uncorrelated and correlated. The uncorrelated background, random decays of pions and kaons and muon 
pairs of mixed origin, can be estimated from the like-sign events and the signal. 

The main correlated background, Drell-Yan production, is unaffected by medium-induced final 
state interactions. These dimuons come directly from the primary nuclear interaction vertex while the 
dimuons from B and D meson decays appear at secondary vertices some distance away from the primary 
vertex. The path length between the primary and secondary vertices is determined by the lifetime and 
Lorentz boosts. A good way to discriminate the dimuons from B mesons from those emitted at the 
primary vertex can be the transverse distance, 6r, defined below. If Pmin is defined as a track point 
with minimal distance to the beam axis, z, then 6r is the distance in the x — y plane between points 
-Pi, min and P2, min belonging to different muon tracks. Muon pairs from bb decays show a rather flat 
distribution 1 25(7, '251'] while those from Drell-Yan production sharply decrease, vanishing at 6r = 70 
/im, see Fig. ED We estimate the accuracy of the track position determination to be ax = cTj^ ~ 10 fim 
and az ~ 100 ^m while the accuracy of the nuclear interaction point determination is ct^: = cjy ~ cJj; ~ 
20 //m for the CMS tracker |256|. We find that, for such a simple simulation of the tracker resolution, a 
5r > 50 fj.m cut suppresses the Drell-Yan rate by two orders of magnitude at the price of a 30% reduction 
in signal. A full GEANT-based simulation, including the real CMS geometry and dimuon reconstruction 
algorithm, is needed 12641 for more detailed conclusions. 



b bbar 

_____ D"' 




1 8G 200 



Fig. 78: The distribution of pairs from bb decays (solid histogram) and from Drell-Yan production (dashed 

histogram) as a function of 5r (see text). 
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Another process which can also carry information about medium-induced rescattering and bottom 
energy loss is secondary J /if) production 12501 12511 . The branching ratio B —>■ J/tpX is 1.15%. The 
J/t/j subsequently decays to dimuons with a 5.9% branching ratio so that gg ^ bb ^ BB X — > 

From a leading order calculation and assuming > 5 GeV/c and \r]^\ < 2.4, we expect « 13000 
dimuons from secondary J/ip decays in a one month Pb+Pb run. Primary J/t/j's produced at the nuclear 
interaction point can be rejected using tracker information on the secondary vertex position. 

Finally, we mention some of the theoretical uncertainties in the QQ cross sections and the cor- 
responding dimuon spectra. Nuclear shadowing, implemented with EKS98 fT6llT7l . is not large in this 
kinematic region, an ^ 15% effect on bb and cc decays and « 25% for Drell-Yan production. In this 
case, the predicted factor of ~ 2 — 4 66 suppression due to final state energy loss can clearly be observed 
over initial-state shadowing. However, shadowing has an 30% effect on secondary J/ip production 
from B decays. 

9.7 Conclusions 

With its w 47r muon acceptance and calorimetric coverage, CMS can make very significant and, in some 
respects, unique contribution to heavy ion physics. Studies of the T family, from pp through Ar-i-Ar and 
Pb-i-Pb, as well as from peripheral to central collisions, is likely to be of great interest at the LHC, just as 
the J ftp has been for the SPS. 

The key issue for CMS is the muon reconstruction efficiency in the tracker with the extreme 
occupancies expected in Pb-i-Pb collisions. The all-silicon tracker results in a 76% dimuon reconstruction 
efficiency in \rj\ < 0.8 for dNci^/drj = 8000. 

The pixel layers in the CMS tracker play a very important role in muon filtering and reconstruction 
efficiency, especially for the T purity, by suppressing tt/K decay muons in combination with a beam- 
line vertex constraint. Furthermore, since the Si tracker occupancies could be somewhat lower, including 
deeper Si layers in the muon reconstruction algorithm than the outermost four could improve the muon 
reconstruction efficiency in the region 0.8 < |r/| < 1.3. 

With a dimuon mass resolution of w 50 MeV, CMS is extremely well suited for quarkonia de- 
tection. The large rapidity aperture of the muon detector, as well as the precise tracking, result in high 
statistics and a very good separation between the T states. Significances between 70 for Pb-i-Pb to near 
1000 for Ar-i-Ar are found for the T. Since pj, > 3.5 GeV/c in the barrel, detection of the higher mass 
T states is favored relative to the J/tp where only J/tp's with pT > 5 GeV/c will pass the muon px 
cut. The threshold decreases from 3.5 to 2.0 GeV/c when 1.3 < rj < 2.4 so that muon detection in the 
endcaps would allow CMS to study J/ip — > over the full px range. 

Above the T mass, detection is a unique feature of CMS. The independent observation of the 
in the outer muon system and in the muon-i-tracker systems will calibrate the tracker muon recon- 
struction efficiency over an extended px range. Theoretical uncertainties in the heavy flavour production 
cross sections in nucleon-nucleon collisions at LHC energies are rather large and significantly affect the 
predicted dimuon production rates before any energy loss. These uncertainties arise from the choice of 
the parton distribution functions, the heavy quark mass, next-to-leading order corrections, heavy quark 
fragmentation, etc. Thus measurements in pp or dd collisions at the same or similar energies per nucleon 
as those for heavy ions are strongly desirable to determine the baseline rate precisely. 
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10. HEAVY FLAVOUR STUDIES IN ATLAS 



The ATLAS detector is designed for high-p-r physics in proton-proton collisions at high luminosity. The 
detector has a wide rapidity coverage and enormous potential for the study of heavy ion collisions. The 
calorimeter is highly segmented both longitudinally and transversely and offers a unique opportunity 
to detect and measure jet properties. The standalone muon system behind the calorimeter offers the 
possibility of detecting dimuons in a quiet environment. The ATLAS detector can be used to study the T 
states and 6-jets. We briefly describe the detector and discuss possible ways heavy flavour physics could 
be studied with ATLAS. 

10,1 The ATLAS Detector 

The ATLAS detector is designed to study proton-proton collisions at the LHC design energy of 14 TeV in 
the center of mass. The physics pursued by the collaboration is vast and includes the Higgs boson search, 
searches for SUSY and other scenarios beyond the Standard Model as well as precision measurements of 
process within (and possibly beyond) the Standard Model. To achieve these goals at full machine lumi- 
nosity of 10^^ cm~^s^^, ATLAS will have a precise tracking system (Inner Detector) for charged particle 
measurements, an as hermetic as possible calorimeter system with extremely fine grained segmentation, 
and a standalone muon system. An overview of the detector is shown in Fig. 110.11 

The Inner Detector is composed of a finely segmented silicon pixel detector, silicon strip detec- 
tors. Semiconductor Tracker (SCT), and the Transition Radiation Tracker (TRT). The segmentation is 
optimized for proton-proton collisions at machine design luminosity. The Inner Detector is designed to 
cover the pseudorapidity interval | r] |< 2.5 and is located inside a 2 T solenoid magnet. 

The ATLAS calorimeter system surrounding the solenoid magnet is divided into electromagnetic 
and hadronic sections and covering | r/ |< 4.9. The EM calorimeter is an accordion liquid argon device 
and is finely segmented longitudinally and transversely for | 77 | < 3.1. The first longitudinal segmentation 
has a granularity of Arj x = 0.003 x 0.1 in the barrel, becoming slightly coarser in the endcaps. 
The second longitudinal segmentation is composed of 0.025 x 0.025 cells while the last segment has 
0.05 X 0.05 cells. In addition, a finely segmented, 0.025 x 0.1, pre-sampler system is placed in front of 
the electromagnetic (EM) calorimeter. The overall energy resolution of the EM calorimeter determined 
experimentally is 10%/\/^ © 0.5%. The calorimeter also has good pointing resolution, 60mrad/\/^ 
for photons and a timing resolution of better than 200 ps for showers with energy greater than 20 GeV. 

The hadronic calorimeter is also segmented longitudinally and transversely. Except for the endcaps 
and the forward calorimeters, the calorimeter is a lead-scintillator tile structure with a granularity of 
0.1 X 0.1. In the endcaps, the hadronic calorimeter uses liquid argon technology for radiation hardness 
and has the same granularity as the barrel hadronic calorimeter. The energy resolution for the hadronic 
calorimeters is 50%/\/^ © 2% for pions. The very forward region, up to r/ = 4.9, is covered by the 
Forward Calorimeter, an axial drift liquid argon calorimeter. The overall performance of the calorimeter 
system is described in Ref. f2651. 

The ATLAS muon spectrometer is located behind the calorimeters, thus shielded from hadronic 
showers. The spectrometer uses several tracking technologies along with a toroidal magnet, providing 
a 4 T field for an independent momentum measurement outside the calorimeter volume. Most of the 
spectrometer volume is covered by Monitored Drift Tubes (MDT). In the forward region, where the rate 
is high. Cathode Strip Chambers are used. The standalone muon spectrometer momentum resolution is 
of the order of 2% for muons with 10 < pr < 100 GeV. The muon spectrometer covers | r] |< 2.7. 

The ATLAS trigger and data acquisition system is multi-level and has to reduce the beam crossing 
rate of 40 MHz to an output to mass storage rate of 0(100) Hz. The first stage (LVLl) is a hardware- 
based trigger, making use of the coarse granularity calorimeter data and the dedicated muon trigger 
chambers only, reducing the output rate to about 75 kHz, within a maximum latency of 2.5 fis. 

^'Authors: B. Cole, H. Takai, S. Tapprogge. 
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Fig. 79: The overall layout of the ATLAS detector. 
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The performance results have been obtained using a detailed full simulation of the ATLAS detector 
response with GEANT and have been validated by an extensive program of test beam measurements of 
all components. 

10.2 Heavy Flavour studies in ATLAS 

Recent theoretical investigations fl4T\ have indicated that a dense partonic medium will suppress the 
gluon radiation from charm and bottom quarks relative to that from light quarks. Thus 6-jet measure- 
ments in ATLAS would provide an important comparison measurement to light quark and gluon jets. 
Such measurements would have important implications for gluon shadowing and saturation effects. Tag- 
ging of 5-jets by the associated muon is possible in the pp environment 12651 . 

It is unlikely that vertexing will be available for these studies due to the high multiplicity Pb-i-Pb 
environment. Therefore we are currently studying the possibility of tagging 6-jets by matching a mea- 
sured muon in the muon spectrometer to the jet measured by the calorimeter. Vertexing will be available 
for lighter ion and pA collisions. 

We will profit from current studies of ATLAS performance for soft lepton tagging of the H —>^ bb 
decay. Initial studies of the detector performance for fe-jet tagging using this technique can be found in 
the TDR for high luminosity pp collisions \266\. We note that the muon momentum should be larger 
than 4 GeV to be detected in the spectrometer. The expected muon background comes from vr, K and c 
decays. 

Quarkonia suppression is expected in a deconfined medium due to the screening of the long-range 
attractive potential. We are beginning studies of the ATLAS capability for T identification. The initial 
evaluation is that the muon system alone will provide marginal resolution for a clear separation of the 
three states. However, use of the SCT and pixel detectors can enhance the mass resolution. Initial full 
simulations of Pb-i-Pb collisions indicate that the occupancies in the Inner Detector (with the exception 
of the Transition Radiation Tracker) are low enough to allow track reconstruction. We expect that all of 
the inner detectors will be available for lighter ion and pA collisions. The expected occupancies obtained 
from preliminary full detector simulations using HIJING are listed in Table l35]for the barrel detectors in 
the Inner Detector. In the endcaps, the Inner Detector covers 1.4 < rj < 2.5 with average occupancies of 
0.2% for the silicon pixels and 6% and 10% for the SCT, depending on the SCT wheel. There are 9 SCT 
tracking stations on each endcap which should allow tracking in the forward direction as well. 

Table 35: Segmentation and occupancies of the Inner Detector (Barrel). Values were obtained from preliminary full detector 
simulations using HIJING with dNch/drj = 8000. The pp occupancies are from Ref. 12671 . The corresponding pp and Pb+Pb 
luminosities are also quoted. 



ID 


r] range 


Radius (cm) 


Channels 


Occupancy 


Occupancy 










pp (lO-'^'^) 


Fh+Ph (1027) 


Si Pixel 1 


±2.5 


5.0 


10'^ 


0.04 


0.6 


Si Pixel 2 


±1.7 


8.8 


8 X lO'^ 


0.01 


0.3 


Si Pixel 3 


±1.7 


12.2 


8 X lO'^ 


0.006 


0.2 


SCT Layer 1 


±1.4 


30.0 


3 X 10^ 


0.6 


15 


SCT Layer 1 


±1.4 


37.0 


3 X 10^ 


0.5 


10 


SCT Layer 1 


±1.4 


44.0 


3 X 10^ 


0.4 


7.5 


SCT Layer 1 


±1.4 


51.5 


3 X 10^ 


0.3 


6.0 



The muon spectrometer will benefit from the low nucleus-nucleus luminosity. In the high lumi- 
nosity pp runs a large number of spectrometer hits come from slow neutrons produced in previous inter- 
actions. Initial simulations using HIJING Pb-i-Pb events indicate that the muon spectrometer is sparsely 
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populated. The muon backgrounds are from tt and K decays in the calorimeter. However, a substantial 
number of these decays are expected to be rejected when matches to Timer Detector tracks are required. 
The minimum muon pr for detection is ~ 5 — 6 GeV. 
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